[petsc-users] TS and snes_grid_sequence

[Jed Brown]

Barry Smith <bsmith at mcs.anl.gov> writes:

>> On Feb 10, 2016, at 4:19 PM, Jed Brown <jed at jedbrown.org> wrote:
>> 
>> Barry Smith <bsmith at mcs.anl.gov> writes:
>>>     We do not currently have library code to do grid sequence within a time-step. I believe that providing this requires a great deal of "replumbing" of the TS solvers to make them "grid sequence aware". This is true I believe for both traditional grid sequencing and FAS (I would be happy if Emil or Jed could prove me wrong).  There are seemingly multiple ways to do the grid sequencing: 
>> 
>> I'm pretty sure you could use TS solved by FAS when Peter and I were
>> working on it.  But my test just failed and it looks like some work
>> vectors aren't being obtained correctly.  I'm afraid I can't dig into it
>> right now.
>
>    Is there an example in the repository?

$ ./ex25 -ts_type beuler -snes_type fas -da_refine 1
[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
[0]PETSC ERROR: Arguments are incompatible
[0]PETSC ERROR: Vector local size 42 is not compatible with DMDA local sizes 22 22

[0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting.
[0]PETSC ERROR: Petsc Development GIT revision: v3.6.3-3237-gc58487d  GIT Date: 2016-01-27 18:34:57 -0600
[0]PETSC ERROR: ./ex25 on a mpich-clang named batura by jed Wed Feb 10 16:13:34 2016
[0]PETSC ERROR: Configure options --download-chaco --download-exodusii --download-sundials --download-triangle --with-fc=0 --with-hdf5-dir=/opt/mpich --with-mpi-dir=/home/jed/usr/ccache/mpich-clang --with-netcdf --with-pcbddc --with-shared-libraries --with-visibility --with-x --with-yaml PETSC_ARCH=mpich-clang
[0]PETSC ERROR: #1 DMDAVecGetArray() line 70 in /home/jed/petsc/src/dm/impls/da/dagetarray.c
[0]PETSC ERROR: #2 FormRHSFunction() line 210 in /home/jed/petsc/src/ts/examples/tutorials/ex25.c
[0]PETSC ERROR: #3 TSComputeRHSFunction() line 627 in /home/jed/petsc/src/ts/interface/ts.c
[0]PETSC ERROR: #4 TSComputeIFunction() line 837 in /home/jed/petsc/src/ts/interface/ts.c
[0]PETSC ERROR: #5 SNESTSFormFunction_Theta() line 508 in /home/jed/petsc/src/ts/impls/implicit/theta/theta.c
[0]PETSC ERROR: #6 SNESTSFormFunction() line 4529 in /home/jed/petsc/src/ts/interface/ts.c
[0]PETSC ERROR: #7 SNESComputeFunction() line 2136 in /home/jed/petsc/src/snes/interface/snes.c
[0]PETSC ERROR: #8 SNESFASCoarseCorrection() line 659 in /home/jed/petsc/src/snes/impls/fas/fas.c
[0]PETSC ERROR: #9 SNESFASCycle_Multiplicative() line 811 in /home/jed/petsc/src/snes/impls/fas/fas.c
[0]PETSC ERROR: #10 SNESSolve_FAS() line 980 in /home/jed/petsc/src/snes/impls/fas/fas.c
[0]PETSC ERROR: #11 SNESSolve() line 3992 in /home/jed/petsc/src/snes/interface/snes.c
[0]PETSC ERROR: #12 TSStep_Theta() line 248 in /home/jed/petsc/src/ts/impls/implicit/theta/theta.c
[0]PETSC ERROR: #13 TSStep() line 3338 in /home/jed/petsc/src/ts/interface/ts.c
[0]PETSC ERROR: #14 TSSolve() line 3544 in /home/jed/petsc/src/ts/interface/ts.c
[0]PETSC ERROR: #15 main() line 119 in /home/jed/petsc/src/ts/examples/tutorials/ex25.c
[0]PETSC ERROR: PETSc Option Table entries:
[0]PETSC ERROR: -da_refine 1
[0]PETSC ERROR: -malloc_test
[0]PETSC ERROR: -snes_type fas
[0]PETSC ERROR: -ts_type beuler
[0]PETSC ERROR: ----------------End of Error Message -------send entire error message to petsc-maint at mcs.anl.gov----------


Problem here is that TSComputeIFunction is getting the Frhs work vec
From the TS rather than something associated with the DM.  I don't know
if there are other issues later.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 818 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160210/af2a4d38/attachment-0001.pgp>