[petsc-users] Solution output of TS

[Michael E Henderson]

Hi,

I'm using TS and seeing a ":flickering" in the output. I use

  ierr=TSCreate(PETSC_COMM_WORLD,&timeStepper);
  ierr=TSSetProblemType(timeStepper,TS_NONLINEAR); 
  ierr=TSSetFromOptions(timeStepper); 
  ierr=TSSetIFunction(timeStepper,     formOperatorImplicitTimeFunction , 
(void*)Ldata);
  ierr=TSSetIJacobian(timeStepper,A,A, formOperatorImplicitTimeJacobian , 
(void*)Ldata);
  ierr=TSMonitorSet(timeStepper,SaveTSSolution,(void*)(&data),NULL); 

I'm using the implicit nonlinear formulation for an ADE, and writing the 
solution in the monitoring routine.

If I instead write out the solution in the IFunction routine when the 
residual is small I see a different (better) solution passed in.

Does this sound like a problem that's been seen before?

Thanks,

MIke Henderson
------------------------------------------------------------------------------------------------------------------------------------
Mathematical Sciences, TJ Watson Research Center
mhender at watson.ibm.com
http://www.research.ibm.com/people/h/henderson/
http://multifario.sourceforge.net/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20101213/d6fbb8e4/attachment.htm>