985 {
987
988 const double vol = getMeasure();
989 auto t_w = getFTensor0IntegrationWeight();
990 auto t_coords = getFTensor1CoordsAtGaussPts();
991 auto t_row_base = row_data.getFTensor0N();
992 auto t_row_diff_base = row_data.getFTensor1DiffN<
SPACE_DIM>();
993
994#ifndef NDEBUG
995 if (row_data.getDiffN().size1() != row_data.getN().size1())
997 if (row_data.getDiffN().size2() != row_data.getN().size2() *
SPACE_DIM) {
999 << "Side " << rowSide << " " << CN::EntityTypeName(rowType);
1000 MOFEM_LOG(
"SELF", Sev::error) << row_data.getN();
1001 MOFEM_LOG(
"SELF", Sev::error) << row_data.getDiffN();
1003 }
1004
1005 if (col_data.getDiffN().size1() != col_data.getN().size1())
1007 if (col_data.getDiffN().size2() != col_data.getN().size2() *
SPACE_DIM) {
1009 << "Side " << rowSide << " " << CN::EntityTypeName(rowType);
1010 MOFEM_LOG(
"SELF", Sev::error) << col_data.getN();
1011 MOFEM_LOG(
"SELF", Sev::error) << col_data.getDiffN();
1013 }
1014#endif
1015
1017
1018 auto t_h = getFTensor0FromVec(*
hPtr);
1019 auto t_grad_g = getFTensor1FromMat<SPACE_DIM>(*
gradGPtr);
1020
1021 for (int gg = 0; gg != nbIntegrationPts; gg++) {
1022
1023 const double r = t_coords(0);
1025
1026 int rr = 0;
1027 for (; rr != nbRows; ++rr) {
1028
1029 auto t_col_base = col_data.getFTensor0N(gg, 0);
1030 auto t_col_diff_base = col_data.getFTensor1DiffN<
SPACE_DIM>(gg, 0);
1031
1032 for (int cc = 0; cc != nbCols; ++cc) {
1033
1034 locMat(rr, cc) += (t_row_base * t_col_base * alpha);
1035
1036 ++t_col_base;
1037 ++t_col_diff_base;
1038 }
1039
1040 ++t_row_base;
1041 ++t_row_diff_base;
1042 }
1043
1044 for (; rr < nbRowBaseFunctions; ++rr) {
1045 ++t_row_base;
1046 ++t_row_diff_base;
1047 }
1048
1049 ++t_h;
1050 ++t_grad_g;
1051 ++t_w;
1052 ++t_coords;
1053 }
1054
1055 } else {
1056
1057 auto t_h = getFTensor0FromVec(*
hPtr);
1058 auto t_grad_g = getFTensor1FromMat<SPACE_DIM>(*
gradGPtr);
1059 auto t_u = getFTensor1FromMat<U_FIELD_DIM>(*
uPtr);
1060
1061 auto ts_a = getTSa();
1062
1063 for (int gg = 0; gg != nbIntegrationPts; gg++) {
1064
1065 const double r = t_coords(0);
1067
1069 int rr = 0;
1070 for (; rr != nbRows; ++rr) {
1071
1072 auto t_col_base = col_data.getFTensor0N(gg, 0);
1073 auto t_col_diff_base = col_data.getFTensor1DiffN<
SPACE_DIM>(gg, 0);
1074
1075 for (int cc = 0; cc != nbCols; ++cc) {
1076
1077 locMat(rr, cc) += (t_row_base * t_col_base * alpha) * ts_a;
1078 locMat(rr, cc) +=
1079 (t_row_base * alpha) * (t_col_diff_base(
i) * t_u(
i));
1080 locMat(rr, cc) +=
1081 (t_row_diff_base(
i) * t_grad_g(
i)) * (t_col_base * m_dh);
1082
1083 ++t_col_base;
1084 ++t_col_diff_base;
1085 }
1086
1087 ++t_row_base;
1088 ++t_row_diff_base;
1089 }
1090
1091 for (; rr < nbRowBaseFunctions; ++rr) {
1092 ++t_row_base;
1093 ++t_row_diff_base;
1094 }
1095
1096 ++t_u;
1097 ++t_h;
1098 ++t_grad_g;
1099 ++t_w;
1100 ++t_coords;
1101 }
1102 }
1103
1105 }
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
@ MOFEM_DATA_INCONSISTENCY
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
#define MOFEM_LOG(channel, severity)
Log.
FTensor::Index< 'i', SPACE_DIM > i
constexpr IntegrationType I