RotationTensorTest.C

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#include "gtest/gtest.h"

#include "RankFourTensor.h"
#include "RotationTensor.h"

TEST(RotationTensorTest, rotation)
{
  std::vector<Real> input(81);
  for (size_t i = 0; i < input.size(); ++i)
    input[i] = i + 1.1;

  RankFourTensor a(input, RankFourTensor::general);

  RotationTensor xrot(RotationTensor::XAXIS, 90);
  RotationTensor yrot(RotationTensor::YAXIS, 90);

  auto xyrot = xrot * yrot;
  auto yxrot = yrot * xrot;

  auto axy1 = a;
  axy1.rotate(xrot);
  axy1.rotate(yrot);

  auto axy2 = a;
  axy2.rotate(yxrot);

  EXPECT_NEAR(0, (axy1 - axy2).L2norm(), 1E-8);
}