On 5/22/06, <b class="gmail_sendername">Randall Mackie</b> <<a href="mailto:randy@geosystem.us">randy@geosystem.us</a>> wrote:<div><span class="gmail_quote"></span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi Matt,<br><br>Thanks, that fixed the compilation problem, but then the program just seemed to<br>hang at the place where it was creating/assembling the matrix.</blockquote><div><br>
My guess is that it is not hanging, but rather allocating all over the place. If you give<br>
-info, it shouldlet you know how manty mallocs were used.<br>
</div><br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Maybe I do not understand how to set bs, the block size for the call<br>to MatCreateMPISBAIJ, any advice?
</blockquote><div><br>
If you do not know the block size, just use 1 to start. <br>
</div><br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Are symmetric MPI BAIJ matrices my only choice? Is there not a symmetric MPIAIJ<br>
matrix format?</blockquote><div><br>
That is really block size 1.<br>
<br>
Matt <br>
</div><br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Thanks, Randy<br><br><br>Matthew Knepley wrote:<br>> 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>> On 5/22/06, *Randall Mackie* <<a href="mailto:randy@geosystem.us">randy@geosystem.us
</a><br>> <mailto:<a href="mailto:randy@geosystem.us">randy@geosystem.us</a>>> wrote:<br>><br>> I'm trying to modify my code to take advantage of the fact that my<br>> system matrix<br>> is symmetric.
<br>><br>> It seemed from reading the manual pages that I could call<br>> 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<br>> format. Ignoring these value<br>> s, run with -mat_ignore_lower_triangular or call<br>> 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<br>> 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<br>> of the fact that<br>> my matrices are symmetric? (they are actually complex symmetric, and
<br>> 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>><br>><br>><br>> --<br>> "Failure has a thousand explanations. Success doesn't need one" -- Sir
<br>> Alec Guiness<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