anonymous_namespace{meshUtil.cpp} Namespace Reference

Functions
util::SymMatrix3	strainFromB (const fe::B &B, const std::vector< util::Point > &u)

Function Documentation

strainFromB()

util::SymMatrix3 anonymous_namespace{meshUtil.cpp}::strainFromB	(	const fe::B &	B,
		const std::vector< util::Point > &	u
	)

Definition at line 31 of file meshUtil.cpp.

                                                            {
  std::vector<double> uflat(B.nDof());
  const int dim = B.dim();
  for (size_t a = 0; a < u.size(); a++) {
    uflat[dim * static_cast<int>(a)] = u[a][0];
    if (dim > 1)
      uflat[dim * static_cast<int>(a) + 1] = u[a][1];
    if (dim > 2)
      uflat[dim * static_cast<int>(a) + 2] = u[a][2];
  }
  std::vector<double> e(B.nStrain(), 0.);
  for (int i = 0; i < B.nStrain(); i++)
    for (int j = 0; j < B.nDof(); j++)
      e[i] += B(i, j) * uflat[j];
  util::SymMatrix3 s;
  s(0, 0) = e[0];
  if (dim > 1) {
    s(1, 1) = e[1];
    s(0, 1) = (dim == 2) ? e[2] : e[5];
  }
  if (dim > 2) {
    s(2, 2) = e[2];
    s(1, 2) = e[3];
    s(0, 2) = e[4];
  }
  return s;
}

References fe::B::dim(), fe::B::nDof(), and fe::B::nStrain().

Here is the call graph for this function: