[MITgcm-devel] updating AD verification outputs

[Jean-Michel Campin]

Hi Patrick,

Getting also some small differences in tlm test, with 2 exp
falling below 13:
Y Y Y Y 16>12<16 FAIL  1D_ocean_ice_column
Y Y Y Y 16>12<16 FAIL  tutorial_tracer_adjsens.som81
(was all 16 last Sunday with TAF 2.3.8)

Did take a short look at the diffences between the code TAF 2.3.10
generates compared with 2.3.8, and I am seeing many more
paranthesis, in multiple multiplication expression, 
(a*b)*c versus a*b*c before, 
and also sum of more than 2 terms: (a+b)+c versus a+b+c before.
This can explain the small differences (we are still getting the same 
1rst 12 digits for any experiment).

Will hardly see any changes in gradient-check (since these diff 
are too small).
 
Probably going to update the output (ad & tlm) that are matching by
less than 13 digits.

Cheers,
Jean-Michel

On Sat, Mar 09, 2013 at 10:06:55AM -0500, Patrick Heimbach wrote:
> 
> Hi Jean-Michel,
> 
> before updating, let's first check the TLM to see if there are similar changes,
> and what the gradient checks indicate in terms of actual accuracy changes.
> 
> May have to send a bug report to Fastopt, which I am not sure I can do
> since I am being blocked from their taf-user at fastopt.de list.
> 
> Cheers
> -Patrick
> 
> On Mar 9, 2013, at 9:58 AM, Jean-Michel Campin <jmc at ocean.mit.edu> wrote:
> 
> > Hi Patrick,
> > 
> > Looks like the recent update of TAF version (from 2.3.8 until yesterday
> > to 2.3.10 today) change some output at the level of machine truncation:
> > 
> > On baudelaire with gfortran (no mpi):
> > Y Y Y Y 16>12<16 FAIL  1D_ocean_ice_column
> > Y Y Y Y 16>12<16 FAIL  bottom_ctrl_5x5
> > Y Y Y Y 16>12<16 FAIL  lab_sea.noseaicedyn
> > Y Y Y Y 16>12<16 FAIL  offline_exf_seaice
> > 
> > whereas yesterday the same 4 experiments:
> > Y Y Y Y 16>13<16 pass  1D_ocean_ice_column
> > Y Y Y Y 16>16<16 pass  bottom_ctrl_5x5
> > Y Y Y Y 16>16<16 pass  lab_sea.noseaicedyn
> > Y Y Y Y 16>16<16 pass  offline_exf_seaice
> > 
> > Should we update these 4 experiments AD output ?
> > 
> > Did not check yet the TLM output but will do later today.
> > 
> > Cheers,
> > Jean-Michel
> > 
> > _______________________________________________
> > MITgcm-devel mailing list
> > MITgcm-devel at mitgcm.org
> > http://mitgcm.org/mailman/listinfo/mitgcm-devel
> 
> 
> ---
> Patrick Heimbach | heimbach at mit.edu | http://www.mit.edu/~heimbach
> MIT | EAPS 54-1420 | 77 Massachusetts Ave | Cambridge MA 02139 USA
> FON +1-617-253-5259 | FAX +1-617-253-4464 | SKYPE patrick.heimbach
> 



> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel