Method computeH

SRC/element/XMUelements/AC3D8HexWithSensitivity.cpp:967–1007

}

int
AC3D8HexWithSensitivity::computeH(void)
{
  if (H == 0 || DH == 0) {
    H = new Matrix*[numGP];
    DH = new Matrix*[numGP];
    if (H == 0 || DH == 0) {
      opserr << "AC3D8HexWithSensitivity::computeH - out of memory!\n";
      return -3;
    }

    double r = 0.0;
    double s = 0.0;
    double t = 0.0;
    short where = 0;
    
    for(short GP_c_r = 1; GP_c_r <= r_integration_order; GP_c_r++) {
      r = get_Gauss_p_c(r_integration_order, GP_c_r);
      for(short GP_c_s = 1; GP_c_s <= s_integration_order; GP_c_s++) {
        s = get_Gauss_p_c(s_integration_order, GP_c_s);
        for(short GP_c_t = 1; GP_c_t <= t_integration_order; GP_c_t++) {
          t = get_Gauss_p_c(t_integration_order, GP_c_t);
          
          H[where] = new Matrix(1, nodes_in_elem);
          DH[where] = new Matrix(dim, nodes_in_elem);
          if(H[where] == 0 || DH[where] == 0) {
            opserr << "AC3D8HexWithSensitivity::computeH - out of memory!\n";
            return -3;
          }
          
          *H[where] = interp_fun(r, s, t);
          *DH[where] = diff_interp_fun(r, s, t);
          
          where++;
        } // for GP_c_t
      } // for GP_c_s
    } // for GP_c_r
  }
  
  return 0;
}

int
AC3D8HexWithSensitivity::computeHH(void)