The Fortran interface fell behind. I will fix it in dev, but for now you can<br>
just use the value, 91. You can also use the command line option.<br>
<br>
Thanks,<br>
<br>
Matt<br><br><div><span class="gmail_quote">On 5/22/06, <b class="gmail_sendername">Randall Mackie</b> <<a href="mailto:randy@geosystem.us">randy@geosystem.us</a>> wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
I'm trying to modify my code to take advantage of the fact that my system matrix<br>is symmetric.<br><br>It seemed from reading the manual pages that I could call MatCreateMPISBAIJ<br>to create a symmetric data matrix.<br>
<br>However, when I tried this, I got the error:<br><br>[1]PETSC ERROR: Lower triangular value cannot be set for sbaij format. Ignoring these value<br>s, run with -mat_ignore_lower_triangular or call MatSetOption(mat,MAT_IGNORE_LOWER_TRIANGUL
<br>AR)!<br><br>So I tried including the line:<br><br> call MatSetOption(A,MAT_IGNORE_LOWER_TRIANGULAR,ierr)<br><br>But when I tried to compile the code, I then got this message:<br><br>fortcom:
Error: csem3dfwd.F, line 771: This name does not have a type, and must
have an explicit type. [MAT_IGNORE_LOWER_TRIANGULAR]<br> call MatSetOption(A,MAT_IGNORE_LOWER_TRIANGULAR,ierr)<br>--------------------------^<br><br><br>So, can anyone point me in the right direction for taking advantage of the fact that
<br>my matrices are symmetric? (they are actually complex symmetric, and NOT Hermitian).<br><br>Thanks, Randy<br><br>--<br>Randall Mackie<br>GSY-USA, Inc.<br>PMB# 643<br>2261 Market St.,<br>San Francisco, CA 94114-1600<br>
Tel (415) 469-8649<br>Fax (415) 469-5044<br><br>California Registered Geophysicist<br>License No. GP 1034<br><br></blockquote></div><br><br clear="all"><br>-- <br>"Failure has a thousand explanations. Success doesn't need one" -- Sir Alec Guiness