void FD1Solver::compute_localSpeed() { SField lowerSpeed; SField upperSpeed; for(int dir=0;dir<m_n;dir++) { upperSpeed = m_eq->get_max_eigenvalue(upper_right_intermediate_un_values[dir], dir); lowerSpeed = m_eq->get_max_eigenvalue(lower_right_intermediate_un_values[dir], dir); right_localSpeed[dir] = upperSpeed.max_field(lowerSpeed); upperSpeed = m_eq->get_max_eigenvalue(upper_left_intermediate_un_values[dir], dir); lowerSpeed = m_eq->get_max_eigenvalue(lower_left_intermediate_un_values[dir], dir); left_localSpeed[dir] = upperSpeed.max_field(lowerSpeed); } }