Method GenerateCode

fem/symbolicintegrator.cpp:216–347

  }

  void ProxyFunction ::
  GenerateCode(Code &code, FlatArray<int> inputs, int index) const
  {
    // auto dims = Dimensions();
    
    string header = "\n\
    {flatmatrix} {values};\n\
    ProxyUserData * {ud} = (ProxyUserData*)mir.GetTransformation().userdata;\n\
    {\n\
      // if (!{ud})\n                                                      \
      // throw Exception (\"cannot evaluate ProxyFunction without userdata\");\n \
          ";

    if(!testfunction) {
      header+=
"      if ({ud}->fel) {\n\
     //     if ({ud}->HasMemory ({this})) {\n";
      if(code.is_simd) {
        header += "auto x = {ud}->GetAMemory ({this});\n";
        header += "{values}.AssignMemory(x.Height(), x.Width(), &x(0,0));\n";
      } else {
        header += "auto x = {ud}->GetMemory ({this});\n";
        header += "{values}.AssignMemory(x.Height(), x.Width(), &x(0,0));\n";
      }
      header+=
"      //    }\n\
       //   else\n\
       //     throw Exception(\"userdata has no memory!\");\n\
      }\n";
    }
    header += "}\n";
    if(code.deriv || !testfunction)
      header += "const bool {has_values} = {ud}->HasMemory({this});\n";
    else
      header += "constexpr bool {has_values} = true;\n"; // always need values in case code.deriv==0 (i.e. evaluation of proxy)

    for (int i = 0; i < this->Dimension(); i++) {
      header += Var("comp", index,i,this->Dimensions()).Declare("{scal_type}", 0.0);
        if(!testfunction && code.deriv==2)
          {
          header += "if(( ({ud}->trialfunction == {this}) && ({ud}->trial_comp=="+ToLiteral(i)+"))\n"+
          " ||  (({ud}->testfunction == {this}) && ({ud}->test_comp=="+ToLiteral(i)+")))\n";
        }
        else
          header += "if({ud}->{comp_string}=="+ToLiteral(i)+" && {ud}->{func_string} == {this})\n";
        header += Var("comp", index,i,this->Dimensions()).S() + string("{get_component}") + " = 1.0;\n";
    }
    
    string body = "";

    if (code_uses_tensors)
      {
        body += "Tens<" + code.res_type;
        for (auto d : this->Dimensions())
          body += ',' + ToLiteral(d);
        body += "> var_" + ToLiteral(index) + ";\n";
      }