Method ContactIntegrator

comp/contact.cpp:671–727

  }

  ContactIntegrator::ContactIntegrator(shared_ptr<CoefficientFunction> _cf,
                                       bool _deformed)
    : cf(_cf), deformed(_deformed)
  {
    cf->TraverseTree
      ([&](CoefficientFunction& nodecf)
       {
         auto proxy = dynamic_cast<ProxyFunction*>(&nodecf);
         if (proxy && !proxy->IsTestFunction() && !trial_proxies.Contains(proxy))
           trial_proxies.Append(proxy);

         if(proxy && proxy->IsTestFunction() && !test_proxies.Contains(proxy))
           test_proxies.Append(proxy);

         auto gf = dynamic_cast<GridFunction*>(&nodecf);
         if (gf && !cf_gridfunctions.Contains(&nodecf))
           cf_gridfunctions.Append(&nodecf);
       });
    if(trial_proxies.Size())
      fes = trial_proxies[0]->GetFESpace();
    else
      {
        if (test_proxies.Size())
          fes = test_proxies[0]->GetFESpace();
        else
          throw Exception("No trial or test function found in ContactIntegrator");
      }

    dcf_dtest.SetSize(test_proxies.Size());
    ddcf_dtest_dtrial.SetSize(test_proxies.Size(), trial_proxies.Size());
    for (int i = 0; i < test_proxies.Size(); i++)
      {
        try
          {
            CoefficientFunction::T_DJC cache;
            dcf_dtest[i] = cf->DiffJacobi(test_proxies[i], cache);
          }
        catch (const Exception& e)
          {
            cout << IM(5) << "dcf_dtest has thrown exception " << e.What() << endl;
          }
          for (auto j : Range(trial_proxies))
            {
              if (dcf_dtest[i])
                try
                  {
                    CoefficientFunction::T_DJC cache2;
                      ddcf_dtest_dtrial(i, j) = dcf_dtest[i]->DiffJacobi(trial_proxies[j], cache2);
                  }
                catch (const Exception& e)
                  {
                    cout << IM(1) << "ddcf_dtest_dtrial has thrown exception " << e.What() << endl;
                  }
            }
      }

  }