On Mon, Jul 16, 2012 at 3:23 AM, Shitij Bhargava <span dir="ltr"><<a href="mailto:shitij.cse@gmail.com" target="_blank">shitij.cse@gmail.com</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi all !!<br><br>I am using Fortan and petsc 3.2-p7.<div><br>I am creating the array as:<br><div><i><font color="#000066">call MatCreateSeqAIJ(PETSC_COMM_SELF,natoms,natoms,PETSC_DECIDE,PETSC_NULL,bonds,ierr)</font></i></div>

</div><div><br></div><div>And then inserting elements as:<br><font color="#000066"><i>call MatSetValue(bonds,temp_int1-1,temp_int2-1,1,INSERT_VALUES,ierr)</i></font><br><br></div><div>I even did:<br><font color="#000066"><i>call MatSetValue(bonds,1,8,35,INSERT_VALUES,ierr)</i></font><br>

<br>then I finalize the matrix by these:</div><div><div><font color="#000066"><i>call MatAssemblyBegin(bonds,MAT_FINAL_ASSEMBLY,ierr)</i></font></div><div><font color="#000066"><i>call MatAssemblyEnd(bonds,MAT_FINAL_ASSEMBLY,ierr)</i></font></div>

<br>But none of these values show up when I view the matrix by this:<br><font color="#000066"><i>call MatView(bonds,PETSC_VIEWER_STDOUT_SELF,ierr)</i></font><br><br>It's output is the following:<br><div><font color="#660000">Matrix Object: 1 MPI processes</font></div>

<div><font color="#660000">  type: seqaij</font></div><div><font color="#660000">row 0: (0, 2.47033e-323)  (1, 2.8034e-85)  (4, 2.8034e-85)  (5, 2.8034e-85)  (6, 2.8034e-85) </font></div><div><font color="#660000">row 1: (0, 2.8034e-85)  (2, 2.8034e-85)  (7, 2.8034e-85)  (8, 1.6976e-313) </font></div>

<div><font color="#660000">row 2: (1, 2.8034e-85)  (3, 2.8034e-85)  (9, 2.8034e-85) </font></div><div><font color="#660000">row 3: (2, 2.8034e-85) </font></div><div><font color="#660000">row 4: (0, 2.8034e-85) </font></div>

<div><font color="#660000">row 5: (0, 2.8034e-85) </font></div><div><font color="#660000">row 6: (0, 2.8034e-85) </font></div><div><font color="#660000">row 7: (1, 2.8034e-85) </font></div><div><font color="#660000">row 8: (1, 2.8034e-85) </font></div>

</div><div><font color="#660000">..............</font></div><div><font color="#660000"><br></font></div><div>Why are the values not being inserted? (Even at (1,8), 35 is not there.)</div></blockquote><div><br></div><div>
Because the values are the wrong type, and Fortran does not know how to coerce them</div><div>correctly (typing is weak in Fortran). Take a look at</div><div><br></div><div>  <a href="http://www.mcs.anl.gov/petsc/petsc-current/src/ksp/ksp/examples/tutorials/ex2f.F.html">http://www.mcs.anl.gov/petsc/petsc-current/src/ksp/ksp/examples/tutorials/ex2f.F.html</a></div>
<div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div>Thank you in advance !!<span class="HOEnZb"><font color="#888888"><br>

<br>Shitij</font></span></div><div><br></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>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<br>