template <int dim,typename NumberType>
void Solid<dim,NumberType>::output_results() const
{
DataOut<dim> data_out;
data_out.attach_dof_handler(dof_handler_ref);
std::vector<std::string> solution_names;
solution_names.emplace_back("u_x");
solution_names.emplace_back("u_y");
solution_names.emplace_back("u_z");
data_out.add_data_vector(solution_n, solution_names);
//// next part calculates Kirchhoff stress using scratch and pertaskdata
FE_DGQ<dim> history_fe (1);
DoFHandler<dim> history_dof_handler (triangulation);
history_dof_handler.distribute_dofs (history_fe);
std::vector< std::vector< Vector<double> > >
history_field (dim, std::vector< Vector<double> >(dim)),
local_history_values_at_qpoints (dim, std::vector< Vector<double> >(dim)),
local_history_fe_values (dim, std::vector< Vector<double> >(dim));
for (unsigned int i=0; i<dim; i++)
for (unsigned int j=0; j<dim; j++)
{
history_field[i][j].reinit(history_dof_handler.n_dofs());
local_history_values_at_qpoints[i][j].reinit(qf_cell.size());
local_history_fe_values[i][j].reinit(history_fe.dofs_per_cell);
}
FullMatrix<double> qpoint_to_dof_matrix (history_fe.dofs_per_cell,
qf_cell.size());
FETools::compute_projection_from_quadrature_points_matrix
(history_fe,
qf_cell, qf_cell,
qpoint_to_dof_matrix);
typename DoFHandler<dim>::active_cell_iterator cell = dof_handler_ref.
begin_active(),
endc = dof_handler_ref.end(),
dg_cell = history_dof_handler.begin_active();
for (; cell!=endc; ++cell, ++dg_cell)
{
/////////////////// added to define scratch and pertaskdata
//////////////////////////////////////
const UpdateFlags uf_cell(update_gradients |
update_JxW_values);
const UpdateFlags uf_face(update_values |
update_JxW_values);
typename Assembler_Base<dim,NumberType>::PerTaskData_ASM data(this);
typename Assembler_Base<dim,NumberType>::ScratchData_ASM scratch(fe, qf_cell,
uf_cell, qf_face, uf_face, solution_n);
const FEValuesExtractors::Vector &u_fe = data.solid->u_fe;
data.reset();
scratch.reset();
scratch.fe_values_ref.reinit(cell);
cell->get_dof_indices(data.local_dof_indices);
const std::vector<std::shared_ptr<const PointHistory<dim,NumberType> > >
lqph =
const_cast<const Solid<dim,NumberType> *>(data.solid)->
quadrature_point_history.get_data(cell);
Assert(lqph.size() == n_q_points, ExcInternalError());
scratch.fe_values_ref[u_fe].get_function_gradients(scratch.solution_total,
scratch.solution_grads_u_total);
///////////////////////////////////////////////////////////////////
for (unsigned int i=0; i<dim; i++)
for (unsigned int j=0; j<dim; j++)
{
for (unsigned int q=0; q<qf_cell.size(); ++q)
{
//////////////// added to compute Kirchhoff stress tau
//////////////////////////////////////
const Tensor<2,dim,NumberType> &grad_u = scratch.solution_grads_u_total[q];
const Tensor<2,dim,NumberType> F = Physics::Elasticity::Kinematics::F
(grad_u);
const NumberType det_F = determinant(F);
const Tensor<2,dim,NumberType> F_bar = Physics::Elasticity::Kinematics::
F_iso(F);
const SymmetricTensor<2,dim,NumberType> b_bar =
Physics::Elasticity::Kinematics::b(F_bar);
const SymmetricTensor<2,dim,NumberType> tau = lqph[q]->get_tau(det_F,b_bar);
//////////////////////////////////////////////////////////////
local_history_values_at_qpoints[i][j](q) = tau[i][j]; // tau_locally stored
}
qpoint_to_dof_matrix.vmult (local_history_fe_values[i][j],
local_history_values_at_qpoints[i][j]); // q_point to node extrapolation
dg_cell->set_dof_values (local_history_fe_values[i][j],
history_field[i][j]); // local to global dof transfer
}
}
data_out.add_data_vector(history_dof_handler,history_field[0][0],
"Kirchhoff_st_xx");
data_out.add_data_vector(history_dof_handler,history_field[1][1],
"Kirchhoff_st_yy");
data_out.add_data_vector(history_dof_handler,history_field[2][2],
"Kirchhoff_st_zz");
data_out.add_data_vector(history_dof_handler,history_field[0][1],
"Kirchhoff_st_xy");
data_out.add_data_vector(history_dof_handler,history_field[0][2],
"Kirchhoff_st_xz");
data_out.add_data_vector(history_dof_handler,history_field[1][2],
"Kirchhoff_st_yz");
// Since we are dealing with a large deformation problem, it would be nice
// to display the result on a displaced grid! The MappingQEulerian class
// linked with the DataOut class provides an interface through which this
// can be achieved without physically moving the grid points in the
// Triangulation object ourselves. We first need to copy the solution to
// a temporary vector and then create the Eulerian mapping. We also
// specify the polynomial degree to the DataOut object in order to produce
// a more refined output data set when higher order polynomials are used.
Vector<double> soln(solution_n.size());
for (unsigned int i = 0; i < soln.size(); ++i)
soln(i) = solution_n(i);
MappingQEulerian<dim> q_mapping(degree, dof_handler_ref, soln);
data_out.build_patches(q_mapping, degree);
std::ostringstream filename;
filename << "solution-" << time.get_timestep() << ".vtk";
std::ofstream output(filename.str().c_str());
data_out.write_vtk(output);
}
..............................................................................................
*error i am getting*
make all
Consolidate compiler generated dependencies of target cook_membrane
[ 50%] Building CXX object CMakeFiles/cook_membrane.dir/cook_membrane.cc.o
/home/newton/dealii-candi/deal.II-v9.4.2/examples/cook_membrane/cook_membrane.cc:
In instantiation of ‘void Cook_Membrane::Solid<dim,
NumberType>::output_results() const [with int dim = 3; NumberType =
Sacado::Fad::Exp::GeneralFad<Sacado::Fad::Exp::DynamicStorage<double,
double> >]’:
/home/newton/dealii-candi/deal.II-v9.4.2/examples/cook_membrane/cook_membrane.cc:1006:5:
required from ‘void Cook_Membrane::Solid<dim, NumberType>::run() [with int
dim = 3; NumberType =
Sacado::Fad::Exp::GeneralFad<Sacado::Fad::Exp::DynamicStorage<double,
double> >]’
/home/newton/dealii-candi/deal.II-v9.4.2/examples/cook_membrane/cook_membrane.cc:2421:23:
required from here
/home/newton/dealii-candi/deal.II-v9.4.2/examples/cook_membrane/cook_membrane.cc:2288:62:
error: cannot convert ‘std::vector<dealii::Tensor<2, 3,
Sacado::Fad::Exp::GeneralFad<Sacado::Fad::Exp::DynamicStorage<double,
double> > >, std::allocator<dealii::Tensor<2, 3,
Sacado::Fad::Exp::GeneralFad<Sacado::Fad::Exp::DynamicStorage<double,
double> > > > >’ to ‘std::vector<dealii::Tensor<2, 3, double>,
std::allocator<dealii::Tensor<2, 3, double> > >&’
2288 | scratch.solution_grads_u_total);
| ~~~~~~~~^~~~~~~~~~~~~~~~~~~~~~
| |
| std::vector<dealii::Tensor<2, 3,
Sacado::Fad::Exp::GeneralFad<Sacado::Fad::Exp::DynamicStorage<double,
double> > >, std::allocator<dealii::Tensor<2, 3,
Sacado::Fad::Exp::GeneralFad<Sacado::Fad::Exp::DynamicStorage<double,
double> > > > >
In file included from
/home/newton/dealii-candi/deal.II-v9.4.2/include/deal.II/grid/grid_tools.h:33,
from
/home/newton/dealii-candi/deal.II-v9.4.2/examples/cook_membrane/cook_membrane.cc:40:
/home/newton/dealii-candi/deal.II-v9.4.2/include/deal.II/fe/fe_values.h:1180:10:
note: initializing argument 2 of ‘void dealii::FEValuesViews::Vector<dim,
spacedim>::get_function_gradients(const InputVector&, std::vector<typename
dealii::ProductType<typename InputVector::value_type, dealii::Tensor<2,
spacedim> >::type>&) const [with InputVector = dealii::BlockVector<double>;
int dim = 3; int spacedim = 3; typename dealii::ProductType<typename
InputVector::value_type, dealii::Tensor<2, spacedim> >::type =
dealii::Tensor<2, 3, double>; typename InputVector::value_type = double]’
1179 | std::vector<solution_gradient_type<typename InputVector::value_type>>
| ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
1180 | &gradients) const;
| ~^~~~~~~~~
/home/newton/dealii-candi/deal.II-v9.4.2/examples/cook_membrane/cook_membrane.cc:2307:44:
error: cannot convert
‘dealii::internal::SymmetricTensorAccessors::Accessor<2, 3, true, 1,
Sacado::Fad::Exp::GeneralFad<Sacado::Fad::Exp::DynamicStorage<double,
double> > >::reference’ {aka
‘Sacado::Fad::Exp::GeneralFad<Sacado::Fad::Exp::DynamicStorage<double,
double> >’} to ‘double’ in assignment
2307 | local_history_values_at_qpoints[i][j](q) = tau[i][j]; // tau_locally
stored
| ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~^~~~~~~~~
/home/newton/dealii-candi/deal.II-v9.4.2/examples/cook_membrane/cook_membrane.cc:
In instantiation of ‘void Cook_Membrane::Solid<dim,
NumberType>::output_results() const [with int dim = 3; NumberType =
Sacado::Rad::ADvar<Sacado::Fad::Exp::GeneralFad<Sacado::Fad::Exp::DynamicStorage<double,
double> > >]’:
/home/newton/dealii-candi/deal.II-v9.4.2/examples/cook_membrane/cook_membrane.cc:1006:5:
required from ‘void Cook_Membrane::Solid<dim, NumberType>::run() [with int
dim = 3; NumberType =
Sacado::Rad::ADvar<Sacado::Fad::Exp::GeneralFad<Sacado::Fad::Exp::DynamicStorage<double,
double> > >]’
/home/newton/dealii-candi/deal.II-v9.4.2/examples/cook_membrane/cook_membrane.cc:2434:23:
required from here
/home/newton/dealii-candi/deal.II-v9.4.2/examples/cook_membrane/cook_membrane.cc:2288:62:
error: cannot convert ‘std::vector<dealii::Tensor<2, 3,
Sacado::Rad::ADvar<Sacado::Fad::Exp::GeneralFad<Sacado::Fad::Exp::DynamicStorage<double,
double> > > >, std::allocator<dealii::Tensor<2, 3,
Sacado::Rad::ADvar<Sacado::Fad::Exp::GeneralFad<Sacado::Fad::Exp::DynamicStorage<double,
double> > > > > >’ to ‘std::vector<dealii::Tensor<2, 3, double>,
std::allocator<dealii::Tensor<2, 3, double> > >&’
2288 | scratch.solution_grads_u_total);
| ~~~~~~~~^~~~~~~~~~~~~~~~~~~~~~
| |
| std::vector<dealii::Tensor<2, 3,
Sacado::Rad::ADvar<Sacado::Fad::Exp::GeneralFad<Sacado::Fad::Exp::DynamicStorage<double,
double> > > >, std::allocator<dealii::Tensor<2, 3,
Sacado::Rad::ADvar<Sacado::Fad::Exp::GeneralFad<Sacado::Fad::Exp::DynamicStorage<double,
double> > > > > >
In file included from
/home/newton/dealii-candi/deal.II-v9.4.2/include/deal.II/grid/grid_tools.h:33,
from
/home/newton/dealii-candi/deal.II-v9.4.2/examples/cook_membrane/cook_membrane.cc:40:
/home/newton/dealii-candi/deal.II-v9.4.2/include/deal.II/fe/fe_values.h:1180:10:
note: initializing argument 2 of ‘void dealii::FEValuesViews::Vector<dim,
spacedim>::get_function_gradients(const InputVector&, std::vector<typename
dealii::ProductType<typename InputVector::value_type, dealii::Tensor<2,
spacedim> >::type>&) const [with InputVector = dealii::BlockVector<double>;
int dim = 3; int spacedim = 3; typename dealii::ProductType<typename
InputVector::value_type, dealii::Tensor<2, spacedim> >::type =
dealii::Tensor<2, 3, double>; typename InputVector::value_type = double]’
1179 | std::vector<solution_gradient_type<typename InputVector::value_type>>
| ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
1180 | &gradients) const;
| ~^~~~~~~~~
/home/newton/dealii-candi/deal.II-v9.4.2/examples/cook_membrane/cook_membrane.cc:2307:44:
error: cannot convert
‘dealii::internal::SymmetricTensorAccessors::Accessor<2, 3, true, 1,
Sacado::Rad::ADvar<Sacado::Fad::Exp::GeneralFad<Sacado::Fad::Exp::DynamicStorage<double,
double> > > >::reference’ {aka
‘Sacado::Rad::ADvar<Sacado::Fad::Exp::GeneralFad<Sacado::Fad::Exp::DynamicStorage<double,
double> > >’} to ‘double’ in assignment
2307 | local_history_values_at_qpoints[i][j](q) = tau[i][j]; // tau_locally
stored
| ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~^~~~~~~~~
make[2]: *** [CMakeFiles/cook_membrane.dir/build.make:76:
CMakeFiles/cook_membrane.dir/cook_membrane.cc.o] Error 1
make[1]: *** [CMakeFiles/Makefile2:90: CMakeFiles/cook_membrane.dir/all]
Error 2
make: *** [Makefile:91: all] Error 2
"make all" terminated with exit code 2. Build might be incomplete.
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/4ee059e6-9311-4f37-b42d-5ff41aac43d4n%40googlegroups.com.
