MinimalSurfaceEquation::MinimalSurfaceElement::OpAssembleTangent Struct Reference

Assemble tangent matrix. More...

#include <users_modules/minimal_surface_equation/src/MinimalSurfaceElement.hpp>

Inheritance diagram for MinimalSurfaceEquation::MinimalSurfaceElement::OpAssembleTangent:

Collaboration diagram for MinimalSurfaceEquation::MinimalSurfaceElement::OpAssembleTangent:

Public Member Functions
	OpAssembleTangent (const string field_name, CommonData &common_data)
PetscErrorCode	doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, DataForcesAndSourcesCore::EntData &row_data, DataForcesAndSourcesCore::EntData &col_data)

Public Attributes
CommonData &	commonData
MatrixDouble	kTangent
VectorDouble	auxVec
MatrixDouble	auxMat

Detailed Description

Assemble tangent matrix.

Examples

minimal_surface_area.cpp.

Definition at line 328 of file MinimalSurfaceElement.hpp.

Constructor & Destructor Documentation

OpAssembleTangent()

MinimalSurfaceEquation::MinimalSurfaceElement::OpAssembleTangent::OpAssembleTangent ( const string  field_name, CommonData &  common_data  )

inline

Definition at line 331 of file MinimalSurfaceElement.hpp.

        : FaceElementForcesAndSourcesCore::UserDataOperator(
              field_name, UserDataOperator::OPROWCOL),
          commonData(common_data) {
      sYmm = false; // matrix is not symmetric
    }

Member Function Documentation

doWork()

PetscErrorCode MinimalSurfaceEquation::MinimalSurfaceElement::OpAssembleTangent::doWork ( int  row_side, int  col_side, EntityType  row_type, EntityType  col_type, DataForcesAndSourcesCore::EntData &  row_data, DataForcesAndSourcesCore::EntData &  col_data  )

inline

Definition at line 340 of file MinimalSurfaceElement.hpp.

                                                                     {
      MoFEMFunctionBegin;
 
      int row_nb_dofs = row_data.getIndices().size();
      if (row_nb_dofs == 0) {
        MoFEMFunctionReturnHot(0);
      }
      int col_nb_dofs = col_data.getIndices().size();
      if (col_nb_dofs == 0) {
        MoFEMFunctionReturnHot(0);
      }
      kTangent.resize(row_nb_dofs, col_nb_dofs, false);
      kTangent.clear();
      auxVec.resize(col_nb_dofs, false);
      auxMat.resize(col_nb_dofs, 2, false);
      int nb_gauss_pts = row_data.getN().size1();
      for (int gg = 0; gg != nb_gauss_pts; gg++) {
        // cerr << "gg " << gg << endl;
        double val = getGaussPts()(2, gg) * getArea();
        MatrixAdaptor row_diff_n = row_data.getDiffN(gg);
        MatrixAdaptor col_diff_n = col_data.getDiffN(gg);
        double an = commonData.aN[gg];
        ublas::matrix_row<MatrixDouble> grad_u(commonData.gradU, gg);
        ublas::matrix_row<MatrixDouble> an_pow3_by_grad_u(
            commonData.aNpow3byGradU, gg);
        // cerr << "grad_u " << grad_u << endl;
        // cerr << "an_pow3_by_grad_u " << an_pow3_by_grad_u << endl;
        noalias(auxVec) = prod(col_diff_n, an_pow3_by_grad_u);
        // cerr << "auxVec " << auxVec << endl;
        noalias(auxMat) = outer_prod(auxVec, grad_u);
        // cerr << "auxMat " << auxMat << endl;
        // cerr << row_diff_n << endl;
        // cerr << col_diff_n << endl;
        // cerr << prod(row_diff_n,trans(val*(an*col_diff_n+auxMat))) << endl;
        noalias(kTangent) +=
            val * prod(row_diff_n, trans(an * col_diff_n - auxMat));
        // cerr << "kTangent " << kTangent << endl;
      }
      CHKERR MatSetValues(getFEMethod()->snes_B, row_nb_dofs,
                          &*row_data.getIndices().data().begin(), col_nb_dofs,
                          &*col_data.getIndices().data().begin(),
                          &*kTangent.data().begin(), ADD_VALUES);
      MoFEMFunctionReturn(0);
    }

Member Data Documentation

auxMat

MatrixDouble MinimalSurfaceEquation::MinimalSurfaceElement::OpAssembleTangent::auxMat

Definition at line 339 of file MinimalSurfaceElement.hpp.

auxVec

VectorDouble MinimalSurfaceEquation::MinimalSurfaceElement::OpAssembleTangent::auxVec

Definition at line 338 of file MinimalSurfaceElement.hpp.

commonData

CommonData& MinimalSurfaceEquation::MinimalSurfaceElement::OpAssembleTangent::commonData

Definition at line 330 of file MinimalSurfaceElement.hpp.

kTangent

MatrixDouble MinimalSurfaceEquation::MinimalSurfaceElement::OpAssembleTangent::kTangent

Definition at line 337 of file MinimalSurfaceElement.hpp.

---

The documentation for this struct was generated from the following file:

• users_modules/minimal_surface_equation/src/MinimalSurfaceElement.hpp