About SNES solver of TS objects

tsjb00 tsjb00 at hotmail.com
Fri Sep 12 16:53:36 CDT 2008

Many thanks for your reply! The code I tried was very similar to the code in:
except that I changed the r.h.s. function.
The code related to TS:
    ierr = TSCreate(PETSC_COMM_WORLD,&ts);CHKERRQ(ierr);
    ierr = TSSetProblemType(ts,TS_NONLINEAR);CHKERRQ(ierr);
    ierr = TSSetRHSFunction(ts,FormFunction,&appctx);CHKERRQ(ierr);
    ierr = DAGetColoring(appctx.da,IS_COLORING_GLOBAL,&iscoloring);CHKERRQ(ierr);
    ierr = MatFDColoringCreate(J,iscoloring,&matfdcoloring);CHKERRQ(ierr);
    ierr = ISColoringDestroy(iscoloring);CHKERRQ(ierr);
    ierr = MatFDColoringSetFunction(matfdcoloring,(PetscErrorCode (*)(void))FormFunction,&appctx);CHKERRQ(ierr);
    ierr = MatFDColoringSetFromOptions(matfdcoloring);CHKERRQ(ierr);
    ierr = TSSetRHSJacobian(ts,J,J,TSDefaultComputeJacobianColor,matfdcoloring);CHKERRQ(ierr);
    dt = appctx.dtmin;
    ierr = TSSetInitialTimeStep(ts,0.0,dt);CHKERRQ(ierr);
     ierr = TSSetType(ts,TS_BEULER);CHKERRQ(ierr);
Then I tried to add the snes monitor:
   ierr = TSGetSNES(ts,&ts_snes);
   ierr = SNESMonitorSet(ts_snes,SNESMonitorDefault,&rviewer,PETSC_NULL);
I got the snes.log with long lines which started with countless blank/space and ended at the very end with the SNES iteration info.  As a result, the file took a lot of memory space even with one line of outputs.
Please let me know if I did something wrong.
Have a nice weekend!

More information about the petsc-users mailing list