<div dir="ltr"><div dir="ltr">On Mon, Nov 14, 2022 at 7:39 PM Mohammad Ali Yaqteen <<a href="mailto:mhyaqteen@sju.ac.kr">mhyaqteen@sju.ac.kr</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class="msg-4839079817109792065">




<div dir="ltr">
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)">
I am using the following procedure to read from Vec, but it keeps giving me the same values! I was told that using VecGetValues gives wrog output. </div></div></div></blockquote><div><br></div><div>By who? It does not give the wrong output.</div><div><br></div><div>You do not show where in the code you define tdof and ei[].</div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class="msg-4839079817109792065"><div dir="ltr"><div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)">If not this, then what function should be used to read the contents of a vector?</div></div></div></blockquote><div><br></div><div>I think it would be simpler for you to use VecGetArrayRead(), unless you want values from other processes.</div><div><br></div><div>  Thanks,</div><div><br></div><div>    Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class="msg-4839079817109792065"><div dir="ltr">
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)">
for (int i = 0; i < nconv; i++)
<div>        {</div>
<div>            PetscCall(EPSGetEigenpair(eps,i,&kr,&ki,xr,xi));</div>
<div>            PetscCall(EPSComputeError(eps,i,EPS_ERROR_RELATIVE,&error1));</div>
<div><br>
</div>
<div>            #if defined(PETSC_USE_COMPLEX)</div>
<div>                re = PetscRealPart(kr);</div>
<div>                im = PetscImaginaryPart(kr);</div>
<div>            #else</div>
<div>                re = kr;</div>
<div>                im = ki;</div>
<div>            #endif</div>
<div>            if (im!=0.0) PetscCall(PetscPrintf(PETSC_COMM_WORLD," %9f%+9fi %12g\n",(double)re,(double)im,(double)error1));</div>
<div>            else PetscCall(PetscPrintf(PETSC_COMM_WORLD,"   %12f       %12g\n",(double)re,(double)error1));</div>
<div>            eval(i) = re;</div>
<div>            VecGetValues(xr, tdof, ei, eveci);</div>
<div>            for (int j = 0; j < tdof; j++)</div>
<div>            {</div>
<div>                evec(j, i) = eveci[j];</div>
<div>            }</div>
<div>        }</div>
        PetscCall(PetscPrintf(PETSC_COMM_WORLD,"\n"));<br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)">
<br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)">
Thank you</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0);background-color:rgb(255,255,255)">
Ali</div>
</div>

</div></blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>