Method BilinearForm

comp/bilinearform.cpp:88–144

  BilinearForm :: 
  BilinearForm (shared_ptr<FESpace> afespace,
                const string & aname,
                const Flags & flags)
    : NGS_Object(afespace->GetMeshAccess(), flags, aname), fespace(afespace)
  {
    fespace2 = NULL;

    multilevel = true;
    symmetric = flags.GetDefineFlag ("symmetric");

    // low_order_bilinear_form = NULL;
    linearform = NULL;


    SetGalerkin( flags.GetDefineFlag( "project" ) );
    SetNonAssemble (flags.GetDefineFlag ("nonassemble"));
    SetDiagonal (flags.GetDefineFlag ("diagonal"));
    if (flags.GetDefineFlag ("nonsym"))  SetSymmetric (0);
    if (flags.GetDefineFlag ("nonmultilevel")) SetMultiLevel (0);
    SetHermitean (flags.GetDefineFlag ("hermitean"));
    SetHermitean (flags.GetDefineFlag ("hermitian"));
    SetUnusedDiag (flags.GetNumFlag ("unuseddiag",0.0));
    SetEpsRegularization (flags.GetNumFlag ("regularization",0));
  
    SetPrint (flags.GetDefineFlag ("print"));
    SetPrintElmat (flags.GetDefineFlag ("printelmat"));
    SetElmatEigenValues (flags.GetDefineFlag ("elmatev")); 
    SetTiming (flags.GetDefineFlag ("timing"));
    SetEliminateInternal (flags.GetDefineFlag ("eliminate_internal") || flags.GetDefineFlag("condense"));

    
    bool has_private_space = false;
    fespace->TraverseTree([&has_private_space](const FESpace& fes) {
      if (dynamic_cast<const HiddenFESpace*>(&fes))
        has_private_space = true;
    });
    SetEliminateHidden(has_private_space);
    if (!flags.GetDefineFlagX("eliminate_hidden").IsMaybe())
      SetEliminateHidden (flags.GetDefineFlag ("eliminate_hidden"));

    
    SetKeepInternal (eliminate_internal &&
                     !flags.GetDefineFlagX ("keep_internal").IsFalse() &&
                     !flags.GetDefineFlag ("nokeep_internal"));
    SetStoreInner (flags.GetDefineFlag ("store_inner"));
    precompute = flags.GetDefineFlag ("precompute");
    checksum = flags.GetDefineFlag ("checksum");
    spd = flags.GetDefineFlag ("spd");
    geom_free = flags.GetDefineFlag("geom_free");
    matrix_free_bdb = flags.GetDefineFlag("matrix_free_bdb");    
    nonlinear_matrix_free_bdb = flags.GetDefineFlag("nonlinear_matrix_free_bdb");    
    if (spd) symmetric = true;
    SetCheckUnused (!flags.GetDefineFlagX("check_unused").IsFalse());
    if (flags.NumFlagDefined("delete_zero_elements"))
      delete_zero_elements = flags.GetNumFlag("delete_zero_elements", 0.0);
  }