[MITgcm-cvs] MITgcm/pkg/seaice CVS Commit
Martin Losch
mlosch at forge.csail.mit.edu
Thu Mar 7 03:32:03 EST 2013
Update of /u/gcmpack/MITgcm/pkg/seaice
In directory forge:/tmp/cvs-serv8005/pkg/seaice
Modified Files:
seaice_lsr.F
Log Message:
Make the subroutines seaice_lsr_tridiagu/v completely unreadable:
modify lsr-algorithm so that it vectorizes AND converges by using a
"zebra" method: split vector loop into two with step size 2; this
solves a general convergence problem for the JFNK solver on a vector
computer. It has the potential to even improve convergence for
non-vector machines, but it does change the results so it is turned
off by default. To turn on define CPP-flags SEAICE_VECTORIZE_LSR and
SEAICE_VECTORIZE_LSR_ZEBRA in SEAICE_OPTIONS.h
More information about the MITgcm-cvs
mailing list