[petsc-users] computing matrix exponentials

[Jose E. Roman]

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