[deal.II] FunctionManifold class

Good morning everyone, I would have a question about the creation of a triangulation. Is it possible to use the 'FunctionManifold' class to create a boundary described by an analytical function? If so, how can I attach this manifold and "substitute" for instance the default straight boundary of

Re: [deal.II] FunctionManifold class

e of examples on how to use FunctionManifold with an analytical > function. Most of them use std::strings but > tests/fe/fe_values_function_manifold.cc and > tests/mappings/mapping_q_manifold_02.cc actually use Function objects. > > Best, > Daniel > > Am Mo., 31. Jan. 2022 um 02

[deal.II] Integral over different lines

Good morning everyone, I would have a question about the computation of integrals. Suppose I have a rectangular domain, is it there any way to compute the integral of a certain function over different heights (not only on the bottom or top boundary)? Thanks in advance, Giuseppe -- The deal.II

[deal.II] Minisymposium at Enumath 2023

would be interested in organizing a mini-symposium on “Computational Fluid Dynamics Applications based on *deal.II*”, possibly in collaboration with some of the developers. Please, let me know if this can be of interest for you. Best, Giuseppe Orlando -- The deal.II project is located at

[deal.II] CUDA facilities

Hello everyone, I am writing this post to ask a couple of questions about the CUDA facilities in deal.II, after a look on step-64 tutorial. 1. Searching for the documentation, I noticed that routines such as CUDAWrappers::FEFaceEvalution or CUDAWrappers::MatrixFree::loop are not defined. Am I

R: [deal.II] Re: CUDA facilities

Dear Bruno,thanks for the kind response. I will certainly have a look to the files that you suggested me about the CUDAWrappers::SparseMatrix (and also the other tests included therein). Just as curiosity (if I am not inappropriate and you can answer me): are these new GPU facilities with Kokkos al

[deal.II] MatrixFree spacedim different than dim

Hello everyone, I am writing this post to ask a curiosity about the MatrixFree infrastructure. Is it possible to employ it for solving PDEs on surface (like step-38 tutorial). Indeed, looking at the documentation of the reinit function, it seems that only a 'DofHanlder' instance is accepted, wh

[deal.II] DG Euler eqautions (step-33 and step-67)

Dear Deal.II community, I'm currently working on DG formulation for Euler equations with implicit time-stepping, trying to combine in some sense the two approaches of ste-33 and step-67. I'm wondering if there is a way to employ automatic differentiation in the matrix-free approach; I think tha

[deal.II] Issue on CUDA code

Hello everyone, I'm writing this post to ask some help for a CUDA code. I have implemented a matrix-based CUDA coda for the so-called solid body rotation, namely a time dependent advection problem. In my first implementation, I encapulased the computation of the matrix and of the right-hand side

[deal.II] Re: Issue on CUDA code

Here, you can find the two codes. 'SBR_GPU' is the one that arrives up to the end, whereas 'SBR_GPU_WRONG' is the one that stops with the error described above once I call 'output_results'. If the code is too long and can be helpful, I could try to reduce it avoiding ParameterHandler and other

Re: [deal.II] Re: Issue on CUDA code

wrong. > > Best, > > Bruno > > Le lun. 17 juil. 2023 à 11:06, giuseppe orlando a > écrit : > >> Here, you can find the two codes. 'SBR_GPU' is the one that arrives up to >> the end, whereas 'SBR_GPU_WRONG' is the one that stops with the error >

Re: [deal.II] Re: Issue on CUDA code

> The new code is not 100% backward compatible but the changes are pretty > small. I advise you to switch to the latest version of deal.II, debug the > code on the CPU, and then run on the GPU. > > Best, > > Bruno > > Le lun. 17 juil. 2023 à 14:12, giuseppe orlando a >

[deal.II] Compilation deal.II 9.5

Hello everyone, I'm writing this post because I'm having some issues to compile the latest release. I would need a version which supports CUDA and MPI being CUDA-aware. Hence, I installed Kokkos using the flags -DKokkos_ENABLE_CUDA=ON and -DKokkos_ENABLE_CUDA_LAMBDA=ON. Then, I tried to build t

[deal.II] Re: Compilation deal.II 9.5

Dear Bruno, thanks for the suggestion. I tried to follow it, but I still the get same error during the creation of the Makefile. Just to sum up I typed "export OMPI_CXX=/project/sm77/bin/kokkos/bin/nvcc_wrapper" "cmake -DCMAKE_INSTALL_PREFIX=/project/sm77/bin/dealII/9.5.1/ -DCMAKE_CXX_COMPILER=/p

Re: [deal.II] Re: Compilation deal.II 9.5

> it might be passed to MPI when it shouldn't. The rest looks good. > > Best, > > Bruno > > > Le mer. 19 juil. 2023 à 13:29, giuseppe orlando a > écrit : > >> Dear Bruno, >> thanks for the suggestion. I tried to follow it, but I still the get same

R: [deal.II] Re: Compilation deal.II 9.5

Yes, I’m inside a clean directory. The output of mpicxx is “/opt/cray/pe/craype/2.7.10/bin/CC”, which is the same that I used to install all the other versions. The cmake version is 3.17.0. Best Giuseppe Orlando Da: Bruno TurcksinInviato: giovedì 20 luglio 2023 17:47A: dealii

Re: [deal.II] Re: Issue on CUDA code

there is an underlying issue in handling periodic bcs in a matrix-free framework for CUDA. Thanks in advance Best regards, Giuseppe Il giorno martedì 18 luglio 2023 alle 17:31:17 UTC+2 giuseppe orlando ha scritto: > Dear Bruno, > thank you very much. I am having some troubles to

[deal.II] Different basis functions and polynomial degrees for each space direction

Hello everyone, I'm writing this post because I would need help to build a "particular" finite element space. Let us suppose to have a scalar field in 2D. I want to consider a family of basis functions along the horizontal direction (e.g. basis functions based on Legendre polynomials) and a diff

Re: [deal.II] Different basis functions and polynomial degrees for each space direction

Dear Wolfgang, thanks for the kind response. As a starting point, I would need to focus on a scalar field, hence a finite element space for scalar values would be sufficient, but obviously the mid-long term goal is to employ vector-value finite element as well. I will try to have a look to the c

Re: [deal.II] Spherical coordinates output

Dear Wolfgang, thanks a lot for you answer. My goal would be able to reproduce images like these ones, which are a cross-section at a fixed latitude. Hence, I coulde use 'slice' type functions if I have spherical coordinates, but not with Cartesian ones. Best, Giuseppe [image: Details are in th

[deal.II] Spherical coordinates output

ent/doxygen/deal.II/namespaceGridGenerator.html#ab2dc2f16fc52cf1e51b9afdc4264ac70> ? Thanks in advance. Best, Giuseppe Orlando -- 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

[deal.II] Re: Supported finite elements in the matrix-free infrastrucutre

Dear Peter, thanks a lot for your answer. Best, Giuseppe Il giorno lunedì 3 marzo 2025 alle 19:53:58 UTC+1 peterr...@gmail.com ha scritto: > The normal RT element is not working. But FE_RaviartThomasNodal is working > with MF. > > Regarding simplices: `FE_Simplex(DG)P` (also with FESystem) is