Hi all,
I had a legacy code were PETScWrappers::MPI::Vectors were used to for
parallel computing. For very large problems, I was told to move to
Trilinos. It was fairly easy to do so for most part of the code but I am
stuck in the implementation of boundary values. Currently the boundary
values are implemented as follows:
template<int dim>
void ElasticityEquation<dim>::apply_boundary_conditions(const double time,
const double velocity)
{
// Begin timer
TimerOutput::Scope t(computing_timer, "Apply Boundary Condition");
boundary_values.clear();
VectorTools::interpolate_boundary_values
(dof_handler,2,ZeroFunction<dim>(dim),boundary_values);
VectorTools::interpolate_boundary_values
(dof_handler,3,BoundaryValues<dim>(time, velocity) , boundary_values);
// The range of owned dofs
std::pair<unsigned int, unsigned int> range = solution.local_range();
// Set the boundary conditions
for (typename std::map<unsigned int, double>::const_iterator i =
boundary_values.begin(); i != boundary_values.end(); ++i){
if (i->first >= range.first && i->first < range.second){
completely_distributed_solution(i->first) = i->second;
}
}
completely_distributed_solution.compress(VectorOperation::insert);
hanging_constraints.distribute( completely_distributed_solution);
solution = completely_distributed_solution;
}
This used to work fine with PETSc but not with Trilinos. The issue seems to
be the fact the Trilinos does not arrange the vector contiguous ranges of
elements (
https://www.dealii.org/current/doxygen/deal.II/classTrilinosWrappers_1_1MPI_1_1Vector.html#abcc22742227d63950cb0995223e6ee32)
. Is there a standard workaround?
Thanks
Vivek
--
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 [email protected].
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/73ea3ef9-b5f8-461d-81c0-e43c0b4d8fbb%40googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
