← Back to team overview

dolfin team mailing list archive

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