Sacado-auxiliary_functions.h File Reference

#include <deal.II/base/symmetric_tensor.h>
#include <iostream>
#include <Sacado.hpp>

Include dependency graph for Sacado-auxiliary_functions.h:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Typedefs
using	fad_double = Sacado::Fad::DFad< double >

Functions
template<int dim>
SymmetricTensor< 2, dim, fad_double >	init_Sacado (const SymmetricTensor< 2, dim, double > &SymTensor)
template<int dim>
SymmetricTensor< 2, dim >	extract_value_from_Sacado (const SymmetricTensor< 2, dim, fad_double > &SymTensor_Sacado)
template<int dim>
Vector< fad_double >	reinit (unsigned int a, bool b)

Typedef Documentation

fad_double

using fad_double = Sacado::Fad::DFad<double>

Function Documentation

extract_value_from_Sacado()

template<int dim>

SymmetricTensor<2,dim> extract_value_from_Sacado

(

const SymmetricTensor< 2, dim, fad_double > &

SymTensor_Sacado

)

{
    SymmetricTensor<2,dim> SymTensor;
    // .val() extracts the \a double value from the Sacado variable
     for ( unsigned int i=0; i<dim; ++i)
        for ( unsigned int j=i; j<dim; ++j )
            SymTensor[i][j] = SymTensor_Sacado[i][j].val();
     
     return SymTensor;
}

init_Sacado()

template<int dim>

SymmetricTensor<2,dim,fad_double> init_Sacado

(

const SymmetricTensor< 2, dim, double > &

SymTensor

)

{
    SymmetricTensor<2,dim,fad_double> SymTensor_Sacado;
     for ( unsigned int i=0; i<dim; ++i)
        for ( unsigned int j=i; j<dim; ++j )
            SymTensor_Sacado[i][j] = SymTensor[i][j];
     
     return SymTensor_Sacado;
}

reinit()

template<int dim>

Vector<fad_double> reinit	(	unsigned int	a,
		bool	b
	)

{
    Vector<fad_double> abc (5);
     return abc;
}