C++ (Cpp) SparseMatrix::MultTranspose Beispiele

Beispiel #1

Datei: fespace.cpp Projekt: YPCC/mfem

void FiniteElementSpace::UpdateAndInterpolate(int num_grid_fns, ...)
{
   if (mesh->GetState() == Mesh::NORMAL)
   {
      MFEM_ABORT("Mesh must be in two-level state, please call "
                 "Mesh::UseTwoLevelState before refining.");
   }

   FiniteElementSpace *cfes = SaveUpdate();

   // obtain the (transpose of) interpolation matrix between mesh levels
   SparseMatrix *R = GlobalRestrictionMatrix(cfes, 0);

   delete cfes;

   // interpolate the grid functions
   std::va_list vl;
   va_start(vl, num_grid_fns);
   for (int i = 0; i < num_grid_fns; i++)
   {
      GridFunction* gf = va_arg(vl, GridFunction*);
      if (gf->FESpace() != this)
      {
         MFEM_ABORT("Cannot interpolate: grid function is not based "
                    "on this space.");
      }

      Vector coarse_gf = *gf;
      gf->Update();
      R->MultTranspose(coarse_gf, *gf);
   }
   va_end(vl);

   delete R;
   mesh->SetState(Mesh::TWO_LEVEL_FINE);
}