Todo List

Class BimetricEvolve

EvolvedBy f, f_t --> add dissipation here?

Implement PIRK?

Global BimetricEvolve::determineGaugeFunctions (Int m)

Go through the order of evaluation for the gauge grid functions

GF_r and GF_rr requires fixed ghost zones

Global BimetricEvolve::integStep_CalcEvolutionRHS (Int m)

fixme: Severely lost accuracy in all second derivatives

Smoothen the time derivatives inside the left grid zone

Global BimetricEvolve::maximalSlice_2 (Int m, Int N, Real gAlp_at_N)

fixme: Should be k_g/2 in the above.

fixme: Current implementation uses rho + J of the perfect fluid (should use bimetric).

Global BimetricEvolve::maximalSlice_PostSteps (Int m, Int N)

fixme: add +TINY in the case if zero

Global gBet (m, n)

gBet, fBet to be grid functions (after calc of gAlp, fAlp)

Global GF_r (f)

A run-time CFDS_ORDER?

Global GridPoint::write (FILE *outf, bool isBinary, const std::vector< GF_Descriptor > output)

fld::output should be dynamic in ID class

Global MoL::quit (Int m=-1)

Distribute the abort msg to all ranks instead of calling MPI_Abort.

Class MPIWorld

Implement the collective I/O with MPI_File_open and MPI_File_write.