Periodic bc for dim=1

39 views
Skip to first unread message

Praveen C

unread,
Feb 8, 2022, 12:34:42 AM2/8/22
to Deal. II Googlegroup
Dear all

I was trying to write a code that works for 1/2/3d navier stokes.

But it seems like collect_periodic_faces does not work for 1d ? Is this the case, and is there any other way to deal with this situation ?

Thanks
praveen

Peter Munch

unread,
Feb 8, 2022, 2:49:58 AM2/8/22
to deal.II User Group
What do you mean by "does not work for 1D"? In hyper.deal (https://github.com/hyperdeal/hyperdeal), we also have PBC in "1D".

Peter

Praveen C

unread,
Feb 8, 2022, 3:01:02 AM2/8/22
to Deal. II Googlegroup
I am getting compilation error

[100%] Linking CXX executable main
Undefined symbols for architecture arm64:
"void dealii::GridTools::collect_periodic_faces<dealii::parallel::distributed::Triangulation<1, 1> >(dealii::parallel::distributed::Triangulation<1, 1> const&, unsigned int, unsigned int, int, std::__1::vector<dealii::GridTools::PeriodicFacePair<dealii::parallel::distributed::Triangulation<1, 1>::cell_iterator>, std::__1::allocator<dealii::GridTools::PeriodicFacePair<dealii::parallel::distributed::Triangulation<1, 1>::cell_iterator> > >&, dealii::Tensor<1, dealii::parallel::distributed::Triangulation<1, 1>::space_dimension, double> const&, dealii::FullMatrix<double> const&)", referenced from:
CGSEM<1>::setup_box_grid() in cgsem.cc.o
ld: symbol(s) not found for architecture arm64
clang: error: linker command failed with exit code 1 (use -v to see invocation)

Maybe I am doing something wrong.

Thanks
praveen
> --
> 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+un...@googlegroups.com.
> To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/3dec29d1-bb19-4f58-8c8f-35568ceaf754n%40googlegroups.com.

Marc Fehling

unread,
Feb 8, 2022, 9:40:06 AM2/8/22
to deal.II User Group
Hello Praveen,

the parallel::distributed::Triangulation class is not instantiated for dim=1, only for dim>1. The reason behind it is that p4est does not allow to create 1D domains. With any other Triangulation class you should be fine.

Marc
Reply all
Reply to author
Forward
0 new messages