36    void update(
const DiffCoeffFunc& computeDiffCoeff)
 
   38        for (
int phaseIdx = 0; phaseIdx < numPhases; ++ phaseIdx)
 
   39            for (
int compIIdx = 0; compIIdx < numComponents; ++compIIdx)
 
   40                for (
int compJIdx = compIIdx+1; compJIdx < numComponents; ++compJIdx)
 
   41                    diffCoeff_[getIndex_(phaseIdx, compIIdx, compJIdx)]
 
   42                        = computeDiffCoeff(phaseIdx, compIIdx, compJIdx);