Giuseppe, Le mar. 22 août 2023 à 13:27, giuseppe orlando <giusor...@gmail.com> a écrit :
> > 1. For 'additional_data.mapping_update_flags', I need to specify > 'update_gradients' even though the resulting matrix for the linear system > is nothing but a mass matrix. If I do not specify 'update_gradients', I get > a segmentation fault during 'reinit' of the MatrixFree structure. Did > anyone find out the same behaviour? > > This should be fixed in master. If it's not, let us know. > > 1. I would need to consider periodic boundary conditions for a > configuration of interest, however this results in a blow-up of the > simulation. Does it mean that the matrix-free framework associated with > CUDA can handle only Dirichlet boundary conditions? I was already able with > the previous version of deal.II to write down a matrix-based version for > GPU with periodic bcs and matrix-based and matrix-free codes for GPU with > Dirichlet bcs provide the same results and, therefore, I was wondering if > there is an underlying issue in handling periodic bcs in a matrix-free > framework for CUDA. > > Periodic boundary conditions are not supported currently. I don't know much about periodic boundaries, so I don't know how much work is required to support that. Maybe someone else knows? Best, Bruno -- 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/CAGVt9eMz7G6yR92fibDjTRra1EgOWoixfzyS%3D8xaWdAB0hxvvw%40mail.gmail.com.
