dolfin team mailing list archive
-
dolfin team
-
Mailing list archive
-
Message #18632
Re: [Bug 594954] Re: uBLAS GMRES solver fails on first run
Is there any specific reason why we have this type of problems in
resizing uBLAS vectors not other type of vectors?
Note that resizing a vector becomes important in the adaptive methods.
Since when a function updates, we need to resize its underlying vector.
If we clear a vector, then we lose all data that we had there.
Mehdi
On Mon, 2010-06-28 at 15:01 +0000, Anders Logg wrote:
> The problem was that x = Vector(b.size()) created a vector which
> sometimes contained strange numbers. This has been changed now (in
> uBLASVector::resize()) so that the vector is always zero. Took me a
> while to work through all the 29 different axpy_prod that call each
> other in operation.hpp in uBLAS...
>
> Should be fixed now.
>
> ** Changed in: dolfin
> Status: Confirmed => Fix Committed
>
Follow ups
References