Method addTimeIntegratedPointSourceNRF

src/Kernels/PointSourceClusterOnDevice.cpp:109–140

constexpr std::size_t MInvJInvPhisAtSourcesSpan = tensor::mInvJInvPhisAtSources::Shape[0];

void PointSourceClusterOnDevice::addTimeIntegratedPointSourceNRF(const std::array<real, 3>& slip,
                                                                 const real* mInvJInvPhisAtSources,
                                                                 const real* tensor,
                                                                 real a,
                                                                 const real* stiffnessTensor,
                                                                 double from,
                                                                 double to,
                                                                 real dofs[tensor::Q::size()]) {
  real rotatedSlip[3] = {real(0.0)};
  for (unsigned i = 0; i < 3; ++i) {
    for (unsigned j = 0; j < 3; ++j) {
      rotatedSlip[j] += tensor[j + i * 3] * slip[i];
    }
  }

  auto mom = [&](unsigned p, unsigned q) {
    real m = 0.0;
    for (unsigned i = 0; i < 3; ++i) {
      for (unsigned j = 0; j < 3; ++j) {
        m += -a * stiffnessTensor[p + 3 * q + 9 * i + 27 * j] * rotatedSlip[i] * tensor[6 + j];
      }
    }
    return m;
  };

  const real moment[6] = {mom(0, 0), mom(1, 1), mom(2, 2), mom(0, 1), mom(1, 2), mom(0, 2)};
  for (unsigned t = 0; t < 6; ++t) {
    for (unsigned k = 0; k < MInvJInvPhisAtSourcesSpan; ++k) {
      dofs[k + t * QSpan] += mInvJInvPhisAtSources[k] * moment[t];
    }
  }
}

void PointSourceClusterOnDevice::addTimeIntegratedPointSourceFSRM(real slip,
                                                                  const real* mInvJInvPhisAtSources,