Barry Smith <[email protected]> writes: >> On Feb 10, 2016, at 4:19 PM, Jed Brown <[email protected]> wrote: >> >> Barry Smith <[email protected]> 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 [email protected] 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.
signature.asc
Description: PGP signature
