Method Update

comp/contact.cpp:985–1227

  }

  void ContactBoundary::
  Update(shared_ptr<GridFunction> displacement_,
         shared_ptr<BilinearForm> bf,
         int intorder, double h, bool both_sides)
  {
    if(!displacement_ && !bf)
      throw Exception("Either displacement or BilinearForm needed in ContactBoundary update!");
    if(!trial_fes)
      {
        if(bf)
          {
            trial_fes = bf->GetTrialSpace();
            test_fes = bf->GetTestSpace();
          }
        else
          {
            trial_fes = displacement_->GetFESpace();
            test_fes = displacement_->GetFESpace();
          }
      }
    if(bf && (bf->GetTrialSpace().get() != trial_fes.get()))
      throw Exception("BilinearForm on different space as given to ContactBoundary!");
    if(displacement_)
      fes_displacement = displacement_->GetFESpace();

    if(draw_pairs)
    {
      secondary_points.SetSize(0);
      primary_points.SetSize(0);
    }

    shared_ptr<GridFunction> displacement = nullptr;
    if(displacement_)
      {
        auto flags = displacement_->GetFlags();
        flags.SetFlag("novisual");
        displacement = CreateGridFunction(displacement_->GetFESpace(), "_cb_displacement", flags);
        displacement->Update();
        displacement->GetVector() = displacement_->GetVector();
      }
    if (displacement)
      gap->Update(displacement, 10*displacement->GetFESpace()->GetOrder(), h,
                  both_sides);
    else
      gap->Update(nullptr, 10, h, both_sides);
      
    auto mesh = trial_fes->GetMeshAccess();
    if(mesh->GetDimension() == 2)
      static_pointer_cast<DisplacedNormal<2>>(normal)->Update(displacement);
    else
      static_pointer_cast<DisplacedNormal<3>>(normal)->Update(displacement);
    if(bf)
      {
        static Timer t1("Build contact pairs");
        RegionTimer regt1(t1);
        LocalHeap lh(1000000, "ContactBoundary-Update", true);
        Iterate<2>
          ([&](auto i)