• Externalpackages: out into astrid-r1b
• Get a wiki together. OpenAD. Get ISSM website together ​https://trac.mcs.anl.gov/projects/ADIC/wiki/ADIC/ISSM
• Switch to metis 5.0

Hard stuff:

• wrap MatSetValues and MatAssembly into one routine.
• To do so, store Ke matrices into elements, do not add to global.
• Once all the Ke are built, loop through them and do all the MatSetValues at once, followed by MatAssembly.
• Do the same for VecSet and VecAssembly.
• Wrap all Petsc data types Matrix and Vector objects, and map all operations onto these new objects. Make this preprocessable for configuration using autotools. We will probably duplicate data, ex:

Matrix{
   ifdef  PETSC_WRAPPER
   Mat pmatrix;
   endif
   ifdef  DOUBLE_WRAPPER
   double* pmatrix;
   endif
   adouble* amatrix;
}

amatrix will be used in adolc. pmatrix will be used in Petsc routines. Matrix will have routines to map one into the other.

Ad:

Eric:

• Create ISSM namespace.
• Make Matrix Assembly one homogeneous block.
• Make Validation suite for Jean for his prognostic capabilities.

Jean:

• typedef all double to ISSM_DOUBLE
• Create an alloc layer with template functions for new and delete.
• Check that the new code is modular --without-kml, options are up to date.
• Replace memcpy, realloc with x layers. Replace all calls to malloc and free in the code.

Questions:

• set of active IssmDouble variables different from set of independent variables activated with <<=?
• same thing dependent variables?