示例#1
0
void RBEIMConstruction::init_data()
{
  // set the coupling matrix to be diagonal
  _coupling_matrix.resize(n_vars());

  // resize zeroed the matrix; we just need to set diagonal entries
  for(unsigned int var=0; var<n_vars(); var++)
    _coupling_matrix(var,var) = 1;

  this->get_dof_map()._dof_coupling = &_coupling_matrix;

  Parent::init_data();
}
示例#2
0
void RBEIMConstruction::init_data()
{
  // set the coupling matrix to be diagonal
  _coupling_matrix.resize(n_vars());
  for(unsigned int var1=0; var1<n_vars(); var1++)
    for(unsigned int var2=0; var2<n_vars(); var2++)
    {
      unsigned char value = (var1==var2) ? 1 : 0;
      _coupling_matrix(var1,var2) = value;
    }

  this->get_dof_map()._dof_coupling = &_coupling_matrix;

  Parent::init_data();
}