I'd like to have some automatic way to use coloring with TS as long as I provide a residual evaluation function. In the current approach, if I don't do anything special, it fails because there are no TS functions on the DM. Barry, how would you like to go about fixing this?<div>
<br><div><br><div><div><div>./ex26 -da_grid_x 20 -da_grid_y 20 -lidvelocity 100 -grashof 1e3 -ts_max_steps 100 -ts_rtol 1e-3 -ts_atol 1e-3 -ts_type rosw -ts_rosw_type ra3pw -ts_monitor -ts_monitor_solution_vtk 'foo-%03D.vts'</div>
<div>20x20 grid, lid velocity = 100, prandtl # = 1, grashof # = 1000</div><div>0 TS dt 0.025 time 0</div></div><div><div>[0]PETSC ERROR: --------------------- Error Message ------------------------------------</div><div>[0]PETSC ERROR: Object is in wrong state!</div>
<div>[0]PETSC ERROR: Must call MatFDColoringSetFunction()!</div><div>[0]PETSC ERROR: ------------------------------------------------------------------------</div><div>[0]PETSC ERROR: Petsc Development HG revision: f25e2c7d799b418675085efb3f389141e0dd9d6e  HG Date: Sat Oct 01 17:17:01 2011 -0500</div>
<div>[0]PETSC ERROR: See docs/changes/index.html for recent updates.</div><div>[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.</div><div>[0]PETSC ERROR: See docs/index.html for manual pages.</div><div>
[0]PETSC ERROR: ------------------------------------------------------------------------</div><div>[0]PETSC ERROR: ./ex26 on a ompi named batura by jed Tue Oct 11 22:25:15 2011</div><div>[0]PETSC ERROR: Libraries linked from /home/jed/petsc/ompi/lib</div>
<div>[0]PETSC ERROR: Configure run at Sun Oct  2 10:41:14 2011</div><div>[0]PETSC ERROR: Configure options --download-blacs --download-hypre --download-ml --download-mumps --download-parmetis --download-scalapack --download-spai --download-sundials --download-superlu --download-superlu_dist --download-umfpack --with-parmetis-dir=/usr --with-rosw --with-shared-libraries --with-single-library=0 --with-sowing PETSC_ARCH=ompi --download-parms</div>
<div>[0]PETSC ERROR: ------------------------------------------------------------------------</div><div>[0]PETSC ERROR: MatFDColoringApply() line 515 in /home/jed/petsc/src/mat/matfd/fdmatrix.c</div><div>[0]PETSC ERROR: DMComputeJacobianDefault() line 813 in /home/jed/petsc/src/dm/interface/dm.c</div>
<div>[0]PETSC ERROR: DMComputeJacobian() line 1290 in /home/jed/petsc/src/dm/interface/dm.c</div><div>[0]PETSC ERROR: SNESDMComputeJacobian() line 23 in /home/jed/petsc/src/snes/interface/snes.c</div><div>[0]PETSC ERROR: SNESComputeJacobian() line 1410 in /home/jed/petsc/src/snes/interface/snes.c</div>
<div>[0]PETSC ERROR: SNESSolve_KSPONLY() line 32 in /home/jed/petsc/src/snes/impls/ksponly/ksponly.c</div><div>[0]PETSC ERROR: SNESSolve() line 2793 in /home/jed/petsc/src/snes/interface/snes.c</div><div>[0]PETSC ERROR: TSStep_RosW() line 669 in /home/jed/petsc/src/ts/impls/rosw/rosw.c</div>
<div>[0]PETSC ERROR: TSStep() line 1813 in /home/jed/petsc/src/ts/interface/ts.c</div><div>[0]PETSC ERROR: TSSolve() line 1929 in /home/jed/petsc/src/ts/interface/ts.c</div><div>[0]PETSC ERROR: main() line 164 in src/ts/examples/tutorials/ex26.c</div>
</div><div><br></div><div><br></div><div><br></div></div></div></div>