adipose_tissue_sommer_holzapfel.h

#pragma once

#include <funcy/examples/volumetric_penalty_functions.h>
#include <funcy/funcy.h>

namespace funcy
{
    namespace detail
    {
        template < linalg::Matrix Mat, int n = linalg::dim< Mat >() >
        auto generate_incompressible_adipose_tissue_sommer_holzapfel( double cCells, double k1,
                                                                      double k2, double kappa,
                                                                      const Mat& A, const Mat& F )
        {
            using namespace linalg;
            auto S = LeftCauchyGreenStrainTensor( F );

            auto aniso = kappa * i1( F ) + ( 1 - 3 * kappa ) * i4( F, A ) - 1;
            auto f = cCells * ( i1( F ) - n ) + ( k1 / k2 ) * ( exp( k2 * squared( aniso ) ) - 1 );
            return f( S );
        }
    } // namespace detail

    template < linalg::Matrix Mat, int offset = linalg::dim< Mat >() >
    auto incompressible_adipose_tissue_sommer_holzapfel( double cCells, double k1, double k2,
                                                         double kappa, const Mat& A, const Mat& F )
    {
        return finalize(
            detail::generate_incompressible_adipose_tissue_sommer_holzapfel< Mat, offset >(
                cCells, k1, k2, kappa, A, F ) );
    }

    template < linalg::Matrix Mat, int offset = linalg::dim< Mat >() >
    auto incompressible_adipose_tissue_sommer_holzapfel( const Mat& A, const Mat& F )
    {
        return incompressible_adipose_tissue_sommer_holzapfel< Mat, offset >( 0.15, 0.8, 47.3, 0.09,
                                                                              A, F );
    }

    template < class Inflation, class Compression, linalg::Matrix Mat,
               int offset = linalg::dim< Mat >() >
    auto compressible_adipose_tissue_sommer_holzapfel( double cCells, double k1, double k2,
                                                       double kappa, double d0, double d1,
                                                       const Mat& M, const Mat& F )
    {
        return finalize(
            detail::generate_incompressible_adipose_tissue_sommer_holzapfel< Mat, offset >(
                cCells, k1, k2, kappa, M, F ) +
            volumetric_penalty< Inflation, Compression >( d0, d1, F ) );
    }

    template < class Inflation, class Compression, linalg::Matrix Mat,
               int offset = linalg::dim< Mat >() >
    auto compressible_adipose_tissue_sommer_holzapfel( double d0, double d1, const Mat& M,
                                                       const Mat& F )
    {
        return compressible_adipose_tissue_sommer_holzapfel< Inflation, Compression, Mat, offset >(
            0.15, 0.8, 47.3, 0.09, d0, d1, M, F );
    }
} // namespace funcy