Function TransposeCF

fem/coefficient.cpp:4350–4367

  }

  shared_ptr<CoefficientFunction> TransposeCF (shared_ptr<CoefficientFunction> coef)
  {
    if (coef->IsZeroCF())
      {
        auto dims = coef->Dimensions();
        if (dims.Size() == 2)
            coef->SetDimensions( Array( {dims[1], dims[0]}) );
        else
            throw Exception("Transpose of non-matrix called");
        return coef;
      }

    // if (coef.use_count() == 1)
    if (dynamic_pointer_cast<IdentityCoefficientFunction> (coef) && !coef->IsVariable())
      return coef;

    return make_shared<TransposeCoefficientFunction> (coef);
  }