dolfin team mailing list archive

[Johan Hake <johan.hake@xxxxxxxxx>]

On Monday 22 March 2010 20:04:47 Murtazo Nazarov wrote:
> New question #105253 on DOLFIN:
> https://answers.launchpad.net/dolfin/+question/105253
> 
> I am trying to use Lu factorized solver, for that did like this.
> 
> LUSolver solver("lu");
> Matrix M;
> Vector b;
> assemble(M, *a);
> assemble(b, *L);
> solver.factorize(M);
> solver.factorized_solve(u.vector(), b);
> 
> The matrix M is a mass matrix. But I get runtime error:
> 
> terminate called after throwing an instance of 'std::runtime_error'
>   what():  *** Error: Unable to return pointers to underlying matrix data.
> 
> Did someone experienced such output? Do you know how to fix it?

I think Garth can give a more comprehensive answer, but it looks like 
solver.factorize only works for the uBLAS backend. So you probably need to 
set:

  parameters["linear_algebra_backend"] = "uBLAS";

in the beginning.

Johan