HcurlEdgeBase Struct Reference

Collaboration diagram for HcurlEdgeBase:

Public Member Functions
template<int DIM, bool CALCULATE_DIRVATIVES>
MoFEMErrorCode	calculate (int p, int nb_integration_pts, int n0_idx, int n1_idx, double n[], FTensor::Tensor1< double, 3 > t_grad_n[], FTensor::Tensor1< FTensor::PackPtr< double *, 3 >, 3 > &t_phi, FTensor::Tensor2< FTensor::PackPtr< double *, 3 *DIM >, 3, DIM > *t_diff_phi_ptr)

Public Attributes
FTensor::Index< 'i', 3 >	i
FTensor::Tensor1< double, 3 >	tGradN0pN1
FTensor::Tensor1< double, 3 >	tPhi0
FTensor::Tensor1< double, 3 >	tDiffb
VectorDouble	fI
MatrixDouble	diffFi

Detailed Description

Definition at line 1978 of file Hcurl.cpp.

Member Function Documentation

calculate()

template<int DIM, bool CALCULATE_DIRVATIVES>

MoFEMErrorCode HcurlEdgeBase::calculate ( int  p, int  nb_integration_pts, int  n0_idx, int  n1_idx, double  n[], FTensor::Tensor1< double, 3 >  t_grad_n[], FTensor::Tensor1< FTensor::PackPtr< double *, 3 >, 3 > &  t_phi, FTensor::Tensor2< FTensor::PackPtr< double *, 3 *DIM >, 3, DIM > *  t_diff_phi_ptr  )

inline

Definition at line 1991 of file Hcurl.cpp.

                                 {
 
    FTensor::Index<'j', DIM> j;
 
    MoFEMFunctionBegin;
 
    FTensor::Tensor1<double, 3> &t_grad_n0 = t_grad_n[n0_idx];
    FTensor::Tensor1<double, 3> &t_grad_n1 = t_grad_n[n1_idx];
    tGradN0pN1(i) = t_grad_n0(i) + t_grad_n1(i);
 
    fI.resize(p + 1);
    diffFi.resize(3, p + 1);
    diffFi.clear();
 
    FTensor::Tensor2<double, 3, DIM> t_diff_phi0;
 
    for (int gg = 0; gg != nb_integration_pts; ++gg) {
 
      const int shift_n = (DIM + 1) * gg;
      const double n0 = n[shift_n + n0_idx];
      const double n1 = n[shift_n + n1_idx];
 
      tPhi0(i) = n0 * t_grad_n1(i) - n1 * t_grad_n0(i);
      t_phi(i) = tPhi0(i);
 
      ++t_phi;
 
      if (CALCULATE_DIRVATIVES) {
 
        t_diff_phi0(i, j) =
            t_grad_n0(j) * t_grad_n1(i) - t_grad_n1(j) * t_grad_n0(i);
        (*t_diff_phi_ptr)(i, j) = t_diff_phi0(i, j);
        ++(*t_diff_phi_ptr);
      }
 
      if (p > 1) {
 
        if (CALCULATE_DIRVATIVES)
          CHKERR Jacobi_polynomials(p, 0, n1, n0 + n1, &t_grad_n1(0),
                                    &tGradN0pN1(0), &*fI.data().begin(),
                                    &*diffFi.data().begin(), DIM);
        else
          CHKERR Jacobi_polynomials(p, 0, n1, n0 + n1, nullptr, nullptr,
                                    &*fI.data().begin(), nullptr, DIM);
 
        FTensor::Tensor1<FTensor::PackPtr<double *, 1>, 3> t_diff_fi(
            &diffFi(0, 1), &diffFi(1, 1), &diffFi(2, 1));
 
        for (int oo = 1; oo <= p - 1; ++oo) {
 
          const double b = pow(n0 + n1, oo);
          t_phi(i) = b * fI[oo] * tPhi0(i);
 
          if (CALCULATE_DIRVATIVES) {
 
            tDiffb(i) =
                oo * pow(n0 + n1, oo - 1) * (t_grad_n0(i) + t_grad_n1(i));
            (*t_diff_phi_ptr)(i, j) = (b * fI[oo]) * t_diff_phi0(i, j) +
                                      (b * t_diff_fi(j)) * tPhi0(i) +
                                      tDiffb(j) * fI[oo] * tPhi0(i);
            ++t_diff_fi;
            ++(*t_diff_phi_ptr);
          }
 
          ++t_phi;
        }
      }
    }
 
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

diffFi

MatrixDouble HcurlEdgeBase::diffFi

Definition at line 1987 of file Hcurl.cpp.

fI

VectorDouble HcurlEdgeBase::fI

Definition at line 1986 of file Hcurl.cpp.

i

FTensor::Index<'i', 3> HcurlEdgeBase::i

Definition at line 1980 of file Hcurl.cpp.

tDiffb

FTensor::Tensor1<double, 3> HcurlEdgeBase::tDiffb

Definition at line 1984 of file Hcurl.cpp.

tGradN0pN1

FTensor::Tensor1<double, 3> HcurlEdgeBase::tGradN0pN1

Definition at line 1982 of file Hcurl.cpp.

tPhi0

FTensor::Tensor1<double, 3> HcurlEdgeBase::tPhi0

Definition at line 1983 of file Hcurl.cpp.

---

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

• src/approximation/impl/Hcurl.cpp