[petsc-users] computing matrix exponentials

[amlan barua]

Hi Jose,
Thanks for the clarification. So if I want to compute the exponential
matrix just once then I should extract the columns one by one and store
them in a PetSc matrix.
Amlan


On Sun, Jan 19, 2014 at 12:12 AM, Jose E. Roman <jroman at dsic.upv.es> wrote:

>
> El 18/01/2014, a las 18:43, amlan barua escribió:
>
> > Also I tried to get the first column of exp(A) by multiplying it with
> e_1 = (1,0,0,\cdots). The code crashes. I have attached the code.
> > Amlan
> >
>
> Please do Reply-all so that the answer goes to the list.
>
> There is a problem when ncv equals the matrix size (will add a fix when I
> have time). Instead of n=10, try a larger matrix. Or alternatively run with
> -mfn_ncv 9
> Use VecView(y,NULL) to see the result.
>
> Jose
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140119/5dbba196/attachment.html>