Sorry I forgot to present the error.
--------------------------------------------------------
An error occurred in line <164> of file
</home/ubuntu/deal.II/dealii/source/base/partitioner.cc> in function
void dealii::Utilities::MPI::Partitioner::set_owned_indices(const
dealii::IndexSet&)
The violated condition was:
locally_owned_indices.is_contiguous() == true
Additional information:
The index set specified in locally_owned_indices is not contiguous.
Stacktrace:
-----------
#0 /home/ubuntu/deal.II/installed/lib/libdeal_II.g.so.9.4.0:
dealii::Utilities::MPI::Partitioner::set_owned_indices(dealii::IndexSet
const&)
#1 /home/ubuntu/deal.II/installed/lib/libdeal_II.g.so.9.4.0:
dealii::Utilities::MPI::Partitioner::Partitioner(dealii::IndexSet const&,
ompi_communicator_t* const&)
...
#9 /home/ubuntu/deal.II/installed/lib/libdeal_II.g.so.9.4.0:
std::shared_ptr<dealii::Utilities::MPI::Partitioner>
std::make_shared<dealii::Utilities::MPI::Partitioner, dealii::IndexSet
const&, ompi_communicator_t* const&>(dealii::IndexSet const&,
ompi_communicator_t* const&)
#10 /home/ubuntu/deal.II/installed/lib/libdeal_II.g.so.9.4.0:
dealii::LinearAlgebra::distributed::Vector<double,
dealii::MemorySpace::Host>::reinit(dealii::IndexSet const&,
ompi_communicator_t* const&)
#11 /home/ubuntu/deal.II/installed/lib/libdeal_II.g.so.9.4.0:
dealii::LinearAlgebra::distributed::Vector<double, dealii::MemorySpace::
在2023年2月24日星期五 UTC+8 19:37:17<yy.wayne> 写道：
> Step-32 and step-50 are important for parallelization of TrilinosWrappered
> data type.
> For someone has to write codes with Trilinos data type, here are necessary
> modification
> from serial to parallel programs:
>
> 1. Change triangulation to parallel::distributed::triangulation
> rather than
> GridTools::partition_triangulation()
>
> <https://www.dealii.org/current/doxygen/deal.II/namespaceGridTools.html#a99eba8e3b388258eda37a2724579dd1d>
>
> or parallel::shared::Triangulation
>
> <https://www.dealii.org/current/doxygen/deal.II/classparallel_1_1shared_1_1Triangulation.html>
> .
> 2. SparsityPattern and solution/rhs Vector initialization specified
> for Trilinos, see Step-42
> 3. Choose locally owned cells or level cells during assembling.
>
> Though they both allow parallel, step-42's DoFTools::make_sparsity_pattern
> <https://www.dealii.org/current/doxygen/deal.II/group__constraints.html#gaf78e864edbfba7e0a7477457bfb96b26>
> has
> better performance
> in matrix LU decomposition (30%-50% less time) than tep-50's.
>
> However another problem pops out. When test a vector-valued problem,
> renumber DoFs
> component wise make reinit of vectors fail. If
> DoFRenumbering::component_wise is called,
>
> 1. LinearAlgebra::distributed::Vector<double>
> v1(dof_handler.locally_owned_dofs(),
> mpi_communicator); v1.reinit(dof_handler.locally_owned_dofs(),
> mpi_communicator);
> Success
> 2. LinearAlgebra::distributed::Vector<double> v2;
> v2.reinit(dof_handler.locally_owned_dofs(), mpi_communicator);
> Fails.
>
> Is that relates to trilinos date type?
> 在2023年2月23日星期四 UTC+8 22:44:57<Wolfgang Bangerth> 写道：
>
>> On 2/23/23 07:42, 'yy.wayne' via deal.II User Group wrote:
>> > So if step-4 doesn't support parallel then what I got is reasonable.
>>
>> Correct.
>> W.
>>
>> --
>> ------------------------------------------------------------------------
>> Wolfgang Bangerth email: bang...@colostate.edu
>> www: http://www.math.colostate.edu/~bangerth/
>>
>>
>>
--
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/4e07362c-e779-4973-ae5e-44b864465409n%40googlegroups.com.
