b2right_stretch_tensor.H

Go to the documentation of this file.

//---------------------------------------------------------------------------
// b2right_stretch_tensor.H --
//
//     Compute the right stretch tensor and its derivatives.
//
// written by Thomas Ludwig.
//
// Copyright (c) 2011-2012 SMR Engineering & Development SA
//               2502 Bienne, Switzerland
//
// All Rights Reserved.  Proprietary source code.  The contents of
// this file may not be disclosed to third parties, copied or
// duplicated in any form, in whole or in part, without the prior
// written permission of SMR.
//---------------------------------------------------------------------------
 
#ifndef __B2_RIGHT_STRETCH_TENSOR_H__
#define __B2_RIGHT_STRETCH_TENSOR_H__
 
#include <algorithm>
#include <cmath>
 
#include "b2singular_value_decomposition.H"
#include "b2tensor_calculus.H"
 
namespace b2000 {
 
inline void polar_decomposition(const double F[3][3], double R[3][3], double U[3][3]) {
    double u[3][3];
    double s[3];
    double v[3][3];
 
    svd_3_3(F, u, s, v);
    inner_product_3_3_NN(u, v, R);
    inner_product_3_3_TN(R, F, U);
}
 
class RightStretchTensor_3_3 {
public:
    RightStretchTensor_3_3() {}
 
    RightStretchTensor_3_3(const double F[3][3], const bool second = false) { compute(F, second); }
 
    void compute(const double F[3][3], const bool second = false);
 
    void compute_numerically(const double F[3][3], const bool second);
 
    static bool check(const double F[3][3], const bool verbose = true);
 
private:
    double u[3][3];
    double s[3];
    double v[3][3];
 
    double Y[6];
    double DY[6];
 
    double YU[6];
    double DYU[6];
 
    double RY[9];
    double DRY[9];
 
    double RYU[9];
    double DRYU[9];
 
public:
    double R_3_3[3][3];
 
    double U_3_3[3][3];
 
    double R[9];
 
    double U[6];
 
    double DU[6][9];
 
    double DR[9][9];
 
    double DDU[6][9][9];
};
}  // namespace b2000
 
#endif /* __B2_RIGHT_STRETCH_TENSOR_H__ */