I tried to modify the right side of example 52 from one equation in
parallel to three equations in parallel, each with own system_matrix,
mass_matrix and source term.
My code for that is
template<int dim>
Table<1, Vector<double>> Diffusion<dim>::evaluate_diffusion_full (const
double time, const Table<1, Vector<double>> &y) const
{
Vector<double> tmp_N(dof_handler.n_dofs());
Vector<double> tmp_TE(dof_handler.n_dofs());
Vector<double> tmp_TL(dof_handler.n_dofs());
tmp_N = 0.;
tmp_TE = 0.;
tmp_TL = 0.;
Vector<double> y_N = y[0];
Vector<double> y_TE = y[1];
Vector<double> y_TL = y[2];
system_matrix_N.vmult(tmp_N,y_N);
system_matrix_TE.vmult(tmp_TE,y_TE);
system_matrix_TL.vmult(tmp_TL,y_TL);
const QGauss<dim> quadrature_formula(fe_degree+1);
FEValues<dim> fe_values(fe,
quadrature_formula,
update_values | update_quadrature_points |
update_JxW_values);
const unsigned int dofs_per_cell = fe.dofs_per_cell;
const unsigned int n_q_points = quadrature_formula.size();
Vector<double> cell_source_N(dofs_per_cell), cell_source_TE(
dofs_per_cell), cell_source_TL(dofs_per_cell);
std::vector<types::global_dof_index> local_dof_indices_N (
dofs_per_cell), local_dof_indices_TE(dofs_per_cell), local_dof_indices_TL (
dofs_per_cell);
for (auto cell = dof_handler.begin_active(); cell!=dof_handler.end
(); ++cell)
{
cell_source_N = 0.;
cell_source_TE = 0.;
cell_source_TL = 0.;
fe_values.reinit (cell);
for (unsigned int q_point=0; q_point<n_q_points; ++q_point)
{
const double source_N = get_source_N(time, fe_values.
quadrature_point(q_point), 0, 0, 0);
const double source_TE = get_source_TE(time, fe_values.
quadrature_point(q_point), 0, 0);
const double source_TL = get_source_TL();
for (unsigned int i=0; i<dofs_per_cell; ++i)
{
cell_source_N(i) += source_N * fe_values.shape_value(i,
q_point) * fe_values.JxW(q_point);//Not sure if that is correct?
cell_source_TE(i) += source_TE * fe_values.shape_value(i
, q_point) * fe_values.JxW(q_point);
cell_source_TL(i) += source_TL * fe_values.shape_value(i
, q_point) * fe_values.JxW(q_point);
}
}
cell->get_dof_indices(local_dof_indices_N);
constraint_matrix.distribute_local_to_global(cell_source_N,
local_dof_indices_N, tmp_N);
cell->get_dof_indices(local_dof_indices_TE);
constraint_matrix.distribute_local_to_global(cell_source_TE,
local_dof_indices_TE,tmp_TE);
cell->get_dof_indices(local_dof_indices_TL);
constraint_matrix.distribute_local_to_global(cell_source_TL,
local_dof_indices_TL, tmp_TL);
}
Vector<double> value_N(dof_handler.n_dofs());
Vector<double> value_TE(dof_handler.n_dofs());
Vector<double> value_TL(dof_handler.n_dofs());
inverse_mass_matrix_N.vmult(value_N,tmp_N);
inverse_mass_matrix_TE.vmult(value_TE,tmp_TE);
inverse_mass_matrix_TL.vmult(value_TL,tmp_TL);
const Vector<double> vector_values[] = {value_N, value_TE, value_TL
};//Have to put all three vectors into one Table in order to return them
return Table<1, Vector<double>>(3, vector_values, true);
}
Now, to integrate that using the DOPRI-solver, I changed the solving method
accordingly:
Table<1, Vector<double>> full_solution = 0.;
TimeStepping::EmbeddedExplicitRungeKutta<Table<1, Vector<double>> >
embedded_explicit_runge_kutta(TimeStepping::DOPRI,
coarsen_param,
refine_param,
min_delta,
max_delta,
refine_tol,
coarsen_tol);
output_results (0, TimeStepping::DOPRI);
// Now for the time loop. The last time step is chosen such that
the final
// time is exactly reached.
unsigned int n_steps=0;
while (time < final_time)
{
if (time+time_step>final_time)
time_step = final_time-time;
time = embedded_explicit_runge_kutta.evolve_one_time_step(
std::bind(&Diffusion::evaluate_diffusion_full,this,
std::placeholders::_1,std::placeholders::_2),
time,time_step, full_solution);
if ((n_steps+1)%10==0)
output_results(n_steps+1,TimeStepping::DOPRI);
time_step = embedded_explicit_runge_kutta.get_status().
delta_t_guess;
++n_steps;
}
When compiling the new code, I get a linker error:
CMakeFiles/nuclear_modified.dir/main.cpp.o: In function
`dealii::TimeStepping::EmbeddedExplicitRungeKutta<dealii::Table<1,
dealii::Vector<double> > >::~EmbeddedExplicitRungeKutta()':
main.cpp:(.text._ZN6dealii12TimeStepping26EmbeddedExplicitRungeKuttaINS_5TableILi1ENS_6VectorIdEEEEED2Ev[_ZN6dealii12TimeStepping26EmbeddedExplicitRungeKuttaINS_5TableILi1ENS_6VectorIdEEEEED5Ev]+0x1c):
undefined reference to `dealii::TimeStepping::EmbeddedExplicitRungeKutta<
dealii::Table<1, dealii::Vector<double> > >::free_memory()'
CMakeFiles/nuclear_modified.dir/main.cpp.o: In function
`dealii::TimeStepping::EmbeddedExplicitRungeKutta<dealii::Table<1,
dealii::Vector<double> > >::~EmbeddedExplicitRungeKutta()':
main.cpp:(.text.
_ZN6dealii12TimeStepping26EmbeddedExplicitRungeKuttaINS_5TableILi1ENS_6VectorIdEEEEED0Ev
[
_ZN6dealii12TimeStepping26EmbeddedExplicitRungeKuttaINS_5TableILi1ENS_6VectorIdEEEEED5Ev
]+0x1c): undefined reference to
`dealii::TimeStepping::EmbeddedExplicitRungeKutta<dealii::Table<1,
dealii::Vector<double> > >::free_memory()'
CMakeFiles/nuclear_modified.dir/main.cpp.o: In function `Step52::Diffusion<3
>::embedded_explicit_method(unsigned int, double, double)':
main.cpp:(.text._ZN6Step529DiffusionILi3EE24embedded_explicit_methodEjdd[_ZN6Step529DiffusionILi3EE24embedded_explicit_methodEjdd]+0x1aa):
undefined reference to
`dealii::TimeStepping::EmbeddedExplicitRungeKutta<dealii::Table<1,
dealii::Vector<double> >
>::EmbeddedExplicitRungeKutta(dealii::TimeStepping::runge_kutta_method,
double, double, double, double, double, double)'
main.cpp:(.text._ZN6Step529DiffusionILi3EE24embedded_explicit_methodEjdd[
_ZN6Step529DiffusionILi3EE24embedded_explicit_methodEjdd]+0x261): undefined
reference to `dealii::TimeStepping::EmbeddedExplicitRungeKutta<dealii::Table<1,
dealii::Vector<double> >
>::evolve_one_time_step(std::function<dealii::Table<1,
dealii::Vector<double> > (double, dealii::Table<1, dealii::Vector<double> >
const&)>, double, double, dealii::Table<1, dealii::Vector<double> >&)'
main.cpp:(.text._ZN6Step529DiffusionILi3EE24embedded_explicit_methodEjdd[_ZN6Step529DiffusionILi3EE24embedded_explicit_methodEjdd]+0x2a8):
undefined reference to `dealii::TimeStepping::EmbeddedExplicitRungeKutta<
dealii::Table<1, dealii::Vector<double> > >::get_status() const'
main.cpp:(.text._ZN6Step529DiffusionILi3EE24embedded_explicit_methodEjdd[_ZN6Step529DiffusionILi3EE24embedded_explicit_methodEjdd]+0x2df):
undefined reference to
`dealii::TimeStepping::EmbeddedExplicitRungeKutta<dealii::Table<1,
dealii::Vector<double> > >::free_memory()'
CMakeFiles/nuclear_modified.dir/main.cpp.o:(.data.rel.ro.
_ZTVN6dealii12TimeStepping10RungeKuttaINS_5TableILi1ENS_6VectorIdEEEEEE[
_ZTVN6dealii12TimeStepping10RungeKuttaINS_5TableILi1ENS_6VectorIdEEEEEE]+
0x20): undefined reference to
`dealii::TimeStepping::RungeKutta<dealii::Table<1,
dealii::Vector<double> >
>::evolve_one_time_step(std::vector<std::function<dealii::Table<1,
dealii::Vector<double> > (double, dealii::Table<1, dealii::Vector<double> >
const&)>, std::allocator<std::function<dealii::Table<1,
dealii::Vector<double> > (double, dealii::Table<1, dealii::Vector<double> >
const&)> > >&, std::vector<std::function<dealii::Table<1,
dealii::Vector<double> > (double, double, dealii::Table<1,
dealii::Vector<double> > const&)>,
std::allocator<std::function<dealii::Table<1, dealii::Vector<double> >
(double, double, dealii::Table<1, dealii::Vector<double> > const&)> > >&,
double, double, dealii::Table<1, dealii::Vector<double> >&)'
CMakeFiles/nuclear_modified.dir/main.cpp.o:(.data.rel.ro._ZTVN6dealii12TimeStepping26EmbeddedExplicitRungeKuttaINS_5TableILi1ENS_6VectorIdEEEEEE[_ZTVN6dealii12TimeStepping26EmbeddedExplicitRungeKuttaINS_5TableILi1ENS_6VectorIdEEEEEE]+0x20):
undefined reference to `dealii::TimeStepping::RungeKutta<dealii::Table<1,
dealii::Vector<double> > >::evolve_one_time_step(std::vector<std::function<
dealii::Table<1, dealii::Vector<double> > (double, dealii::Table<1, dealii::
Vector<double> > const&)>, std::allocator<std::function<dealii::Table<1,
dealii::Vector<double> > (double, dealii::Table<1, dealii::Vector<double> >
const&)> > >&, std::vector<std::function<dealii::Table<1, dealii::Vector
<double> > (double, double, dealii::Table<1, dealii::Vector<double> > const
&)>, std::allocator<std::function<dealii::Table<1, dealii::Vector<double> >
(double, double, dealii::Table<1, dealii::Vector<double> > const&)> > >&,
double, double, dealii::Table<1, dealii::Vector<double> >&)'
CMakeFiles/nuclear_modified.dir/main.cpp.o:(.data.rel.ro._ZTVN6dealii12TimeStepping26EmbeddedExplicitRungeKuttaINS_5TableILi1ENS_6VectorIdEEEEEE[_ZTVN6dealii12TimeStepping26EmbeddedExplicitRungeKuttaINS_5TableILi1ENS_6VectorIdEEEEEE]+0x28):
undefined reference to
`dealii::TimeStepping::EmbeddedExplicitRungeKutta<dealii::Table<1,
dealii::Vector<double> > >::get_status() const'
CMakeFiles/nuclear_modified.dir/main.cpp.o:(.data.rel.ro.
_ZTVN6dealii12TimeStepping26EmbeddedExplicitRungeKuttaINS_5TableILi1ENS_6VectorIdEEEEEE
[
_ZTVN6dealii12TimeStepping26EmbeddedExplicitRungeKuttaINS_5TableILi1ENS_6VectorIdEEEEEE
]+0x30): undefined reference to
`dealii::TimeStepping::EmbeddedExplicitRungeKutta<dealii::Table<1,
dealii::Vector<double> >
>::initialize(dealii::TimeStepping::runge_kutta_method)'
CMakeFiles/nuclear_modified.dir/main.cpp.o:(.data.rel.ro._ZTVN6dealii12TimeStepping26EmbeddedExplicitRungeKuttaINS_5TableILi1ENS_6VectorIdEEEEEE[_ZTVN6dealii12TimeStepping26EmbeddedExplicitRungeKuttaINS_5TableILi1ENS_6VectorIdEEEEEE]+0x38):
undefined reference to `dealii::TimeStepping::EmbeddedExplicitRungeKutta<
dealii::Table<1, dealii::Vector<double> > >::evolve_one_time_step(std::
function<dealii::Table<1, dealii::Vector<double> > (double, dealii::Table<1,
dealii::Vector<double> > const&)>, std::function<dealii::Table<1, dealii::
Vector<double> > (double, double, dealii::Table<1, dealii::Vector<double> >
const&)>, double, double, dealii::Table<1, dealii::Vector<double> >&)'
collect2: error: ld returned 1 exit status
Does that mean that my approach to solve three equations at once is not
valid? If yes, how can I fix that problem?
--
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.
For more options, visit https://groups.google.com/d/optout.
