<div style="line-height:1.7;color:#000000;font-size:14px;font-family:arial"><DIV>(1)  using AMG<BR>[0]PCSetData_AGG bs=1 MM=7601<BR>Linear solve converged due to CONVERGED_RTOL iterations 445<BR>Norm of error  0.2591E+04 iterations   445<BR>  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000<BR> -2.105776715959587E-017  0.000000000000000E+000  0.000000000000000E+000<BR>   26.4211453778391      -3.262172452839194E-017 -2.114490133288630E-017<BR>************************************************************************************************************************<BR>***             WIDEN YOUR WINDOW TO 120 CHARACTERS.  Use 'enscript -r -fCourier9' to print this document            ***<BR>********************************************************************
 ****************************************************</DIV>
<DIV>---------------------------------------------- PETSc Performance Summary: ----------------------------------------------</DIV>
<DIV>./ex4f on a arch-linux2-c-debug named ubuntu with 4 processors, by ubu Sat Nov  3 09:29:28 2012<BR>Using Petsc Release Version 3.3.0, Patch 4, Fri Oct 26 10:46:51 CDT 2012 </DIV>
<DIV>                         Max       Max/Min        Avg      Total <BR>Time (sec):           3.198e+02      1.00002   3.198e+02<BR>Objects:              4.480e+02      1.00000   4.480e+02<BR>Flops:                2.296e+09      1.08346   2.172e+09  8.689e+09<BR>Flops/sec:            7.181e+06      1.08344   6.792e+06  2.717e+07<BR>Memory:       
         2.374e+07      1.04179              9.297e+07<BR>MPI Messages:         6.843e+03      1.87582   5.472e+03  2.189e+04<BR>MPI Message Lengths:  2.660e+07      2.08884   3.446e+03  7.542e+07<BR>MPI Reductions:       6.002e+04      1.00000</DIV>
<DIV>Flop counting convention: 1 flop = 1 real number operation of type (multiply/divide/add/subtract)<BR>                            e.g., VecAXPY() for real vectors of length N --> 2N flops<BR>                            and VecAXPY() for complex vectors of length N --> 8N flops</DIV>
<DIV>Summary of Stages:   ----- Time ------  ----- Flops -----  --- Messages ---  -- Message Lengths --  -- Reductions --<BR>                        Avg     %Total     Avg     %Total   counts   %Total     Avg         %Total   counts   %Total <BR> 0:      Main Stage: 3.1981e+02 100.0%  8.6886e+09 100.0%  2.189e+04 100.0%  3.446e+03      100.0%  6.001e+04 100.0% <BR>------------------------------------------------------------------------------------------------------------------------<BR>See the 'Profiling' chapter of the users' manual for details on interpreting output.<BR>Phase summary info:<BR> &n
 bsp; Count: number of times phase was executed<BR>   Time and Flops: Max - maximum over all processors<BR>                   Ratio - ratio of maximum to minimum over all processors<BR>   Mess: number of messages sent<BR>   Avg. len: average message length<BR>   Reduct: number of global reductions<BR>   Global: entire computation<BR>   Stage: stages of a computation. Set stages with PetscLogStagePush() and PetscLogStagePop().<BR>      %T - percent time in this phase         %f - percent flops in this phase<BR>      %M - percent messages in this phase     %L - percent message lengths in this phase<BR>      %R - percent reductions in this phase<BR>   Total Mflop/s: 10e-6 * (sum of flops o
 ver all processors)/(max time over all processors)<BR>------------------------------------------------------------------------------------------------------------------------</DIV>
<DIV><BR>      ##########################################################<BR>      #                                                        #<BR>      #                          WARNING!!!                    #<BR>      #                             &
 nbsp;                          #<BR>      #   This code was compiled with a debugging option,      #<BR>      #   To get timing results run ./configure                #<BR>      #   using --with-debugging=no, the performance will      #<BR>      #   be generally two or three times faster.              #<BR>      #                                &n
 bsp;                       #<BR>      ##########################################################</DIV>
<DIV><BR>Event                Count      Time (sec)     Flops                             --- Global ---  --- Stage ---   Total<BR>                   Max Ratio  Max     Ratio   Max  Ratio  Mess   Avg len Reduct  %T %f %M %L %R  %T %f %M %L %R Mflop/s<BR>------------------------------------------------------------------------------------------------------------------------</DIV>
<DIV>--- Event Stage 0: Main Stage</DIV>
<DIV>MatMult             6291 1.0 2.1431e+01 4.2 1.01e+09 1.2 1.9e+04 3.4e+03 0.0e+00  4 42 86 85  0   4 42 86 85  0   170<BR>MatMultAdd           896 1.0 4.8204e-01 1.1 2.79e+07 1.2 1.3e+03 3.4e+03 0.0e+00  0  1  6  6  0   0  1  6  6  0   208<BR>MatMultTranspose     896 1.0 2.2052e+00 1.3 2.79e+07 1.2 1.3e+03 3.4e+03 1.8e+03  1  1  6  6  3   1  1  6  6  3    45<BR>MatSolve             896 0.0 1.4953e-02 0.0 2.44e+06 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0   163<BR>MatLUFactorSym         1 1.0 1.7595e
 -04 4.7 0.00e+00 0.0 0.0e+00 0.0e+00 5.0e+00  0  0  0  0  0   0  0  0  0  0     0<BR>MatLUFactorNum         1 1.0 1.3995e-0423.5 1.85e+04 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0   132<BR>MatConvert             2 1.0 3.1026e-02 1.7 0.00e+00 0.0 0.0e+00 0.0e+00 1.2e+01  0  0  0  0  0   0  0  0  0  0     0<BR>MatScale               6 1.0 8.9679e-03 5.7 3.98e+05 1.2 6.0e+00 3.4e+03 4.0e+00  0  0  0  0  0   0  0  0  0  0   158<BR>MatAssemblyBegin      37 1.0 2.1544e-01 1.7 0.00e+00 0.0 5.4e+01 6
 .4e+03 4.2e+01  0  0  0  0  0   0  0  0  0  0     0<BR>MatAssemblyEnd        37 1.0 2.5336e-01 1.4 0.00e+00 0.0 9.0e+01 6.8e+02 3.1e+02  0  0  0  0  1   0  0  0  0  1     0<BR>MatGetRow          26874 1.2 9.8243e-02 1.6 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0<BR>MatGetRowIJ            1 0.0 2.5988e-05 0.0 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0<BR>MatGetOrdering         1 0.0 1.5616e-04 0.0 0.00e+00 0.0 0.0e+00 0.0e+00 5.0e-01  0  0
   0  0  0   0  0  0  0  0     0<BR>MatCoarsen             2 1.0 2.9671e-02 1.0 0.00e+00 0.0 2.4e+01 7.0e+03 3.8e+01  0  0  0  0  0   0  0  0  0  0     0<BR>MatAXPY                2 1.0 6.1393e-04 2.1 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0<BR>MatTranspose           2 1.0 1.8161e-01 1.1 0.00e+00 0.0 3.0e+01 7.1e+03 9.4e+01  0  0  0  0  0   0  0  0  0  0     0<BR>MatMatMult             2 1.0 1.1968e-01 1.0 3.31e+05 1.2
  3.6e+01 1.7e+03 1.1e+02  0  0  0  0  0   0  0  0  0  0    10<BR>MatMatMultSym          2 1.0 9.8982e-02 1.0 0.00e+00 0.0 3.0e+01 1.4e+03 9.6e+01  0  0  0  0  0   0  0  0  0  0     0<BR>MatMatMultNum          2 1.0 2.1248e-02 1.1 3.31e+05 1.2 6.0e+00 3.4e+03 1.2e+01  0  0  0  0  0   0  0  0  0  0    56<BR>MatPtAP                2 1.0 1.7070e-01 1.1 2.36e+06 1.2 5.4e+01 3.3e+03 1.1e+02  0  0  0  0  0   0  0  0  0  0    50<BR>MatPtAPSymbolic        2 1.0 1.3786e-01 1.1 0.00e+00 0.0 4.8e+01 3.1e+03 1.0
 e+02  0  0  0  0  0   0  0  0  0  0     0<BR>MatPtAPNumeric         2 1.0 4.7638e-02 2.2 2.36e+06 1.2 6.0e+00 4.8e+03 1.2e+01  0  0  0  0  0   0  0  0  0  0   180<BR>MatTrnMatMult          2 1.0 7.9914e-01 1.0 1.14e+07 1.3 3.6e+01 2.1e+04 1.2e+02  0  0  0  1  0   0  0  0  1  0    48<BR>MatGetLocalMat        10 1.0 6.7852e-03 1.1 0.00e+00 0.0 0.0e+00 0.0e+00 2.4e+01  0  0  0  0  0   0  0  0  0  0     0<BR>MatGetBrAoCol          6 1.0 3.6962e-02 2.4 0.00e+00 0.0 4.2e+01 4.5e+03 1.6e+01  0  0  0  0  0 &nbsp
 ; 0  0  0  0  0     0<BR>MatGetSymTrans         4 1.0 4.4394e-04 1.2 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0<BR>VecMDot              913 1.0 2.8472e+00 3.0 5.28e+08 1.0 0.0e+00 0.0e+00 9.1e+02  1 24  0  0  2   1 24  0  0  2   741<BR>VecNorm             1367 1.0 2.7202e+00 1.6 7.12e+06 1.0 0.0e+00 0.0e+00 1.4e+03  1  0  0  0  2   1  0  0  0  2    10<BR>VecScale            4950 1.0 6.0693e-02 1.2 1.96e+07 1.1 0.0e+00 0.0e+00 0.0e+00  0  1  0  0  0   0  
 1  0  0  0  1169<BR>VecCopy             1349 1.0 7.8685e-03 1.4 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0<BR>VecSet              4972 1.0 1.2852e-02 1.3 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0<BR>VecAXPY             7624 1.0 1.3610e-01 1.1 6.44e+07 1.2 0.0e+00 0.0e+00 0.0e+00  0  3  0  0  0   0  3  0  0  0  1676<BR>VecAYPX             7168 1.0 1.5877e-01 1.2 4.01e+07 1.2 0.0e+00 0.0e+00 0.0e+00  0  2  0  0  0 &nbsp
 ; 0  2  0  0  0   896<BR>VecMAXPY            1366 1.0 6.3739e-01 1.3 5.35e+08 1.0 0.0e+00 0.0e+00 0.0e+00  0 25  0  0  0   0 25  0  0  0  3353<BR>VecAssemblyBegin      21 1.0 4.7891e-02 1.2 0.00e+00 0.0 0.0e+00 0.0e+00 6.0e+01  0  0  0  0  0   0  0  0  0  0     0<BR>VecAssemblyEnd        21 1.0 4.5776e-05 1.4 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0<BR>VecPointwiseMult    5398 1.0 1.2278e-01 1.3 2.42e+07 1.2 0.0e+00 0.0e+00 0.0e+00  0  1  0  0  0   0  1  0  0  0   698<BR>VecScatterBegin     8107 1.0 3.9436
 e-02 1.3 0.00e+00 0.0 2.2e+04 3.4e+03 0.0e+00  0  0 99 98  0   0  0 99 98  0     0<BR>VecScatterEnd       8107 1.0 1.5414e+01344.4 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  2  0  0  0  0   2  0  0  0  0     0<BR>VecSetRandom           2 1.0 1.1868e-03 1.6 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0<BR>VecNormalize        1366 1.0 2.7266e+00 1.6 1.07e+07 1.0 0.0e+00 0.0e+00 1.4e+03  1  0  0  0  2   1  0  0  0  2    15<BR>KSPGMRESOrthog       913 1.0 8.5743e+00 1.1 1.06e+09 1.0 0.0e+00 0.0e+00 3.6e+04  3 49  0  0 60 &nb
 sp; 3 49  0  0 60   492<BR>KSPSetUp               7 1.0 2.4805e-02 1.1 0.00e+00 0.0 0.0e+00 0.0e+00 6.0e+00  0  0  0  0  0   0  0  0  0  0     0<BR>KSPSolve               1 1.0 5.7656e+01 1.0 2.30e+09 1.1 2.2e+04 3.4e+03 6.0e+04 18100100100100  18100100100100   151<BR>PCSetUp                2 1.0 2.1524e+00 1.0 2.03e+07 1.3 3.8e+02 6.0e+03 1.1e+03  1  1  2  3  2   1  1  2  3  2    33<BR>PCSetUpOnBlocks      448 1.0 1.5607e-03 1.8 1.85e+04 0.0 0.0e+00 0.0e+00 8.0e+00  0  0  0  0  0   0  0  0  0  0  &
 nbsp; 12<BR>PCApply              448 1.0 4.4331e+01 1.1 1.08e+09 1.2 1.9e+04 3.4e+03 2.3e+04 14 44 86 85 38  14 44 86 85 38    87<BR>PCGAMGgraph_AGG        2 1.0 6.1282e-01 1.0 3.36e+05 1.2 6.6e+01 5.7e+03 1.9e+02  0  0  0  0  0   0  0  0  0  0     2<BR>PCGAMGcoarse_AGG       2 1.0 8.8854e-01 1.0 1.14e+07 1.3 9.0e+01 1.2e+04 2.1e+02  0  0  0  1  0   0  0  0  1  0    44<BR>PCGAMGProl_AGG         2 1.0 6.3711e-02 1.1 0.00e+00 0.0 4.2e+01 2.3e+03 1.0e+02  0  0  0  0  0   0  0  0  0  0     0<BR>PCGAMGPOpt_AGG         2 1.0 3.4247e-01 1.0
  6.22e+06 1.2 9.6e+01 2.8e+03 3.3e+02  0  0  0  0  1   0  0  0  0  1    65<BR>------------------------------------------------------------------------------------------------------------------------</DIV>
<DIV>Memory usage is given in bytes:</DIV>
<DIV>Object Type          Creations   Destructions     Memory  Descendants' Mem.<BR>Reports information only for process 0.</DIV>
<DIV>--- Event Stage 0: Main Stage</DIV>
<DIV>              Matrix    68             68     27143956     0<BR>      Matrix Coarsen     2              2          704     0<BR>              Vector   296            296     13385864     0<BR>      Vector Scatter    18             18        11304     0<BR>        &
 nbsp;  Index Set    47             47        34816     0<BR>       Krylov Solver     7              7       554688     0<BR>      Preconditioner     7              7         3896     0<BR>         PetscRandom     2              2          704     0<BR>            &nbsp
 ; Viewer     1              0            0     0<BR>========================================================================================================================<BR>Average time to get PetscTime(): 1.3113e-06<BR>Average time for MPI_Barrier(): 9.62257e-05<BR>Average time for zero size MPI_Send(): 0.00019449<BR>#PETSc Option Table entries:<BR>-ksp_converged_reason<BR>-ksp_gmres_restart 170<BR>-ksp_rtol 1.0e-15<BR>-ksp_type gmres<BR>-log_summary<BR>-pc_gamg_agg_nsmooths 1<BR>-pc_type gamg<BR>#End of PETSc Option Table entries<BR>Compiled without FORTRAN kernels<BR>Compiled with full precision matrices (default)<BR>sizeof(short) 2 sizeof(int) 4 sizeof(long) 4 sizeof(void*) 4 sizeof(PetscScalar) 8 sizeof(PetscInt) 4<BR>Configure run at: Thu Nov  1 05:54:48 2012<BR>Configure options: --with-mpi-dir=/home/ub
 u/soft/mpich2/ --download-f-blas-lapack =1<BR>-----------------------------------------<BR>Libraries compiled on Thu Nov  1 05:54:48 2012 on ubuntu <BR>Machine characteristics: Linux-2.6.32-38-generic-i686-with-Ubuntu-10.04-lucid<BR>Using PETSc directory: /home/ubu/soft/petsc/petsc-3.3-p4<BR>Using PETSc arch: arch-linux2-c-debug<BR>-----------------------------------------</DIV>
<DIV>Using C compiler: /home/ubu/soft/mpich2/bin/mpicc  -wd1572 -g  ${COPTFLAGS} ${CFLAGS}<BR>Using Fortran compiler: /home/ubu/soft/mpich2/bin/mpif90  -g   ${FOPTFLAGS} ${FFLAGS} <BR>-----------------------------------------</DIV>
<DIV>Using include paths: -I/home/ubu/soft/petsc/petsc-3.3-p4/arch-linux2-c-debug/include -I/home/ubu/soft/petsc/petsc-3.3-p4/include -I/home/ubu/soft/petsc/petsc-3.3-p4/include -I/home/ubu/soft/petsc/petsc-3.3-p4/arch-linux2-c-debug/include -I/home/ubu/soft/mpich2/include<BR>-----------------------------------------</DIV>
<DIV>Using C linker: /home/ubu/soft/mpich2/bin/mpicc<BR>Using Fortran linker: /home/ubu/soft/mpich2/bin/mpif90<BR>Using libraries: -Wl,-rpath,/home/ubu/soft/petsc/petsc-3.3-p4/arch-linux2-c-debug/lib -L/home/ubu/soft/petsc/petsc-3.3-p4/arch-linux2-c-debug/lib -lpetsc -lpthread -Wl,-rpath,/home/ubu/soft/petsc/petsc-3.3-p4/arch-linux2-c-debug/lib -L/home/ubu/soft/petsc/petsc-3.3-p4/arch-linux2-c-debug/lib -lflapack -lfblas -L/home/ubu/soft/mpich2/lib -L/opt/intel/composer_xe_2011_sp1.10.319/compiler/lib/ia32 -L/opt/intel/composer_xe_2011_sp1.10.319/ipp/lib/ia32 -L/opt/intel/composer_xe_2011_sp1.10.319/mkl/lib/ia32 -L/opt/intel/composer_xe_2011_sp1.10.319/tbb/lib/ia32/cc4.1.0_libc2.4_kernel2.6.16.21 -L/opt/intel/composer_xe_2011_sp1.9.293/compiler/lib/ia32 -L/opt/intel/composer_xe_2011_sp1.9.293/mkl/lib/ia32 -L/usr/lib/gcc/i486-linux-gnu/4.4.3 -L/usr/lib/i486-linux-gnu -lmpichf90 -lifport -lifcore -lm -lm -ldl -lmpich -lopa -lmpl -lrt -lpthread -limf -lsvml -lipgo -ldecimal -lci
 lkrts -lstdc++ -lgcc_s -lirc -lirc_s -ldl <BR>-----------------------------------------</DIV>
<DIV> time   14.4289010000000     <BR> time   14.4289020000000     <BR> time   14.4449030000000     <BR> time   14.4809050000000     <BR></DIV>
<DIV> </DIV>
<DIV> </DIV>
<DIV>(2) using asm</DIV>
<DIV>Linear solve converged due to CONVERGED_RTOL iterations 483<BR>Norm of error  0.2591E+04 iterations   483<BR>  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000<BR>  4.866092420969481E-018  0.000000000000000E+000  0.000000000000000E+000<BR>   26.4211453778395      -4.861214483821431E-017  5.379151535696287E-018<BR>************************************************************************************************************************<BR>***             WIDEN YOUR WINDOW TO 120 CHARACTERS.  Use 'enscript -r -fCourier9' to print this document            ***<BR>************************************************************************************************************************</DIV>
<DIV>---------------------------------------------- PETSc Performance Summary: ----------------------------------------------</DIV>
<DIV>./ex4f on a arch-linux2-c-debug named ubuntu with 4 processors, by ubu Sat Nov  3 10:00:43 2012<BR>Using Petsc Release Version 3.3.0, Patch 4, Fri Oct 26 10:46:51 CDT 2012 </DIV>
<DIV>                         Max       Max/Min        Avg      Total <BR>Time (sec):           2.952e+02      1.00006   2.952e+02<BR>Objects:              2.040e+02      1.00000   2.040e+02<BR>Flops:                1.502e+09      1.00731   1.496e+09  5.983e+09<BR>Flops/sec:            5.088e+06      1.00734   5.067e+06  2.027e+07<BR>Memory:       
         2.036e+07      1.01697              8.073e+07<BR>MPI Messages:         1.960e+03      2.00000   1.470e+03  5.880e+03<BR>MPI Message Lengths:  7.738e+06      3.12820   3.474e+03  2.042e+07<BR>MPI Reductions:       4.236e+04      1.00000</DIV>
<DIV>Flop counting convention: 1 flop = 1 real number operation of type (multiply/divide/add/subtract)<BR>                            e.g., VecAXPY() for real vectors of length N --> 2N flops<BR>                            and VecAXPY() for complex vectors of length N --> 8N flops</DIV>
<DIV>Summary of Stages:   ----- Time ------  ----- Flops -----  --- Messages ---  -- Message Lengths --  -- Reductions --<BR>                        Avg     %Total     Avg     %Total   counts   %Total     Avg         %Total   counts   %Total <BR> 0:      Main Stage: 2.9517e+02 100.0%  5.9826e+09 100.0%  5.880e+03 100.0%  3.474e+03      100.0%  4.236e+04 100.0% </DIV>
<DIV>------------------------------------------------------------------------------------------------------------------------<BR>See the 'Profiling' chapter of the users' manual for details on interpreting output.<BR>Phase summary info:<BR>   Count: number of times phase was executed<BR>   Time and Flops: Max - maximum over all processors<BR>                   Ratio - ratio of maximum to minimum over all processors<BR>   Mess: number of messages sent<BR>   Avg. len: average message length<BR>   Reduct: number of global reductions<BR>   Global: entire computation<BR>   Stage: stages of a computation. Set stages with PetscLogStagePush() and PetscLogStagePop().<BR>      %T - percent time in this phase         %f - percent flops in this phase<BR>  &
 nbsp;   %M - percent messages in this phase     %L - percent message lengths in this phase<BR>      %R - percent reductions in this phase<BR>   Total Mflop/s: 10e-6 * (sum of flops over all processors)/(max time over all processors)<BR>------------------------------------------------------------------------------------------------------------------------</DIV>
<DIV><BR>      ##########################################################<BR>      #                                                        #<BR>      #                          WARNING!!!                    #<BR>      #                             &
 nbsp;                          #<BR>      #   This code was compiled with a debugging option,      #<BR>      #   To get timing results run ./configure                #<BR>      #   using --with-debugging=no, the performance will      #<BR>      #   be generally two or three times faster.              #<BR>      #                                &n
 bsp;                       #<BR>      ##########################################################</DIV>
<DIV><BR>Event                Count      Time (sec)     Flops                             --- Global ---  --- Stage ---   Total<BR>                   Max Ratio  Max     Ratio   Max  Ratio  Mess   Avg len Reduct  %T %f %M %L %R  %T %f %M %L %R Mflop/s<BR>------------------------------------------------------------------------------------------------------------------------</DIV>
<DIV>--- Event Stage 0: Main Stage</DIV>
<DIV>MatMult              485 1.0 3.1170e+00 5.9 1.36e+08 1.0 2.9e+03 3.4e+03 0.0e+00  1  9 49 48  0   1  9 49 48  0   173<BR>MatSolve             486 1.0 6.8313e-01 1.3 1.49e+08 1.1 0.0e+00 0.0e+00 0.0e+00  0 10  0  0  0   0 10  0  0  0   842<BR>MatLUFactorNum         1 1.0 6.0117e-02 1.2 1.54e+06 1.1 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0    99<BR>MatILUFactorSym        1 1.0 1.5973e-01 2.5 0.00e+00 0.0 0.0e+00 0.0e+00 1.0e+00  0  0  0  0  0   0  0  0  0  0     0<BR>MatAssemblyBegin       2 1
 .0 6.0572e-02 9.5 0.00e+00 0.0 0.0e+00 0.0e+00 2.0e+00  0  0  0  0  0   0  0  0  0  0     0<BR>MatAssemblyEnd         2 1.0 1.7764e-02 1.6 0.00e+00 0.0 1.2e+01 8.5e+02 1.9e+01  0  0  0  0  0   0  0  0  0  0     0<BR>MatGetRowIJ            1 1.0 3.8147e-06 1.3 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0<BR>MatGetSubMatrice       1 1.0 1.9829e-01 2.1 0.00e+00 0.0 3.0e+01 2.1e+04 1.0e+01  0  0  1  3  0   0  0  1  3  0     0<BR>MatGetOrdering         1 1.0 1.2739e-02 5.8 0.00e+00 0.0 0.0e+00 0
 .0e+00 4.0e+00  0  0  0  0  0   0  0  0  0  0     0<BR>MatIncreaseOvrlp       1 1.0 1.8877e-02 1.2 0.00e+00 0.0 0.0e+00 0.0e+00 2.0e+00  0  0  0  0  0   0  0  0  0  0     0<BR>VecMDot              483 1.0 4.5957e+00 2.5 5.98e+08 1.0 0.0e+00 0.0e+00 4.8e+02  1 40  0  0  1   1 40  0  0  1   521<BR>VecNorm              487 1.0 2.0843e+00 1.2 7.40e+06 1.0 0.0e+00 0.0e+00 4.9e+02  1  0  0  0  1   1  0  0  0  1    14<BR>VecScale             486 1.0 1.2140e-02 1.1 3.69e+06 1.0 0.0e+00 0.0e+00 0.0
 e+00  0  0  0  0  0   0  0  0  0  0  1217<BR>VecCopy                3 1.0 4.2915e-05 1.2 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0<BR>VecSet               979 1.0 8.1432e-03 1.1 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0<BR>VecAXPY                6 1.0 2.3413e-04 1.3 9.12e+04 1.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0  1558<BR>VecMAXPY             486 1.0 7.
 0027e-01 1.2 6.06e+08 1.0 0.0e+00 0.0e+00 0.0e+00  0 41  0  0  0   0 41  0  0  0  3460<BR>VecAssemblyBegin       1 1.0 5.7101e-03 2.5 0.00e+00 0.0 0.0e+00 0.0e+00 3.0e+00  0  0  0  0  0   0  0  0  0  0     0<BR>VecAssemblyEnd         1 1.0 3.0994e-06 1.6 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0<BR>VecScatterBegin     1457 1.0 9.1357e-02 2.6 0.00e+00 0.0 5.8e+03 3.4e+03 0.0e+00  0  0 99 97  0   0  0 99 97  0     0<BR>VecScatterEnd       1457 1.0 2.5327e+00323.7 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0&n
 bsp; 0     0<BR>VecNormalize         486 1.0 2.0858e+00 1.2 1.11e+07 1.0 0.0e+00 0.0e+00 4.9e+02  1  1  0  0  1   1  1  0  0  1    21<BR>KSPGMRESOrthog       483 1.0 1.1152e+01 1.2 1.20e+09 1.0 0.0e+00 0.0e+00 4.0e+04  3 80  0  0 94   3 80  0  0 94   429<BR>KSPSetUp               2 1.0 1.1989e-02 1.4 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00  0  0  0  0  0   0  0  0  0  0     0<BR>KSPSolve               1 1.0 2.0990e+01 1.0 1.50e+09 1.0 5.9e+03 3.5e+03 4.2e+04  7100100100100   7100100100100   285<BR>PCSetUp      &
 nbsp;         2 1.0 3.5928e-01 1.0 1.54e+06 1.1 4.2e+01 1.5e+04 3.5e+01  0  0  1  3  0   0  0  1  3  0    17<BR>PCSetUpOnBlocks        1 1.0 2.2166e-01 1.7 1.54e+06 1.1 0.0e+00 0.0e+00 6.0e+00  0  0  0  0  0   0  0  0  0  0    27<BR>PCApply              486 1.0 5.9831e+00 2.8 1.49e+08 1.1 2.9e+03 3.4e+03 1.5e+03  1 10 50 48  3   1 10 50 48  3    96<BR>------------------------------------------------------------------------------------------------------------------------</DIV>
<DIV>Memory usage is given in bytes:</DIV>
<DIV>Object Type          Creations   Destructions     Memory  Descendants' Mem.<BR>Reports information only for process 0.</DIV>
<DIV>--- Event Stage 0: Main Stage</DIV>
<DIV>              Matrix     5              5      8051336     0<BR>              Vector   182            182     11109304     0<BR>      Vector Scatter     2              2         1256     0<BR>           Index Set    10             10       121444     0<BR>   &
 nbsp;   Krylov Solver     2              2       476844     0<BR>      Preconditioner     2              2         1088     0<BR>              Viewer     1              0            0     0<BR>========================================================================================================================<BR>Average time to get PetscTime(): 9.05991e-07<BR>Average time for MPI_Barrier(): 0.000297785<BR>Average time for zero size MPI_Send()
 : 0.000174284<BR>#PETSc Option Table entries:<BR>-ksp_converged_reason<BR>-ksp_gmres_restart 170<BR>-ksp_rtol 1.0e-15<BR>-ksp_type gmres<BR>-log_summary<BR>-pc_type asm<BR>#End of PETSc Option Table entries<BR>Compiled without FORTRAN kernels<BR>Compiled with full precision matrices (default)<BR>sizeof(short) 2 sizeof(int) 4 sizeof(long) 4 sizeof(void*) 4 sizeof(PetscScalar) 8 sizeof(PetscInt) 4<BR>Configure run at: Thu Nov  1 05:54:48 2012<BR>Configure options: --with-mpi-dir=/home/ubu/soft/mpich2/ --download-f-blas-lapack =1<BR>-----------------------------------------<BR>Libraries compiled on Thu Nov  1 05:54:48 2012 on ubuntu <BR>Machine characteristics: Linux-2.6.32-38-generic-i686-with-Ubuntu-10.04-lucid<BR>Using PETSc directory: /home/ubu/soft/petsc/petsc-3.3-p4<BR>Using PETSc arch: arch-linux2-c-debug<BR>-----------------------------------------</DIV>
<DIV>Using C compiler: /home/ubu/soft/mpich2/bin/mpicc  -wd1572 -g  ${COPTFLAGS} ${CFLAGS}<BR>Using Fortran compiler: /home/ubu/soft/mpich2/bin/mpif90  -g   ${FOPTFLAGS} ${FFLAGS} <BR>-----------------------------------------</DIV>
<DIV>Using include paths: -I/home/ubu/soft/petsc/petsc-3.3-p4/arch-linux2-c-debug/include -I/home/ubu/soft/petsc/petsc-3.3-p4/include -I/home/ubu/soft/petsc/petsc-3.3-p4/include -I/home/ubu/soft/petsc/petsc-3.3-p4/arch-linux2-c-debug/include -I/home/ubu/soft/mpich2/include<BR>-----------------------------------------</DIV>
<DIV>Using C linker: /home/ubu/soft/mpich2/bin/mpicc<BR>Using Fortran linker: /home/ubu/soft/mpich2/bin/mpif90<BR>Using libraries: -Wl,-rpath,/home/ubu/soft/petsc/petsc-3.3-p4/arch-linux2-c-debug/lib -L/home/ubu/soft/petsc/petsc-3.3-p4/arch-linux2-c-debug/lib -lpetsc -lpthread -Wl,-rpath,/home/u<BR>bu/soft/petsc/petsc-3.3-p4/arch-linux2-c-debug/lib -L/home/ubu/soft/petsc/petsc-3.3-p4/arch-linux2-c-debug/lib -lflapack -lfblas -L/home/ubu/soft/mpich2/lib -L/opt/intel/composer_xe_2011_sp1.10.319/compiler/lib/ia32 -L/opt/intel/composer_xe_2011_sp1.10.319/ipp/lib/ia32 -L/opt/intel/composer_xe_2011_sp1.10.319/mkl/lib/ia32 -L/opt/intel/composer_xe_2011_sp1.10.319/tbb/lib/ia32/cc4.1.0_libc2.4_kernel2.6.16.21 -L/opt/intel/composer_xe_2011_sp1.9.293/compiler/lib/ia32 -L/opt/intel/composer_xe_2011_sp1.9.293/mkl/lib/ia32 -L/usr/lib/gcc/i486-linux-gnu/4.4.3 -L/usr/lib/i486-linux-gnu -lmpichf90 -lifport -lifcore -lm -lm -ldl -lmpich -lopa -lmpl -lrt -lpthread -limf -lsvml -lipgo -ldecimal 
 -lcilkrts -lstdc++ -lgcc_s -lirc -lirc_s -ldl <BR>-----------------------------------------</DIV>
<DIV> time   5.72435800000000     <BR> time   5.40433800000000     <BR> time   5.39233700000000     <BR> time   5.51634499999999   <BR><BR><BR><BR></DIV>
<DIV></DIV>
<DIV id="divNeteaseMailCard"></DIV>
<DIV><BR></DIV>>At 2012-11-03 23:53:42,"Jed Brown" <jedbrown@mcs.anl.gov> wrote:<BR>
<BLOCKQUOTE style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" id="isReplyContent">>Just pass it as a command line option. It gives profiling output in PetscFinalize().
<DIV class="gmail_extra"><BR><BR>
<DIV class="gmail_quote">>On Sat, Nov 3, 2012 at 10:52 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>
<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>>Is there something that need attention when setting up PETSc? The -log_summary<BR>>is no use in my system.<BR><BR><BR></DIV>
<DIV class="im">>At 2012-11-03 23:31:52,"Jed Brown" <<A href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@mcs.anl.gov</A>> wrote:<BR></DIV>
<BLOCKQUOTE style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex">
<DIV class="im">>1. <B>Always</B> send -log_summary when asking about performance. 
<DIV>>2. AMG setup costs more, the solve should be faster, especially for large problems.</DIV>
<DIV>>3. 30k degrees of freedom is not large.</DIV></DIV>
<DIV class="gmail_extra"><BR><BR>
<DIV class="gmail_quote">
<DIV class="im">>>On Sat, Nov 3, 2012 at 10:27 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>
<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 class="im">
<DIV>>>Hello,</DIV>
<DIV>>>    I have tried AMG, but there are some problems. I use the command:<BR>>>    mpiexec -n 4 ./ex4f -ksp_type gmres -pc_type gamg -pc_gamg_agg_nsmooths 1 -ksp_gmres_restart 170 -ksp_rtol 1.0e-15 ->>ksp_converged_reason.<BR>>>    The matrix has a size of 30000. However, compared with -pc_type asm,<BR>>>the amg need more time:asm needs 4.9s, amg needs 13.7s. I did several tests<BR>>>and got the same conclusion. When it begins, the screen shows the information: <BR>>>[0]PCSetData_AGG bs=1 MM=7601. I do not know the meaning. And if there is some<BR>>>parameters that affect the performance of AMG?<BR>>>    Besides, I want to confirm a conception. In my view, AMG itself can be a solver<BR>>>like gmres. It can also be used as a preconditioner like jacobi and is used by combining<BR>>>with other solver. Is it right? If it is right, how use AMG solver?<BR
 >>>    My codes are attached.<BR>>>    Thanks.<BR>>>                                             Jim</DIV>
<DIV><BR><BR> </DIV></DIV>
<DIV class="im">>At 2012-11-01 22:00:28,"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">
<DIV class="gmail_extra">>Yes, it's faster to understand this error message than to have "mysteriously slow performance".</DIV>
<DIV class="gmail_extra"><BR></DIV>
<DIV class="gmail_extra">><EM>* Preallocation routines now automatically set MAT_NEW_NONZERO_ALLOCATION_ERR, if you intentionally preallocate less than necessary then </EM>><EM>use MatSetOption(mat,MAT_NEW_NONZERO_ALLOCATION_ERR,PETSC_FALSE) to disable the error generation.<BR></EM><BR>><A href="http://www.mcs.anl.gov/petsc/documentation/changes/33.html" target="_blank">http://www.mcs.anl.gov/petsc/documentation/changes/33.html</A><BR><BR>
<DIV class="gmail_quote">>On Thu, Nov 1, 2012 at 8:57 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>
<BLOCKQUOTE style="BORDER-LEFT: rgb(204,204,204) 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" class="gmail_quote">
<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></BLOCKQUOTE></DIV><BR></DIV></BLOCKQUOTE></DIV></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></DIV></BLOCKQUOTE></div><br><br><span title="neteasefooter"><span id="netease_mail_footer"></span></span>