Skip to content
This repository has been archived by the owner on Mar 8, 2021. It is now read-only.

Commit

Permalink
Browse files Browse the repository at this point in the history
Use a more conservative ARPACK workspace estimate.
  • Loading branch information
inducer committed Nov 12, 2009
1 parent 7d388d8 commit d214410
Showing 1 changed file with 2 additions and 6 deletions.
8 changes: 2 additions & 6 deletions src/cpp/pyublasext/arpack.hpp
Expand Up @@ -268,12 +268,8 @@ namespace boost { namespace numeric { namespace bindings { namespace arpack {

boost::scoped_array<value_type> workd(new value_type[3*n]);
int lworkl;
if (detail::is_complex(value_type()))
lworkl = 3 * number_of_arnoldi_vectors * number_of_arnoldi_vectors
+ 5 * number_of_arnoldi_vectors;
else
lworkl = 3 * number_of_arnoldi_vectors * number_of_arnoldi_vectors
+ 6 * number_of_arnoldi_vectors;
lworkl = 3 * number_of_arnoldi_vectors * number_of_arnoldi_vectors
+ 6 * number_of_arnoldi_vectors;

boost::scoped_array<value_type> workl(new value_type[lworkl]);
boost::scoped_array<double> rwork(new double[number_of_arnoldi_vectors]);
Expand Down

0 comments on commit d214410

Please sign in to comment.