[petsc-users] Issue of mg_coarse_ksp not converge
Wenbo Zhao
zhaowenbo.npic at gmail.com
Tue Oct 3 02:55:24 CDT 2017
On Tue, Oct 3, 2017 at 1:49 AM, Mark Adams <mfadams at lbl.gov> wrote:
> Well that is strange, the PETSc tests work.
>
> Wenbo, could you please:
>
> > git clone -b gamg-fix-eig-err https://bitbucket.org/petsc/petsc petsc2
> > cd petsc2
>
> and reconfigure, make, and then run your test without the
> -st_gamg_est_ksp_error_if_not_converged 0 fix, and see if this fixes the
> problem.
>
> Don't forget to set PETSC_DIR=..../petsc2
>
> >git clone -b mark/gamg-fix-eig-err https://bitbucket.org/petsc/petsc
petsc2
It works well without the -st_gamg_est_ksp_error_if_not_converged 0.
mpirun -n 1 ./step-41 \
-st_ksp_type gmres -st_ksp_view -st_ksp_monitor \
-st_pc_type gamg -st_pc_gamg_type agg -st_pc_gamg_agg_nsmooths 1 \
-mata AMAT.dat -matb BMAT.dat \
-st_gamg_est_ksp_view -st_gamg_est_ksp_monitor \
-st_gamg_est_ksp_converged_reason \
-eps_nev 1 -eps_ncv 10 -eps_monitor -log_view > log_smooth 2>&1
Thanks,
Wenbo
> If you have time and this works, you could do a 'git checkout master' and
> remake, and retest. You should not have to reconfigure. I have tested
> master on petsc tests. I don't understand how this happened.
>
> Thanks,
> Mark
>
>
> On Mon, Oct 2, 2017 at 12:32 PM, Wenbo Zhao <zhaowenbo.npic at gmail.com>
> wrote:
>
>>
>>
>> On Tue, Oct 3, 2017 at 12:23 AM, Matthew Knepley <knepley at gmail.com>
>> wrote:
>>
>>> On Mon, Oct 2, 2017 at 11:56 AM, Wenbo Zhao <zhaowenbo.npic at gmail.com>
>>> wrote:
>>>
>>>>
>>>>
>>>> On Mon, Oct 2, 2017 at 11:49 PM, Mark Adams <mfadams at lbl.gov> wrote:
>>>>
>>>>> Whenbo, do you build your PETSc?
>>>>>
>>>>> Yes.
>>>> My configure option is listed below
>>>> ./configure --with-mpi=1 --with-shared-libraries=1 \
>>>> --with-64-bit-indices=1 --with-debugging=1
>>>>
>>>> And I set PETSC_DIR, PETSC_ARCH and SLEPC_DIR in my ~/.bashrc.
>>>>
>>>>
>>>> The Makefile for my problem is listed below,
>>>>
>>>> PETSC_ARCH = arch-linux2-c-debug
>>>> PETSC_DIR = /home/zhaowenbo/research/petsc/petsc_git
>>>> SLEPC_DIR = /home/zhaowenbo/research/slepc/slepc_git
>>>> #PETSC_DIR = /home/zhaowenbo/research/petsc/petsc-3.7.4
>>>> #SLEPC_DIR = /home/zhaowenbo/research/slepc/slepc-3.7.3
>>>> HYPRE_DIR = /usr/local/hypre
>>>> #
>>>> DEBUG_OPT = -g
>>>> COMP_FLAGS = -fPIC -Wall \
>>>> -I${SLEPC_DIR}/include -I${SLEPC_DIR}/${PETSC_ARCH}/include \
>>>> -I${PETSC_DIR}/include -I${PETSC_DIR}/${PETSC_ARCH}/include \
>>>> -Isrc
>>>>
>>>> LINK_FLAGS = -fPIC -Wall \
>>>> -Wl,-rpath,${SLEPC_DIR}/${PETSC_ARCH}/lib
>>>> -L${SLEPC_DIR}/${PETSC_ARCH}/lib -lslepc \
>>>> -Wl,-rpath,${PETSC_DIR}/${PETSC_ARCH}/lib
>>>> -L${PETSC_DIR}/${PETSC_ARCH}/lib -lpetsc \
>>>> -llapack -lblas -lhwloc -lm -lgfortran -lquadmath
>>>>
>>>> step-41: src/main.o src/readinp.o src/base.o src/sp3.o src/diffu.o
>>>> mpicxx -o step-41 $^ ${LINK_FLAGS} ${DEBUG_OPT}
>>>>
>>>> src/main.o: src/main.c
>>>> mpicxx -o src/main.o -c $^ ${COMP_FLAGS} ${DEBUG_OPT}
>>>>
>>>> src/readinp.o: src/readinp.c
>>>> mpicxx -o src/readinp.o -c $^ ${COMP_FLAGS} ${DEBUG_OPT}
>>>>
>>>> src/sp3.o: src/sp3.c
>>>> mpicxx -o src/sp3.o -c $^ ${COMP_FLAGS} ${DEBUG_OPT}
>>>>
>>>> src/diffu.o: src/diffu.c
>>>> mpicxx -o src/diffu.o -c $^ ${COMP_FLAGS} ${DEBUG_OPT}
>>>>
>>>> src/base.o: src/base.c
>>>> mpicxx -o src/base.o -c $^ ${COMP_FLAGS} ${DEBUG_OPT}
>>>>
>>>>
>>>> clean:
>>>> rm step-41 src/main.o src/readinp.o src/sp3.o src/diffu.o src/base.o
>>>>
>>>> runkr_smooth:
>>>> mpirun -n ${NCORE} ./step-41 \
>>>> -st_ksp_type gmres -st_ksp_view -st_ksp_monitor \
>>>> -st_pc_type gamg -st_pc_gamg_type agg -st_pc_gamg_agg_nsmooths 1 \
>>>> -mata AMAT.dat -matb BMAT.dat \
>>>> -st_gamg_est_ksp_view -st_gamg_est_ksp_monitor \
>>>> -st_gamg_est_ksp_converged_reason \
>>>>
>>>
>>> Add -st_gamg_est_ksp_error_if_not_converged 0
>>>
>>> Thanks,
>>>
>>> Matt
>>>
>>
>> It works after adding -st_gamg_est_ksp_error_if_not_converged 0.
>>
>> Thanks,
>> Wenbo
>>
>>
>>>
>>> -eps_nev 1 -eps_ncv 10 -eps_monitor -log_view > log_smooth 2>&1
>>>>
>>>> runkr_nonsmooth:
>>>> mpirun -n ${NCORE} ./step-41 \
>>>> -st_ksp_type gmres -st_ksp_view -st_ksp_monitor \
>>>> -st_pc_type gamg -st_pc_gamg_type agg -st_pc_gamg_agg_nsmooths 0 \
>>>> -mata AMAT.dat -matb BMAT.dat \
>>>> -st_gamg_est_ksp_view -st_gamg_est_ksp_monitor \
>>>> -st_gamg_est_ksp_converged_reason \
>>>> -eps_nev 1 -eps_ncv 10 -eps_monitor -log_view > log_nonsmooth 2>&1
>>>>
>>>>
>>>> Thanks,
>>>> Wenbo
>>>>
>>>>
>>>>
>>>>
>>>>> On Mon, Oct 2, 2017 at 11:45 AM, Mark Adams <mfadams at lbl.gov> wrote:
>>>>>
>>>>>> This is normal:
>>>>>>
>>>>>> Linear st_gamg_est_ solve did not converge due to DIVERGED_ITS
>>>>>> iterations 10
>>>>>>
>>>>>> It looks like ksp->errorifnotconverged got set somehow. If the
>>>>>> default changed in KSP then (SAGG) GAMG would not ever work.
>>>>>>
>>>>>> I assume you don't have a .petscrc file with more (crazy) options in
>>>>>> it ...
>>>>>>
>>>>>>
>>>>>> On Mon, Oct 2, 2017 at 11:39 AM, Wenbo Zhao <zhaowenbo.npic at gmail.com
>>>>>> > wrote:
>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> On Mon, Oct 2, 2017 at 11:30 PM, Matthew Knepley <knepley at gmail.com>
>>>>>>> wrote:
>>>>>>>
>>>>>>>> On Mon, Oct 2, 2017 at 11:15 AM, Mark Adams <mfadams at lbl.gov>
>>>>>>>> wrote:
>>>>>>>>
>>>>>>>>> non-smoothed aggregation is converging very fast. smoothed fails
>>>>>>>>> in the eigen estimator.
>>>>>>>>>
>>>>>>>>> Run this again with -st_gamg_est_ksp_view and
>>>>>>>>> -st_gamg_est_ksp_monitor, and see if you get more output (I'm not 100% sure
>>>>>>>>> about these args).
>>>>>>>>>
>>>>>>>>
>>>>>>>> I also want -st_gamg_est_ksp_converged_reason
>>>>>>>>
>>>>>>>> Thanks,
>>>>>>>>
>>>>>>>> Matt
>>>>>>>>
>>>>>>> $make NCORE=1 runkr_smooth
>>>>>>> mpirun -n 1 ./step-41 \
>>>>>>> -st_ksp_type gmres -st_ksp_view -st_ksp_monitor \
>>>>>>> -st_pc_type gamg -st_pc_gamg_type agg -st_pc_gamg_agg_nsmooths 1 \
>>>>>>> -mata AMAT.dat -matb BMAT.dat \
>>>>>>> -st_gamg_est_ksp_view -st_gamg_est_ksp_monitor \
>>>>>>> -st_gamg_est_ksp_converged_reason \
>>>>>>> -eps_nev 1 -eps_ncv 10 -eps_monitor -log_view > log_smooth 2>&1
>>>>>>> makefile:43: recipe for target 'runkr_smooth' failed
>>>>>>> make: *** [runkr_smooth] Error 91
>>>>>>>
>>>>>>> Thanks
>>>>>>> Wenbo
>>>>>>>
>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>>>
>>>>>>>>> On Mon, Oct 2, 2017 at 11:06 AM, Wenbo Zhao <
>>>>>>>>> zhaowenbo.npic at gmail.com> wrote:
>>>>>>>>>
>>>>>>>>>> Matt,
>>>>>>>>>>
>>>>>>>>>> Test 1 nonsmooth
>>>>>>>>>> zhaowenbo at ubuntu:~/test_slepc/SPARK/spark$ make NCORE=1
>>>>>>>>>> runkr_nonsmooth
>>>>>>>>>> mpirun -n 1 ./step-41 \
>>>>>>>>>> -st_ksp_type gmres -st_ksp_view -st_ksp_monitor \
>>>>>>>>>> -st_pc_type gamg -st_pc_gamg_type agg -st_pc_gamg_agg_nsmooths
>>>>>>>>>> 0 \
>>>>>>>>>> -mata AMAT.dat -matb BMAT.dat \
>>>>>>>>>> -eps_nev 1 -eps_ncv 10 -eps_monitor -log_view >
>>>>>>>>>> log_nonsmooth 2>&1
>>>>>>>>>>
>>>>>>>>>> Test 2 smooth
>>>>>>>>>> zhaowenbo at ubuntu:~/test_slepc/SPARK/spark$ make NCORE=1
>>>>>>>>>> runkr_smooth
>>>>>>>>>> mpirun -n 1 ./step-41 \
>>>>>>>>>> -st_ksp_type gmres -st_ksp_view -st_ksp_monitor \
>>>>>>>>>> -st_pc_type gamg -st_pc_gamg_type agg -st_pc_gamg_agg_nsmooths
>>>>>>>>>> 1 \
>>>>>>>>>> -mata AMAT.dat -matb BMAT.dat \
>>>>>>>>>> -eps_nev 1 -eps_ncv 10 -eps_monitor -log_view > log_smooth
>>>>>>>>>> 2>&1
>>>>>>>>>> makefile:43: recipe for target 'runkr_smooth' failed
>>>>>>>>>> make: *** [runkr_smooth] Error 91
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>> Thanks,
>>>>>>>>>>
>>>>>>>>>> Wenbo
>>>>>>>>>>
>>>>>>>>>> On Mon, Oct 2, 2017 at 10:48 PM, Matthew Knepley <
>>>>>>>>>> knepley at gmail.com> wrote:
>>>>>>>>>>
>>>>>>>>>>> On Mon, Oct 2, 2017 at 10:43 AM, Wenbo Zhao <
>>>>>>>>>>> zhaowenbo.npic at gmail.com> wrote:
>>>>>>>>>>>
>>>>>>>>>>>> Mark,
>>>>>>>>>>>>
>>>>>>>>>>>> Thanks for your reply.
>>>>>>>>>>>>
>>>>>>>>>>>> On Mon, Oct 2, 2017 at 9:51 PM, Mark Adams <mfadams at lbl.gov>
>>>>>>>>>>>> wrote:
>>>>>>>>>>>>
>>>>>>>>>>>>> Please send the output with -st_ksp_view and -st_ksp_monitor
>>>>>>>>>>>>> and we can start to debug it.
>>>>>>>>>>>>>
>>>>>>>>>>>>> Test 1 with nonsmooth and preonly is OK
>>>>>>>>>>>> zhaowenbo at ubuntu:~/test_slepc/SPARK/spark$ make NCORE=1
>>>>>>>>>>>> runkr_nonsmooth
>>>>>>>>>>>> mpirun -n 1 ./step-41 \
>>>>>>>>>>>> -st_ksp_type gmres -st_ksp_view -st_ksp_monitor \
>>>>>>>>>>>> -st_pc_type gamg -st_pc_gamg_type agg
>>>>>>>>>>>> -st_pc_gamg_agg_nsmooths 0 \
>>>>>>>>>>>> -st_ksp_view -mata AMAT.dat -matb BMAT.dat \
>>>>>>>>>>>> -st_mg_coarse_ksp_type preonly -st_mg_coarse_ksp_monitor \
>>>>>>>>>>>> -eps_nev 1 -eps_ncv 10 -eps_monitor -log_view >
>>>>>>>>>>>> log_nonsmooth 2>&1
>>>>>>>>>>>>
>>>>>>>>>>>> Test 2 smooth and preonly is not OK
>>>>>>>>>>>> zhaowenbo at ubuntu:~/test_slepc/SPARK/spark$ make NCORE=1
>>>>>>>>>>>> runkr_smooth
>>>>>>>>>>>> mpirun -n 1 ./step-41 \
>>>>>>>>>>>> -st_ksp_type gmres -st_ksp_view -st_ksp_monitor \
>>>>>>>>>>>> -st_pc_type gamg -st_pc_gamg_type agg
>>>>>>>>>>>> -st_pc_gamg_agg_nsmooths 1 \
>>>>>>>>>>>> -st_ksp_view -mata AMAT.dat -matb BMAT.dat \
>>>>>>>>>>>> -st_mg_coarse_ksp_type preonly -st_mg_coarse_ksp_monitor \
>>>>>>>>>>>> -eps_nev 1 -eps_ncv 10 -eps_monitor -log_view > log_smooth
>>>>>>>>>>>> 2>&1
>>>>>>>>>>>> makefile:43: recipe for target 'runkr_smooth' failed
>>>>>>>>>>>> make: *** [runkr_smooth] Error 91
>>>>>>>>>>>>
>>>>>>>>>>>> Test 3 nonsmooth and gmres is not OK
>>>>>>>>>>>> zhaowenbo at ubuntu:~/test_slepc/SPARK/spark$ make NCORE=1
>>>>>>>>>>>> runkr_gmres
>>>>>>>>>>>> mpirun -n 1 ./step-41 \
>>>>>>>>>>>> -st_ksp_type gmres -st_ksp_view -st_ksp_monitor \
>>>>>>>>>>>> -st_pc_type gamg -st_pc_gamg_type agg
>>>>>>>>>>>> -st_pc_gamg_agg_nsmooths 0 \
>>>>>>>>>>>> -st_ksp_view -mata AMAT.dat -matb BMAT.dat \
>>>>>>>>>>>> -st_mg_coarse_ksp_type gmres -st_mg_coarse_ksp_monitor
>>>>>>>>>>>> -st_mg_coarse_ksp_rtol 1.0e-6 \
>>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>> DO NOT DO THIS. Please send the output where you do NOTHING to
>>>>>>>>>>> the coarse solver.
>>>>>>>>>>>
>>>>>>>>>>> Thanks,
>>>>>>>>>>>
>>>>>>>>>>> Matt
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>>> -eps_nev 1 -eps_ncv 10 -eps_monitor -log_view > log_gmres
>>>>>>>>>>>> 2>&1
>>>>>>>>>>>> makefile:59: recipe for target 'runkr_gmres' failed
>>>>>>>>>>>> make: *** [runkr_gmres] Error 91
>>>>>>>>>>>>
>>>>>>>>>>>> log-files is attached.
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>> You mentioned that B is not symmetric. I assume it is elliptic
>>>>>>>>>>>>> (diffusion). Where does the asymmetry come from?
>>>>>>>>>>>>>
>>>>>>>>>>>>>
>>>>>>>>>>>> It is a two-group diffusion equations, where group denotes
>>>>>>>>>>>> neutron enegry discretisation.
>>>>>>>>>>>> Matrix B consists of neutron diffusion/leakage term, removal
>>>>>>>>>>>> term and minus neutron scatter source term between different energies, when
>>>>>>>>>>>> matrix A denotes neutron fission source.
>>>>>>>>>>>>
>>>>>>>>>>>> Diffusion term(Laplace operator) is elliptic and symmetric.
>>>>>>>>>>>> Removal term is diagonal only. However scatter term is asymmetry since
>>>>>>>>>>>> scatter term from high energy to low energy is far greater than the term
>>>>>>>>>>>> from low to high.
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>> Wenbo
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>>> On Mon, Oct 2, 2017 at 9:39 AM, Wenbo Zhao <
>>>>>>>>>>>>> zhaowenbo.npic at gmail.com> wrote:
>>>>>>>>>>>>>
>>>>>>>>>>>>>> Matt,
>>>>>>>>>>>>>> Thanks for your reply.
>>>>>>>>>>>>>> For the defalt option doesnt work firstly( -st_ksp_type gmres
>>>>>>>>>>>>>> -st_pc_type gamg -st_pc_gamg_type agg -st_pc_gamg_agg_nsmooths 1), I tried
>>>>>>>>>>>>>> to test those options.
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> Wenbo
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> On Mon, Oct 2, 2017 at 9:08 PM, Matthew Knepley <
>>>>>>>>>>>>>> knepley at gmail.com> wrote:
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> On Mon, Oct 2, 2017 at 8:30 AM, Wenbo Zhao <
>>>>>>>>>>>>>>> zhaowenbo.npic at gmail.com> wrote:
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> Matt
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> Because I am not clear about what will happen using
>>>>>>>>>>>>>>>> 'preonly' for large scale problem.
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> The size of the problem has nothing to do with 'preonly'.
>>>>>>>>>>>>>>> All it means is to apply a preconditioner without a Krylov solver.
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> It seems to use a direct solver from below,
>>>>>>>>>>>>>>>> http://www.mcs.anl.gov/petsc/p
>>>>>>>>>>>>>>>> etsc-current/docs/manualpages/KSP/KSPPREONLY.html
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> However, I still cannot understand why you would change the
>>>>>>>>>>>>>>> default?
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> Matt
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> Thanks!
>>>>>>>>>>>>>>>> Wenbo
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> On Mon, Oct 2, 2017 at 5:09 PM, Matthew Knepley <
>>>>>>>>>>>>>>>> knepley at gmail.com> wrote:
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> On Sun, Oct 1, 2017 at 9:53 PM, Wenbo Zhao <
>>>>>>>>>>>>>>>>> zhaowenbo.npic at gmail.com> wrote:
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>> Matt,
>>>>>>>>>>>>>>>>>> Thanks for your reply.
>>>>>>>>>>>>>>>>>> It DOES make no sense for this problem.
>>>>>>>>>>>>>>>>>> But I am not clear about the 'preonly' option. Which
>>>>>>>>>>>>>>>>>> solver is used in preonly? I wonder if 'preonly' is suitable for large
>>>>>>>>>>>>>>>>>> scale problem such as 400,000,000 unknowns.
>>>>>>>>>>>>>>>>>> So I tried 'gmres' option and found these error messages.
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> I mean, why are you setting this at all. Just do not set
>>>>>>>>>>>>>>>>> the coarse solver. The default should work fine.
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> Thanks,
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> Matt
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>> Could you give me some suggestions?
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>> Thanks.
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>> Wenbo
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>> On Mon, Oct 2, 2017 at 12:34 AM, Matthew Knepley <
>>>>>>>>>>>>>>>>>> knepley at gmail.com> wrote:
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> On Sun, Oct 1, 2017 at 6:49 AM, Wenbo Zhao <
>>>>>>>>>>>>>>>>>>> zhaowenbo.npic at gmail.com> wrote:
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> Hi,
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> I met some questions when I use PETSC/SLEPC to solve
>>>>>>>>>>>>>>>>>>>> two-group neutron diffusion equations with finite difference method. The
>>>>>>>>>>>>>>>>>>>> grid is 3*3*3, when DOF on each points is 2. So the matrix size is 54*54.
>>>>>>>>>>>>>>>>>>>> It is generalized eigenvalue problem Ax=\lamda Bx,
>>>>>>>>>>>>>>>>>>>> where B is diagonally dominant matrix but not symmetry.
>>>>>>>>>>>>>>>>>>>> EPS is set as below,
>>>>>>>>>>>>>>>>>>>> ierr = EPSSetProblemType(eps,EPS_GNHE
>>>>>>>>>>>>>>>>>>>> P);CHKERRQ(ierr);¬
>>>>>>>>>>>>>>>>>>>> ierr = EPSSetWhichEigenpairs(eps,EPS_
>>>>>>>>>>>>>>>>>>>> LARGEST_REAL);CHKERRQ(ierr);¬
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> Krylovschur is used as eps sovler. GAMG is used as PC.
>>>>>>>>>>>>>>>>>>>> I tried agg_nsmooths and mg_coarse_ksp_type. Only
>>>>>>>>>>>>>>>>>>>> non-smooths and preonly is OK.
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> Why are you setting the coarse solver. This makes no
>>>>>>>>>>>>>>>>>>> sense.
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> Thanks,
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> Matt
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> Test 1
>>>>>>>>>>>>>>>>>>>> $ make NCORE=1 runkr_nonsmooth
>>>>>>>>>>>>>>>>>>>> mpirun -n 1 ./step-41 \
>>>>>>>>>>>>>>>>>>>> -st_ksp_type gmres \
>>>>>>>>>>>>>>>>>>>> -st_pc_type gamg -st_pc_gamg_type agg
>>>>>>>>>>>>>>>>>>>> -st_pc_gamg_agg_nsmooths 0 \
>>>>>>>>>>>>>>>>>>>> -st_ksp_view -mata AMAT.dat -matb BMAT.dat \
>>>>>>>>>>>>>>>>>>>> -st_mg_coarse_ksp_type preonly
>>>>>>>>>>>>>>>>>>>> -st_mg_coarse_ksp_monitor \
>>>>>>>>>>>>>>>>>>>> -eps_nev 1 -eps_ncv 10 -eps_monitor -log_view >
>>>>>>>>>>>>>>>>>>>> log_nonsmooth 2>&1
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> Test 2
>>>>>>>>>>>>>>>>>>>> $ make NCORE=1 runkr_smooth
>>>>>>>>>>>>>>>>>>>> mpirun -n 1 ./step-41 \
>>>>>>>>>>>>>>>>>>>> -st_ksp_type gmres \
>>>>>>>>>>>>>>>>>>>> -st_pc_type gamg -st_pc_gamg_type agg
>>>>>>>>>>>>>>>>>>>> -st_pc_gamg_agg_nsmooths 1 \
>>>>>>>>>>>>>>>>>>>> -st_ksp_view -mata AMAT.dat -matb BMAT.dat \
>>>>>>>>>>>>>>>>>>>> -st_mg_coarse_ksp_type preonly
>>>>>>>>>>>>>>>>>>>> -st_mg_coarse_ksp_monitor \
>>>>>>>>>>>>>>>>>>>> -eps_nev 1 -eps_ncv 10 -eps_monitor -log_view >
>>>>>>>>>>>>>>>>>>>> log_smooth 2>&1
>>>>>>>>>>>>>>>>>>>> makefile:43: recipe for target 'runkr_smooth' failed
>>>>>>>>>>>>>>>>>>>> make: *** [runkr_smooth] Error 91
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> Test 3
>>>>>>>>>>>>>>>>>>>> $ make NCORE=1 runkr_gmres
>>>>>>>>>>>>>>>>>>>> mpirun -n 1 ./step-41 \
>>>>>>>>>>>>>>>>>>>> -st_ksp_type gmres \
>>>>>>>>>>>>>>>>>>>> -st_pc_type gamg -st_pc_gamg_type agg
>>>>>>>>>>>>>>>>>>>> -st_pc_gamg_agg_nsmooths 0 \
>>>>>>>>>>>>>>>>>>>> -st_ksp_view -mata AMAT.dat -matb BMAT.dat \
>>>>>>>>>>>>>>>>>>>> -st_mg_coarse_ksp_type gmres
>>>>>>>>>>>>>>>>>>>> -st_mg_coarse_ksp_monitor -st_mg_coarse_ksp_rtol 1.0e-6 \
>>>>>>>>>>>>>>>>>>>> -eps_nev 1 -eps_ncv 10 -eps_monitor -log_view >
>>>>>>>>>>>>>>>>>>>> log_gmres 2>&1
>>>>>>>>>>>>>>>>>>>> makefile:59: recipe for target 'runkr_gmres' failed
>>>>>>>>>>>>>>>>>>>> make: *** [runkr_gmres] Error 91
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> Log files were attched.
>>>>>>>>>>>>>>>>>>>> The matrix file were also attched as AMAT.dat and
>>>>>>>>>>>>>>>>>>>> BMAT.dat.
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> Is it correct? Or something wrong with my code or
>>>>>>>>>>>>>>>>>>>> commad-line?
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> Thanks!
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> Wenbo
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> --
>>>>>>>>>>>>>>>>>>> What most experimenters take for granted before they
>>>>>>>>>>>>>>>>>>> begin their experiments is infinitely more interesting than any results to
>>>>>>>>>>>>>>>>>>> which their experiments lead.
>>>>>>>>>>>>>>>>>>> -- Norbert Wiener
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> https://www.cse.buffalo.edu/~knepley/
>>>>>>>>>>>>>>>>>>> <http://www.caam.rice.edu/%7Emk51/>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> --
>>>>>>>>>>>>>>>>> What most experimenters take for granted before they begin
>>>>>>>>>>>>>>>>> their experiments is infinitely more interesting than any results to which
>>>>>>>>>>>>>>>>> their experiments lead.
>>>>>>>>>>>>>>>>> -- Norbert Wiener
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> https://www.cse.buffalo.edu/~knepley/
>>>>>>>>>>>>>>>>> <http://www.caam.rice.edu/%7Emk51/>
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> --
>>>>>>>>>>>>>>> What most experimenters take for granted before they begin
>>>>>>>>>>>>>>> their experiments is infinitely more interesting than any results to which
>>>>>>>>>>>>>>> their experiments lead.
>>>>>>>>>>>>>>> -- Norbert Wiener
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> https://www.cse.buffalo.edu/~knepley/
>>>>>>>>>>>>>>> <http://www.caam.rice.edu/%7Emk51/>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>> --
>>>>>>>>>>> What most experimenters take for granted before they begin their
>>>>>>>>>>> experiments is infinitely more interesting than any results to which their
>>>>>>>>>>> experiments lead.
>>>>>>>>>>> -- Norbert Wiener
>>>>>>>>>>>
>>>>>>>>>>> https://www.cse.buffalo.edu/~knepley/
>>>>>>>>>>> <http://www.caam.rice.edu/%7Emk51/>
>>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>> --
>>>>>>>> What most experimenters take for granted before they begin their
>>>>>>>> experiments is infinitely more interesting than any results to which their
>>>>>>>> experiments lead.
>>>>>>>> -- Norbert Wiener
>>>>>>>>
>>>>>>>> https://www.cse.buffalo.edu/~knepley/
>>>>>>>> <http://www.caam.rice.edu/%7Emk51/>
>>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>
>>>>>
>>>>
>>>
>>>
>>> --
>>> What most experimenters take for granted before they begin their
>>> experiments is infinitely more interesting than any results to which their
>>> experiments lead.
>>> -- Norbert Wiener
>>>
>>> https://www.cse.buffalo.edu/~knepley/
>>> <http://www.caam.rice.edu/%7Emk51/>
>>>
>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20171003/04b01945/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: log_new5.tgz
Type: application/x-gzip
Size: 12481 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20171003/04b01945/attachment-0001.tgz>
More information about the petsc-users
mailing list