Method BDDCMatrix

comp/bddc.cpp:52–256

    void SetCoarseGridPreconditioner (bool ah = true) { coarse = ah; }
    
    BDDCMatrix (shared_ptr<BilinearForm> abfa,
                shared_ptr<BitArray> freedofs,
                Flags flags,
		const string & ainversetype, 
		const string & acoarsetype, 
                bool ablock, 
                bool ahypre
      )
      : bfa(abfa), block(ablock), inversetype(ainversetype), coarsetype(acoarsetype)
    {
      static Timer timer ("BDDC Constructor");

      fes = bfa->GetFESpace();
      
      coarse = (coarsetype != "none");

      hypre = ahypre;

      local = flags.GetDefineFlag("local");
      if (block)
        blocktype = flags.GetStringFlag ("blocktype", "");
      // pwbmat = NULL;
      inv = NULL;

      inv_coarse = NULL;
      // tmp = NULL;
      // tmp2 = NULL;
      RegionTimer reg(timer);

      // auto fes = bfa -> GetFESpace();
      shared_ptr<MeshAccess> ma = fes->GetMeshAccess();

      Array<int> wbdcnt(ma->GetNE()+ma->GetNSE()+ma->GetNCD2E()+bfa->GetSpecialElements().Size());
      Array<int> ifcnt(ma->GetNE()+ma->GetNSE()+ma->GetNCD2E()+bfa->GetSpecialElements().Size());
      wbdcnt = 0;
      ifcnt = 0;
      // const BitArray & freedofs = *fes->GetFreeDofs();
      

      LocalHeap lh(10000, "BDDC-constr, dummy heap");
      
      for (auto vb : { VOL, BND, BBND })
        IterateElements 
          (*fes, vb, lh, 
           [&] (FESpace::Element el, LocalHeap & lh)
           {
             int base = (vb == VOL) ? 0 : ((vb == BND) ? ma->GetNE() : ma->GetNE() + ma->GetNSE());
             for (auto d : el.GetDofs())
               {
                 if (!IsRegularDof(d)) continue;
                 if (!freedofs->Test(d)) continue;
                 COUPLING_TYPE ct = fes->GetDofCouplingType(d);
                 if ( ((ct & CONDENSABLE_DOF) != 0) && bfa->UsesEliminateInternal()) continue;
		
                 int ii = base + el.Nr();
                 if (ct == WIREBASKET_DOF)
                   wbdcnt[ii]++;
                 else