void
IsotropicElasticityTensor::calculateEntries(unsigned int /*qp*/)
{
  calculateLameCoefficients();

  unsigned int i, j, k, l;
  i = j = k = l = 0;

  // Walk down the columns of the 9x9 matrix
  for (unsigned int q = 0; q < 81; ++q)
  {
    // This algorithm was developed by Derek Gaston and Cody Permann
    // it's based on page 29 of Michael Tonk's notes
    j += i / 3;
    k += j / 3;
    l += k / 3;

    i %= 3;
    j %= 3;
    k %= 3;
    l %= 3;

    _values[q] = isotropicEntry(i, j, k, l);

    ++i;
  }
}
void
SymmIsotropicElasticityTensor::calculateEntries(unsigned int /*qp*/)
{
  calculateLameCoefficients();

  const Real C12(_lambda);
  const Real C44(_mu);
  const Real C11(2*C44+C12);

  setEntries( C11, C12, C44 );
}