Many thanks for your reply! The code I tried was very similar to the code in:
ts/examples/tutorials/ex7.c
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:
PetscViewerASCIIOpen(PETSC_COMM_WORLD,"snes.log",&rviewer);
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!
_________________________________________________________________
????MSN?????????????????????????
http://mobile.msn.com.cn/
