FFaTensorTransforms.H

Go to the documentation of this file.

// SPDX-FileCopyrightText: 2023 SAP SE
//
// SPDX-License-Identifier: Apache-2.0
//
// This file is part of FEDEM - https://openfedem.org
 
#ifndef FFA_TENSOR_TRANSFORMS_H
#define FFA_TENSOR_TRANSFORMS_H
 
 
namespace FFaTensorTransforms 
{
  double vonMises(int N, const double* S);
  double vonMises(double s11, double s22, double s12);
  double vonMises(double s11, double s22, double s33,
                  double s12, double s13, double s23);
 
  int principalDirs(int N, const double* S, double* pVals, double* pDirs);
  int principalDirs(const double* S,
                    double* pVal, double* p1Dir, double* p2Dir);
  int principalDirs(const double* S,
                    double* pVal, double* p1Dir, double* p2Dir, double* p3Dir);
 
  bool principalValues(int N, const double* S, double* P);
  bool principalVals2D(double s11, double s22, double s12, double* pVal);
  bool principalVals3D(double s11, double s22, double s33,
                       double s12, double s13, double s23, double* pVal);
 
  // For backward compatibility, purge these two later..
  inline bool principalValues(double s11, double s22, double s12,
                              double& p1, double& p2)
  {
    double p[3] = { 0.0, 0.0, 0.0 };
    bool ok = principalVals2D(s11,s22,s12,p);
    p1 = p[0];
    p2 = p[1];
    return ok;
  }
  inline bool principalValues(double s11, double s22, double s33,
                              double s12, double s13, double s23,
                              double& p1, double& p2, double& p3)
  {
    double p[3] = { 0.0, 0.0, 0.0 };
    bool ok = principalVals3D(s11,s22,s33,s12,s13,s23,p);
    p1 = p[0];
    p2 = p[1];
    p3 = p[2];
    return ok;
  }
 
  double maxShearValue(double pMax, double pMin);
  void   maxShearDir  (int N, const double* pMaxDir, const double* pMinDir,
                       double* maxShearDir);
 
  void rotate2D(const double* inTensor,
                const double* rotMx,
                double* outTensor);
  void rotate  (const double* inTensor,
                const double* eX, const double* eY,
                double* outTensor);
  void rotate3D(const double* inTensor,
                const double* rotMx,
                double* outTensor);
  void rotate  (const double* inTensor,
                const double* eX, const double* eY, const double* eZ,
                double* outTensor);
 
  void from2Dto3D(const double* tensor2D, double* tensor3D);
  void from3Dto2D(const double* tensor3D, double* tensor2D);
}
 
#endif