Re: [deal.II] Issue with PETSc vector

2023-12-20 Thread Giselle Sosa Jones
13:23, Giselle Sosa Jones wrote: > > // Initialize vector > > dof_handler_dg0.distribute_dofs(fe_dg0); > > const std::vector locally_owned_dofs_per_proc_dg0 = > > DoFTools::locally_owned_dofs_per_subdomain(dof_handler_dg0); > > locally_owned_dofs_dg0 = > locally_o

[deal.II] Issue with PETSc vector

2023-12-18 Thread Giselle Sosa Jones
Hello, I have an issue that seems very simple but I cannot figure out. I have a piecewise constant parallel vector that I defined using a DG-0 space. I have an auxiliary function that calculates the value according to the cell, and then in my main file, I would like to assign the value to the

Re: [deal.II] Configuration error

2023-10-04 Thread Giselle Sosa Jones
try to compile the deal.II version 9.5.1. > > I hope this gives sufficient information, and we are able to identify the > problem (probably either with PETSc or the linking of libraries). > > Thank you in advance! > > Best, > Florian > > > Am 04.10.2023

Re: [deal.II] Configuration error

2023-10-03 Thread Giselle Sosa Jones
I do not have this file in my CMakeFiles directory. I was looking for it in the other directories, but could not find it... Thank you for your help. On Tue, Oct 3, 2023 at 8:12 PM Wolfgang Bangerth wrote: > On 10/3/23 16:15, Giselle Sosa Jones wrote: > > > > I am trying to co

[deal.II] Configuration error

2023-10-03 Thread Giselle Sosa Jones
Hello, I am trying to configure deal II with options: cmake -DDEAL_II_WITH_PETSC=ON -DCMAKE_INSTALL_PREFIX=/Users/gisellesosajones/Documents/deal_ii_95 -DPETSC_DIR=$PETSC_DIR -DPETSC_ARCH=$PETSC_ARCH -DDEAL_II_WITH_MPI=ON -DDEAL_II_WITH_HDF5=OFF I am getting this error CMake Error at

[deal.II] Install deal II with PETSc and MPICH

2023-09-01 Thread Giselle Sosa Jones
Hello, I have a PETSc installation that uses mpich. I am trying to install deal II with: cmake -DDEAL_II_WITH_PETSC=ON -DCMAKE_INSTALL_PREFIX=/Users/gisellesosajones/Documents/deal_ii_95 -DPETSC_DIR=$PETSC_DIR -DPETSC_ARCH=$PETSC_ARCH -DDEAL_II_WITH_MPICH=ON -DDEAL_II_WITH_HDF5=OFF

Re: [deal.II] Tensor product polynomials

2022-06-22 Thread Giselle Sosa Jones
o make sure you use the same compilers both > for your code and the library. In practice, though, the ABI change occurred > in GCC 5 so combining any two newer versions is typically okay. > > Best, > David > -- > *From:* dealii@googlegroups.com on beh

Re: [deal.II] Tensor product polynomials

2022-05-31 Thread Giselle Sosa Jones
d from FE_DGRaviartThomas). I am not sure if this is a C++ issue and/or if there is an easier way of doing this. Any help is highly appreciated! Best, Giselle Sosa Jones On Tuesday, August 31, 2021 at 9:59:55 a.m. UTC-5 Giselle Sosa Jones wrote: > Thank you, Dr. Bangerth, I think these extra

Re: [deal.II] Tensor product polynomials

2021-08-31 Thread Giselle Sosa Jones
Thank you, Dr. Bangerth, I think these extra comments clarify it for me. Best, Giselle On Monday, August 30, 2021 at 8:15:16 p.m. UTC-5 Wolfgang Bangerth wrote: > On 8/28/21 4:52 PM, Giselle Sosa Jones wrote: > > > > I need to create a local space of polynomials in 3D of the for

[deal.II] Tensor product polynomials

2021-08-28 Thread Giselle Sosa Jones
Hello, I need to create a local space of polynomials in 3D of the form Q_{k-1,k,k} \times Q_{k,k-1,k} \times Q_{k,k,k-1}. This space is very similar to the DG Raviart-Thomas space, with the difference that the polynomial order of each component is one order lower in the corresponding

[deal.II] Local Raviart-Thomas space

2021-06-16 Thread Giselle Sosa Jones
Hello, I am trying to implement a projection of a DG solution to a Raviart-Thomas space. For this, I need a local space Q_{p-1,p}(K) \times Q_{p,p-1}(K). I have been using FE_DGRaviartThomas of degree p-1, and then I use shape_value_component(i, q, 1) when I want to use the x component of the

Re: [deal.II] DG weighted average

2021-03-12 Thread Giselle Sosa Jones
Thank you Dr. Bangerth, this is exactly what I was missing. Best, Giselle On Thursday, March 11, 2021 at 3:36:53 p.m. UTC-6 Wolfgang Bangerth wrote: > On 3/11/21 10:33 PM, Giselle Sosa Jones wrote: > > > > I have a diffusion problem where the diffusion coeffi

[deal.II] DG weighted average

2021-03-11 Thread Giselle Sosa Jones
Hello, I have a diffusion problem where the diffusion coefficient K depends on the solution at the previous time step and is therefore discontinuous. I am trying to implement a weighted average for DG, which is similar to the standard average, with the difference that the weights are not 0.5,

Re: [deal.II] Re: Configuring deal.II with LAPACK

2021-02-20 Thread Giselle Sosa Jones
Hi Bruno, Thank you for your reply. I tried this but got the same error as before. I am not sure what's happening since I am using PETSc with Intel mkl in other software, and it works fine there. Do you happen to have any other suggestions? Maybe it is an issue with Intel mkl? Thanks again.