Yes, like
here
https://github.com/dealii/dealii/blob/b84270a1d4099292be5b3d43c2ea65f3ee005919/tests/matrix_free/pre_and_post_loops_01.cc#L100-L121
On Saturday, 18 January 2020 12:57:24 UTC+1, Maxi Miller wrote:
>
> In step-48 the inverse mass matrix is applied by moving the inverse data
> into a vector and applying the function scale(), i.e. as in the following
> code:
> data.cell_loop(&LaplaceOperator::local_apply_cell,
> this,
> dst,
> src,
> true);
> computing_times[0] += timer.wall_time();
> timer.restart();
>
> dst.scale(inv_mass_matrix);
>
>
>
> Now I would like to do the same, but use a cell_loop instead of scale().
> Thus, I created an additional function called
> "local_apply_inverse_mass_matrix" as
> template <int dim, int degree, int n_points_1d>
> void LaplaceOperator<dim, degree, n_points_1d>::
> local_apply_inverse_mass_matrix(
> const MatrixFree<dim, Number> & data,
> LinearAlgebra::distributed::Vector<Number> & dst,
> const LinearAlgebra::distributed::Vector<Number> &src,
> const std::pair<unsigned int, unsigned int> & cell_range)
> const
> {
> (void) data;
> (void) cell_range;
> dst = src;
> dst.scale(inv_mass_matrix);
> }
>
> When running that code, though, using
> LinearAlgebra::distributed::Vector<Number> tmp(src);
>
> data.initialize_dof_vector(tmp);
> data.initialize_dof_vector(dst);
> data.cell_loop(&LaplaceOperator::local_apply_cell,
> this,
> tmp,
> src,
> true);
> computing_times[0] += timer.wall_time();
> timer.restart();
>
> data.cell_loop(&LaplaceOperator::local_apply_inverse_mass_matrix,
> this,
> dst,
> tmp,
> true);
> computing_times[1] += timer.wall_time();
>
> computing_times[3] += 1.;
>
>
> I get the error
> An error occurred in line <3338> of file </opt/dealii/include/deal.II/
> matrix_free/matrix_free.h> in function
> void dealii::internal::VectorDataExchange<dim, Number,
> VectorizedArrayType>::compress_start(unsigned int, VectorType&) [with
> VectorType = dealii::LinearAlgebra::distributed::Vector<double>; typename
> std::enable_if<(dealii::internal::has_compress_start<VectorType>::value &&
> dealii::internal::has_exchange_on_subset<VectorType>::value), VectorType
> >::type* <anonymous> = 0; int dim = 2; Number = double;
> VectorizedArrayType = dealii::VectorizedArray<double, 4>]
> The violated condition was:
> vec.has_ghost_elements() == false
> Additional information:
> You are trying to use functionality in deal.II that is currently not
> implemented. In many cases, this indicates that there simply didn't
> appear much of a need for it, or that the author of the original code did
> not have the time to implement a particular case. If you hit this
> exception, it is therefore worth the time to look into the code to find out
> whether you may be able to implement the missing functionality. If you do,
> please consider providing a patch to the deal.II development sources (see
> the deal.II website on how to contribute).
>
> Stacktrace:
> -----------
> #0 ./MF_FES_RK4-Test: void dealii::internal::VectorDataExchange<2,
> double, dealii::VectorizedArray<double, 4>
> >::compress_start<dealii::LinearAlgebra::distributed::Vector<double,
> dealii::MemorySpace::Host>,
> (dealii::LinearAlgebra::distributed::Vector<double,
> dealii::MemorySpace::Host>*)0>(unsigned int,
> dealii::LinearAlgebra::distributed::Vector<double,
> dealii::MemorySpace::Host>&)
> #1 ./MF_FES_RK4-Test: void dealii::internal::compress_start<2,
> dealii::LinearAlgebra::distributed::Vector<double,
> dealii::MemorySpace::Host>, double, dealii::VectorizedArray<double, 4>,
> (dealii::LinearAlgebra::distributed::Vector<double,
> dealii::MemorySpace::Host>*)0>(dealii::LinearAlgebra::distributed::Vector<double,
>
> dealii::MemorySpace::Host>&, dealii::internal::VectorDataExchange<2,
> double, dealii::VectorizedArray<double, 4> >&, unsigned int)
> #2 ./MF_FES_RK4-Test: dealii::internal::MFWorker<dealii::MatrixFree<2,
> double, dealii::VectorizedArray<double, 4> >,
> dealii::LinearAlgebra::distributed::Vector<double,
> dealii::MemorySpace::Host>,
> dealii::LinearAlgebra::distributed::Vector<double,
> dealii::MemorySpace::Host>, Step40::LaplaceOperator<2, 2, 4>,
> true>::vector_compress_start()
> #3 /opt/dealii/lib/libdeal_II.g.so.9.2.0-pre:
> dealii::internal::MatrixFreeFunctions::MPICommunication::execute()
> #4
>
> /opt/intel/compilers_and_libraries_2019.5.281/linux/tbb/lib/intel64_lin/gcc4.7/libtbb_debug.so.2:
>
>
> #5
>
> /opt/intel/compilers_and_libraries_2019.5.281/linux/tbb/lib/intel64_lin/gcc4.7/libtbb_debug.so.2:
>
>
> #6
>
> /opt/intel/compilers_and_libraries_2019.5.281/linux/tbb/lib/intel64_lin/gcc4.7/libtbb_debug.so.2:
>
>
> #7
>
> /opt/intel/compilers_and_libraries_2019.5.281/linux/tbb/lib/intel64_lin/gcc4.7/libtbb_debug.so.2:
>
>
> #8
>
> /opt/intel/compilers_and_libraries_2019.5.281/linux/tbb/lib/intel64_lin/gcc4.7/libtbb_debug.so.2:
>
>
> #9 /lib64/libpthread.so.0:
> #10 /lib64/libc.so.6: clone
>
> Why that, and how can I fix it?
> Thanks!
>
--
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/44155e0b-4ca1-41b2-9875-97db0951b48d%40googlegroups.com.
