예제 #1
0
파일: FD1Solver.cpp 프로젝트: Yukee/flume
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);
    }
}