<div dir="ltr"><div><div>Hi Jose,<br></div>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.<br></div>Amlan <br>
</div><div class="gmail_extra"><br><br><div class="gmail_quote">On Sun, Jan 19, 2014 at 12:12 AM, Jose E. Roman <span dir="ltr"><<a href="mailto:jroman@dsic.upv.es" target="_blank">jroman@dsic.upv.es</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
El 18/01/2014, a las 18:43, amlan barua escribió:<br>
<div class="im"><br>
> 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.<br>
> Amlan<br>
><br>
<br>
</div><div class="im">Please do Reply-all so that the answer goes to the list.<br>
<br>
</div>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<br>
Use VecView(y,NULL) to see the result.<br>
<span class="HOEnZb"><font color="#888888"><br>
Jose<br>
<br>
</font></span></blockquote></div><br></div>