[vsipl++] [patch] double support for SAL LU solver

[Jules Bergmann]

Don, Thanks for getting this checked in! -- Jules


Don McCoy wrote:
>   VSIP_THROW((std::bad_alloc))
>     : length_ (lu.length_),
>       ipiv_   (length_),
> +     recip_  (length_),
>       data_   (length_, length_)
>   {
>     data_ = lu.data_;
>     for (index_type i=0; i<length_; ++i)
> +   {
>       ipiv_[i] = lu.ipiv_[i];
> +     recip_ = lu.recip_;

The recip_ assignment should go outside the loop, right?

> +   }
>   }


-- 
Jules Bergmann
CodeSourcery
jules at codesourcery.com
(650) 331-3385 x705