<div style="line-height:1.7;color:#000000;font-size:14px;font-family:arial"><DIV>Hello, Matthew </DIV>
<DIV> Do you mean that the two versions have a different in this point? If I use the new version, I have to</DIV>
<DIV>make some modifications on my codes?</DIV>
<DIV> Thanks.</DIV>
<DIV> Jim<BR><BR><BR><BR><BR></DIV>
<DIV></DIV>
<DIV id="divNeteaseMailCard"></DIV>
<DIV><BR></DIV>>ÔÚ 2012-11-01 21:49:35£¬"Matthew Knepley" <knepley@gmail.com> дµÀ£º<BR>
<BLOCKQUOTE style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" id="isReplyContent">>On Thu, Nov 1, 2012 at 9:47 AM, w_ang_temp <SPAN dir="ltr"><<A href="mailto:w_ang_temp@163.com" target="_blank">w_ang_temp@163.com</A>></SPAN> wrote:<BR>
<DIV class="gmail_quote">
<BLOCKQUOTE style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" class="gmail_quote">
<DIV style="LINE-HEIGHT: 1.7; FONT-FAMILY: arial; FONT-SIZE: 14px">
<DIV>>Hello,</DIV>
<DIV>> I have just used the latest version 'petsc-3.3-p4.tar.gz'. The codes which are ok under version 'petsc-3.2-p7.tar.gz'</DIV>
<DIV>>are not ok now. The error infomation is as follows. So if can I still use the original codes without any</DIV>
<DIV>>modification and also use multigrid with -pc_type gamg -pc_gamg_agg_nsmooths 1.</DIV>
<DIV>> Thanks.</DIV>
<DIV> Jim</DIV>
<DIV>>[0]PETSC ERROR: --------------------- Error Message ------------------------------------<BR>>[0]PETSC ERROR: Object is in wrong state!<BR>>[0]PETSC ERROR: Must call MatXXXSetPreallocation() or MatSetUp() on argument 1 "mat" before MatGetOwnershipRange()!<BR></DIV></DIV></BLOCKQUOTE>
<DIV><BR></DIV>
<DIV>>This says what is wrong. Now you must either preallocate your matrix or turn of this error using MatSetOption()</DIV>
<DIV><BR></DIV>
<DIV> > Matt</DIV>
<DIV> </DIV>
<BLOCKQUOTE style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" class="gmail_quote">
<DIV style="LINE-HEIGHT: 1.7; FONT-FAMILY: arial; FONT-SIZE: 14px">
<DIV>>>[0]PETSC ERROR: ------------------------------------------------------------------------<BR>>>[0]PETSC ERROR: Petsc Release Version 3.3.0, Patch 4, Fri Oct 26 10:46:51 CDT 2012 <BR>>>[0]PETSC ERROR: See docs/changes/index.html for recent updates.<BR>>>[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.<BR>>>[0]PETSC ERROR: See docs/index.html for manual pages.<BR>>>[0]PETSC ERROR: ------------------------------------------------------------------------<BR>>>[0]PETSC ERROR: ./ex4f on a arch-linu named ubuntu by ubu Thu Nov 1 06:10:23 2012<BR>>>[0]PETSC ERROR: Libraries linked from /home/ubu/soft/petsc/petsc-3.3-p4/arch-linux2-c-debug/lib<BR>>>[0]PETSC ERROR: Configure run at Thu Nov 1 05:54:48 2012<BR>>>[0]PETSC ERROR: Configure options --with-mpi-dir=/home/ubu/soft/mpich2/ --download-f-blas-lapack =1<BR>>>[0]PETSC ERROR: ------------------------------------------------------------------------<BR>>>[0]PETSC ERROR: MatGetOwnershipRange() line 5992 in src/mat/interface/matrix.c<BR>>>[1]PETSC ERROR: [2]PETSC ERROR: --------------------- Error Message ------------------------------------<BR>>>[2]PETSC ERROR: Object is in wrong state!<BR>>>[2]PETSC ERROR: Must call MatXXXSetPreallocation() or MatSetUp() on argument 1 "mat" before MatGetOwnershipRange()!<BR>>>[2]PETSC ERROR: ------------------------------------------------------------------------<BR>>>[2]PETSC ERROR: Petsc Release Version 3.3.0, Patch 4, Fri Oct 26 10:46:51 CDT 2012 <BR>>>[2]PETSC ERROR: See docs/changes/index.html for recent updates.<BR>>>[2]PETSC ERROR: See docs/faq.html for hints about trouble shooting.<BR>>>[2]PETSC ERROR: See docs/index.html for manual pages.<BR>>>[2]PETSC ERROR: ------------------------------------------------------------------------<BR>>>[2]PETSC ERROR: ./ex4f on a arch-linu named ubuntu by ubu Thu Nov 1 06:10:23 2012<BR>>>[2]PETSC ERROR: Libraries linked from /home/ubu/soft/petsc/petsc-3.3-p4/arch-linux2-c-debug/lib<BR>>>[2]PETSC ERROR: Configure run at Thu Nov 1 05:54:48 2012<BR>>>[2]PETSC ERROR: Configure options --with-mpi-dir=/home/ubu/soft/mpich2/ --download-f-blas-lapack =1<BR>>>[2]PETSC ERROR: ------------------------------------------------------------------------<BR>>>[2]PETSC ERROR: MatGetOwnershipRange() line 5992 in src/mat/interface/matrix.c<BR>>>[3]PETSC ERROR: --------------------- Error Message ------------------------------------<BR>>>[3]PETSC ERROR: Object is in wrong state!<BR>>>[3]PETSC ERROR: Must call MatXXXSetPreallocation() or MatSetUp() on argument 1 "mat" before MatGetOwnershipRange()!<BR>>>[3]PETSC ERROR: ------------------------------------------------------------------------<BR>>>[3]PETSC ERROR: Petsc Release Version 3.3.0, Patch 4, Fri Oct 26 10:46:51 CDT 2012 <BR>>>[3]PETSC ERROR: See docs/changes/index.html for recent updates.<BR>>>[3]PETSC ERROR: See docs/faq.html for hints about trouble shooting.<BR>>>[3]PETSC ERROR: See docs/index.html for manual pages.<BR>>>[3]PETSC ERROR: ------------------------------------------------------------------------<BR>>>[3]PETSC ERROR: ./ex4f on a arch-linu named ubuntu by ubu Thu Nov 1 06:10:23 2012<BR>>>[3]PETSC ERROR: Libraries linked from /home/ubu/soft/petsc/petsc-3.3-p4/arch-linux2-c-debug/lib<BR>>>[3]PETSC ERROR: Configure run at Thu Nov 1 05:54:48 2012<BR>>>[3]PETSC ERROR: Configure options --with-mpi-dir=/home/ubu/soft/mpich2/ --download-f-blas-lapack =1<BR>>>[3]PETSC ERROR: ------------------------------------------------------------------------<BR>>>[3]PETSC ERROR: MatGetOwnershipRange() line 5992 in src/mat/interface/matrix.c<BR>>>--------------------- Error Message ------------------------------------<BR>>>[1]PETSC ERROR: Object is in wrong state!<BR>>>[1]PETSC ERROR: Must call MatXXXSetPreallocation() or MatSetUp() on argument 1 "mat" before MatGetOwnershipRange()!<BR>>>[1]PETSC ERROR: ------------------------------------------------------------------------<BR>>>[1]PETSC ERROR: Petsc Release Version 3.3.0, Patch 4, Fri Oct 26 10:46:51 CDT 2012 <BR>>>[1]PETSC ERROR: See docs/changes/index.html for recent updates.<BR>>>[1]PETSC ERROR: See docs/faq.html for hints about trouble shooting.<BR>>>[1]PETSC ERROR: See docs/index.html for manual pages.<BR>>>[1]PETSC ERROR: ------------------------------------------------------------------------<BR>>>[1]PETSC ERROR: ./ex4f on a arch-linu named ubuntu by ubu Thu Nov 1 06:10:23 2012<BR>>>[1]PETSC ERROR: Libraries linked from /home/ubu/soft/petsc/petsc-3.3-p4/arch-linux2-c-debug/lib<BR>>>[1]PETSC ERROR: Configure run at Thu Nov 1 05:54:48 2012<BR>>>[1]PETSC ERROR: Configure options --with-mpi-dir=/home/ubu/soft/mpich2/ --download-f-blas-lapack =1<BR>>>[1]PETSC ERROR: ------------------------------------------------------------------------<BR>>>[1]PETSC ERROR: MatGetOwnershipRange() line 5992 in src/mat/interface/matrix.c<BR>>>[1]PETSC ERROR: [2]PETSC ERROR: [3]PETSC ERROR: --------------------- Error Message ------------------------------------<BR>>>[3]PETSC ERROR: Object is in wrong state!<BR>>>[3]PETSC ERROR: Matrix is missing diagonal entry 0!<BR>>>[3]PETSC ERROR: ------------------------------------------------------------------------<BR>>>[3]PETSC ERROR: Petsc Release Version 3.3.0, Patch 4, Fri Oct 26 10:46:51 CDT 2012 <BR>>>[3]PETSC ERROR: See docs/changes/index.html for recent updates.<BR>>>[3]PETSC ERROR: See docs/faq.html for hints about trouble shooting.<BR>>>[3]PETSC ERROR: See docs/index.html for manual pages.<BR>>>[3]PETSC ERROR: ------------------------------------------------------------------------<BR>>>[3]PETSC ERROR: ./ex4f on a arch-linu named ubuntu by ubu Thu Nov 1 06:10:23 2012<BR>>>[3]PETSC ERROR: Libraries linked from /home/ubu/soft/petsc/petsc-3.3-p4/arch-linux2-c-debug/lib<BR>>>[3]PETSC ERROR: Configure run at Thu Nov 1 05:54:48 2012<BR>>>[3]PETSC ERROR: Configure options --with-mpi-dir=/home/ubu/soft/mpich2/ --download-f-blas-lapack =1<BR>>>[3]PETSC ERROR: ------------------------------------------------------------------------<BR>>>[3]PETSC ERROR: MatILUFactorSymbolic_SeqAIJ_ilu0() line 1641 in src/mat/impls/aij/seq/aijfact.c<BR>>>[3]PETSC ERROR: MatILUFactorSymbolic_SeqAIJ() line 1745 in src/mat/impls/aij/seq/aijfact.c<BR>>>[3]PETSC ERROR: MatILUFactorSymbolic() line 6130 in src/mat/interface/matrix.c<BR>>>[3]PETSC ERROR: PCSetUp_ILU() line 216 in src/ksp/pc/impls/factor/ilu/ilu.c<BR>>>[3]PETSC ERROR: PCSetUp() line 832 in src/ksp/pc/interface/precon.c<BR>>>[3]PETSC ERROR: KSPSetUp() line 278 in src/ksp/ksp/interface/itfunc.c<BR>>>[3]PETSC ERROR: PCSetUpOnBlocks_BJacobi_Singleblock() line 715 in src/ksp/pc/impls/bjacobi/bjacobi.c<BR>>>[3]PETSC ERROR: PCSetUpOnBlocks() line 865 in src/ksp/pc/interface/precon.c<BR>>>[3]PETSC ERROR: KSPSetUpOnBlocks() line 154 in src/ksp/ksp/interface/itfunc.c<BR>>>[3]PETSC ERROR: KSPSolve() line 403 in src/ksp/ksp/interface/itfunc.c<BR>>>Fatal error in MPI_Send: Invalid count, error stack:<BR>>>MPI_Send(173): MPI_Send(buf=0x100, count=-199040697, MPI_REAL, dest=0, tag=3, MPI_COMM_WORLD) failed<BR>>>MPI_Send(97).: Negative count, value is -199040697</DIV>
<DIV> </DIV>
<DIV><BR><BR><BR><BR><BR> </DIV>
<DIV></DIV>
<DIV></DIV>
<DIV><BR></DIV>>ÔÚ 2012-10-29 21:57:49£¬"Mark F. Adams" <<A href="mailto:mark.adams@columbia.edu" target="_blank">mark.adams@columbia.edu</A>> дµÀ£º<BR>
<BLOCKQUOTE style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex">>You need an updated PETSc. (I thought 3.2 had an early version of gamg ¡ but you need 3.3 or dev)
<DIV><BR></DIV>
<DIV>>Mark</DIV>
<DIV><BR>
<DIV>
<DIV>>On Oct 29, 2012, at 8:49 AM, w_ang_temp <<A href="mailto:w_ang_temp@163.com" target="_blank">w_ang_temp@163.com</A>> wrote:</DIV><BR>
<BLOCKQUOTE type="cite">
<DIV style="LINE-HEIGHT: 1.7; FONT-FAMILY: arial; FONT-SIZE: 14px">
<DIV>>Hello, Jed</DIV>
<DIV>> I use the command:</DIV>
<DIV>> mpiexec -n 4 ./ex4f -ksp_type cgs -pc_type gamg -pc_gamg_agg_nsmooths 1 -ksp_rtol 1.0e-15 -ksp_converged_reason</DIV>
<DIV>> The error is as follows:<BR>>[0]PETSC ERROR: --------------------- Error Message ------------------------------------<BR>>[0]PETSC ERROR: Unknown type. Check for miss-spelling or missing external package needed for type<BR>> <A>seehttp://www.mcs.anl.gov/petsc/petsc-as/documentation/installation.html#external!</A><BR>>[0]PETSC ERROR: Unable to find requested PC type gamg!<BR>>[0]PETSC ERROR: ------------------------------------------------------------------------<BR>>[0]PETSC ERROR: Petsc Release Version 3.2.0, Patch 7, Thu Mar 15 09:30:51 CDT 2012 <BR>>[0]PETSC ERROR: See docs/changes/index.html for recent updates.<BR>>[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.<BR>>[0]PETSC ERROR: See docs/index.html for manual pages.<BR>>[0]PETSC ERROR: ------------------------------------------------------------------------<BR>>[0]PETSC ERROR: ./ex4f on a arch-linu named ubuntu by geo Mon Oct 29 05:40:13 2012<BR>>[0]PETSC ERROR: Libraries linked from /home /geo/soft/petsc/petsc-3.2-p7/arch-linux2-c-opt/lib<BR>>[0]PETSC ERROR: Configure run at Mon Jul 2 20:33:17 2012<BR>>[0]PETSC ERROR: Configure options --with-mpi-dir=/home/geo/soft/mpich2 --download-f-blas-lapack=1 --with-x=1 --with-debugging=0 -->download-parmetis --download-mumps --download-scalapack --download-blacs<BR>>[0]PETSC ERROR: ------------------------------------------------------------------------<BR>>[0]PETSC ERROR: PCSetType() line 67 in src/ksp/pc/interface/pcset.c<BR>>[0]PETSC ERROR: PCSetFromOptions() line 184 in src/ksp/pc/interface/pcset.c<BR>>[0]PETSC ERROR: KSPSetFromOptions() line 286 in src/ksp/ksp/interface/itcl.c</DIV>
<DIV>> When I use 'mpiexec -n 4 ./ex4f -ksp_type cgs -pc_type sor -ksp_rtol 1.0e-15 -ksp_converged_reason', it is ok. So</DIV>
<DIV>>what is the possible reason?</DIV>
<DIV>> Thanks.</DIV>
<DIV>> Jim<BR><BR><BR><BR><BR></DIV>
<DIV></DIV>
<DIV></DIV>
<DIV><BR></DIV>>>At 2012-10-28 22:25:31,"Jed Brown" <<A href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@mcs.anl.gov</A>> wrote:<BR>
<BLOCKQUOTE style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex">
<P>>>Try the simple option I just sent.</P>
<DIV class="gmail_quote">>>On Oct 28, 2012 6:38 AM, "w_ang_temp" <<A href="mailto:w_ang_temp@163.com" target="_blank">w_ang_temp@163.com</A>> wrote:<BR type="attribution">
<BLOCKQUOTE style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" class="gmail_quote">
<DIV style="LINE-HEIGHT: 1.7; FONT-FAMILY: arial; FONT-SIZE: 14px">
<DIV>>>Hello, Jed</DIV>
<DIV>>> Thanks for your timely reply. I deal with the soil-water coupled problem in geotechnical engineering, </DIV>
<DIV>>>whose stiffness matrix is ill-conditioned. I have did some work about it, mainly finding the effective</DIV>
<DIV>>>solvers and preconditioners. I used the command line option like this:</DIV>
<DIV>>>mpiexec -n 4 ./ex4f -ksp_type cgs -pc_type sor -ksp_rtol 1.0e-15 -ksp_converged_reason.</DIV>
<DIV>>> So, I also want to use the multigrid just like the simple command. There is only a little introduction</DIV>
<DIV>>>about multigrid in the manual. Multigrid is complex and not a easy thing for me, so I just need to konw how</DIV>
<DIV>>>to use it simply in PETSc to solve the Ax=b system.</DIV>
<DIV>>> Thanks.</DIV>
<DIV> Jim<BR><BR><BR><BR><BR></DIV>
<DIV></DIV>
<DIV></DIV>
<DIV><BR></DIV>>At 2012-10-28 21:17:00,"Jed Brown" <<A href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@mcs.anl.gov</A>> wrote:<BR>
<BLOCKQUOTE style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex">
<P>>Algebraic multigrid can be used directly, -pc_type gamg -pc_gamg_agg_nsmooths 1. Geometric either required that you use the PCMG interface to set >interpolation (and provide a coarse operator for non-Galerkin) or use a DM that provides coarsening capability.</P>
<P>>What kind of problem are you solving?</P>
<DIV class="gmail_quote">>>On Oct 28, 2012 6:09 AM, "w_ang_temp" <<A href="mailto:w_ang_temp@163.com" target="_blank">w_ang_temp@163.com</A>> wrote:<BR type="attribution">
<BLOCKQUOTE style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" class="gmail_quote">
<DIV style="LINE-HEIGHT: 1.7; FONT-FAMILY: arial; FONT-SIZE: 14px">
<DIV>>>Hello,</DIV>
<DIV>>> I want to use the multigrid as a preconditioner. The introduction about it in the manual is little.</DIV>
<DIV>>>So are there some typical examples or details about multigrid? Is it used just like other preconditioners</DIV>
<DIV>>>like jacobi, sor, which can be simply used in the cammand line options?</DIV>
<DIV>>> Thanks.</DIV>
<DIV>>> Jim</DIV></DIV><BR><BR><SPAN title="neteasefooter"><SPAN></SPAN></SPAN></BLOCKQUOTE></DIV></BLOCKQUOTE></DIV><BR><BR><SPAN title="neteasefooter"><SPAN></SPAN></SPAN></BLOCKQUOTE></DIV></BLOCKQUOTE></DIV><BR><BR><SPAN title="neteasefooter"><SPAN></SPAN></SPAN></BLOCKQUOTE></DIV><BR></DIV></BLOCKQUOTE></DIV><BR><BR><SPAN title="neteasefooter"><SPAN></SPAN></SPAN></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></BLOCKQUOTE></div><br><br><span title="neteasefooter"><span id="netease_mail_footer"></span></span>