Re: [deal.II] Tensor product polynomials

2022-06-22 Thread Wells, David
Hi Giselle,

What do you mean by "the hh file is not recognized in the cc file" - is there a 
specific error message this refers to?

In general, if you need to implement a new finite element, you will need to 
write some code. Can you compile any program that depends on deal.II right now?

Best,
David

From: dealii@googlegroups.com  on behalf of Giselle 
Sosa Jones 
Sent: Wednesday, June 22, 2022 10:45 AM
To: dealii@googlegroups.com 
Subject: Re: [deal.II] Tensor product polynomials

Hi David,

I've tried updating my C++ compiler and recompiling deal.II, but I have the 
same error. It is as if the hh file is not recognized in the cc file. I know 
this is a long shot, but do you happen to have any other suggestions? Maybe 
there is an easier way to create this space that I need?

Thank you!

Best,
Giselle

On Tue, May 31, 2022 at 4:28 PM Wells, David 
mailto:drwe...@email.unc.edu>> wrote:
Hi Gisselle,

This error message occurs when you use a very old (probably GCC 4.8 or 4.9) C++ 
compiler with a version of deal.II built with a newer C++ compiler (or vice 
versa). A sure fix is to 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 
mailto:dealii@googlegroups.com>> on behalf of Giselle 
Sosa Jones mailto:giselle1...@gmail.com>>
Sent: Tuesday, May 31, 2022 11:55 AM
To: deal.II User Group mailto:dealii@googlegroups.com>>
Subject: Re: [deal.II] Tensor product polynomials

You don't often get email from 
giselle1...@gmail.com. Learn why this is 
important
Hello,

I sent this message a while ago and things became clearer for me, but I am 
still not able to create the spaces that I want. My goal is to create local 
spaces that look like this: Q_{k-1,k,k} \times Q_{k,k-1,k} \times Q_{k,k,k-1}, 
for k \geq 1.This is what I am currently doing:

- Created my own class of polynomials similar to the deal.II class 
PolynomialsRaviartThomas, but with the polynomial degrees as described above.
- Created my own class (FE_RTProj) for a finite element space similar to the 
one called FE_DGRaviartThomas. This is a derived class from FE_DGVector and the 
polynomial type is the one created in the item above.

My code does not even compile since I get error messages like: error: undefined 
reference to 'FE_RTProj::FE_RTProj<2, 2>::get_name[abi:cxx11]() const. The 
function get_name is defined in the FE_RTProj h and cc files (I basically 
copied 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 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 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 direction. I was looking at the file
> polynomials_raviart_thomas.cc to get an idea of how to create them but I
> am a bit lost. Can someone offer a suggestion?

Giselle -- can you be specific about which part you have trouble with?
I've tried to improve the documentation of the key function of that
class a bit:
https://github.com/dealii/dealii/pull/12726

Best
Wolfgang

--

Wolfgang Bangerth email: bang...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/

--
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/f3388a43-fe77-4c01-a346-1a7fb89d243dn%40googlegroups.com.

--
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 

[deal.II] deal.II Newsletter #216

2022-06-22 Thread 'Rene Gassmoeller' via deal.II User Group
Hello everyone!

This is deal.II newsletter #216.
It automatically reports recently merged features and discussions about the 
deal.II finite element library.


## Below you find a list of recently proposed or merged features:

#14033: [9.4] Takeover #14032: Disable HYPRE in the Tvmult test. (proposed by 
drwells) https://github.com/dealii/dealii/pull/14033

#14032: Disable HYPRE in the Tvmult test. (proposed by drwells) 
https://github.com/dealii/dealii/pull/14032

#14031: Remove deprecated atomicAdd_wrapper (proposed by masterleinad) 
https://github.com/dealii/dealii/pull/14031

#14030: [9.4] Take over #14029: Work around a GCC-12 warning about array 
accesses. (proposed by drwells; merged) 
https://github.com/dealii/dealii/pull/14030

#14029: Work around a GCC-12 warning about array accesses. (proposed by 
drwells; merged) https://github.com/dealii/dealii/pull/14029

#14028: [9.4] Take over #14015: Make get_cell_range_category() more consistent 
with get_cell_category() (proposed by peterrum; merged) 
https://github.com/dealii/dealii/pull/14028

#14027: [9.4] Take over #14026: Fix warning of gcc-12 with gather_pd/ps 
instructions (proposed by tamiko; merged) 
https://github.com/dealii/dealii/pull/14027

#14026: Fix warning of gcc-12 with gather_pd/ps instructions (proposed by 
kronbichler; merged) https://github.com/dealii/dealii/pull/14026

#14025: Query ghost index set only once in Partitioner::global_to_local 
(proposed by kronbichler; merged) https://github.com/dealii/dealii/pull/14025

#14024: Do not add locally owned dofs to ghost set during construction 
(proposed by kronbichler) https://github.com/dealii/dealii/pull/14024

#14023: Remove obsolete header file (proposed by kronbichler; merged) 
https://github.com/dealii/dealii/pull/14023

#14022: Reduce header inclusions of index_set.h in evaluation_kernels (proposed 
by kronbichler) https://github.com/dealii/dealii/pull/14022

#14021: Restructure inline functions of IndexSet (proposed by kronbichler) 
https://github.com/dealii/dealii/pull/14021

#14020: Avoid unnecessary locally_relevant_dofs in mg_transfer_internal 
(proposed by kronbichler) https://github.com/dealii/dealii/pull/14020

#14019: MGConstrainedDoFs::clear_user_constraints (proposed by zjiaqi2018; 
merged) https://github.com/dealii/dealii/pull/14019

#14018: Better describe what FullMatrix is and does. (proposed by bangerth; 
merged) https://github.com/dealii/dealii/pull/14018

#14017: Minor text improvements in a tutorial documentation. (proposed by 
bangerth; merged) https://github.com/dealii/dealii/pull/14017

#14016: Add reference to fast hanging-node algo and global coarsening (proposed 
by peterrum; merged) https://github.com/dealii/dealii/pull/14016

#14015: Make get_cell_range_category() more consistent with get_cell_category() 
(proposed by peterrum; merged) https://github.com/dealii/dealii/pull/14015

#14014: Add test (MG+DG+hp) (proposed by peterrum) 
https://github.com/dealii/dealii/pull/14014

#14013: MatrixFree: allow categorization of ghost cells (proposed by peterrum; 
merged) https://github.com/dealii/dealii/pull/14013

#14012: [9.4] Take over #14010 and #14011 (proposed by peterrum; merged) 
https://github.com/dealii/dealii/pull/14012

#14011: Can't use ReferenceCells in switch-case. (proposed by bangerth; merged) 
https://github.com/dealii/dealii/pull/14011

#14010: Can't use ReferenceCells in switch-case. (proposed by bangerth; merged) 
https://github.com/dealii/dealii/pull/14010

#14009: [9.4] Fix compilation with clang-13. (proposed by drwells; merged) 
https://github.com/dealii/dealii/pull/14009

#14007: Improve documentation of MF::AD::cell_vectorization_category (proposed 
by peterrum; merged) https://github.com/dealii/dealii/pull/14007

#14006: Fix compilation with clang-13. (proposed by drwells; merged) 
https://github.com/dealii/dealii/pull/14006

#14005: MatrixFree: improve assert (proposed by peterrum; merged) 
https://github.com/dealii/dealii/pull/14005

#14003: [9.4] Take over #14002, Fix failing symengine tests (proposed by 
tamiko; merged) https://github.com/dealii/dealii/pull/14003

#14002: Fix failing symengine tests (proposed by tamiko; merged) 
https://github.com/dealii/dealii/pull/14002

#14001: Update AUTHORS.md (proposed by NiklasWik; merged) 
https://github.com/dealii/dealii/pull/14001

#13999: [9.4] CMake: Remove taskflow configure (proposed by tamiko; merged) 
https://github.com/dealii/dealii/pull/13999

#13998: [9.4] remove p4est_bug_01 test (proposed by tamiko; merged) 
https://github.com/dealii/dealii/pull/13998

#13995: switch all DEAL_II_DEPRECATED_EARLY to DEAL_II_DEPRECATED (proposed by 
tamiko; merged) https://github.com/dealii/dealii/pull/13995

#13994: update VERSION (proposed by tamiko; merged) 
https://github.com/dealii/dealii/pull/13994

#13993: Update tag file references for release. (proposed by tamiko; merged) 
https://github.com/dealii/dealii/pull/13993

#13992: Update changelog for release (proposed by tamiko; merged) 

Re: [deal.II] Tensor product polynomials

2022-06-22 Thread Giselle Sosa Jones
Hi David,

I've tried updating my C++ compiler and recompiling deal.II, but I have the
same error. It is as if the hh file is not recognized in the cc file. I
know this is a long shot, but do you happen to have any other suggestions?
Maybe there is an easier way to create this space that I need?

Thank you!

Best,
Giselle

On Tue, May 31, 2022 at 4:28 PM Wells, David  wrote:

> Hi Gisselle,
>
> This error message occurs when you use a very old (probably GCC 4.8 or
> 4.9) C++ compiler with a version of deal.II built with a newer C++ compiler
> (or vice versa). A sure fix is to 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 behalf of
> Giselle Sosa Jones 
> *Sent:* Tuesday, May 31, 2022 11:55 AM
> *To:* deal.II User Group 
> *Subject:* Re: [deal.II] Tensor product polynomials
>
> You don't often get email from giselle1...@gmail.com. Learn why this is
> important 
> Hello,
>
> I sent this message a while ago and things became clearer for me, but I am
> still not able to create the spaces that I want. My goal is to create local
> spaces that look like this: Q_{k-1,k,k} \times Q_{k,k-1,k} \times
> Q_{k,k,k-1}, for k \geq 1.This is what I am currently doing:
>
> - Created my own class of polynomials similar to the deal.II
> class PolynomialsRaviartThomas, but with the polynomial degrees as
> described above.
> - Created my own class (FE_RTProj) for a finite element space similar to
> the one called FE_DGRaviartThomas. This is a derived class from FE_DGVector
> and the polynomial type is the one created in the item above.
>
> My code does not even compile since I get error messages like: error:
> undefined reference to 'FE_RTProj::FE_RTProj<2, 2>::get_name[abi:cxx11]()
> const. The function get_name is defined in the FE_RTProj h and cc files (I
> basically copied 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 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 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 direction. I was looking at the file
> > polynomials_raviart_thomas.cc to get an idea of how to create them but I
> > am a bit lost. Can someone offer a suggestion?
>
> Giselle -- can you be specific about which part you have trouble with?
> I've tried to improve the documentation of the key function of that
> class a bit:
> https://github.com/dealii/dealii/pull/12726
>
> Best
> Wolfgang
>
> --
> 
> Wolfgang Bangerth email: bang...@colostate.edu
> www: http://www.math.colostate.edu/~bangerth/
>
> --
> 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/f3388a43-fe77-4c01-a346-1a7fb89d243dn%40googlegroups.com
> 
> .
>
> --
> 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/MN2PR03MB4879589A5424DA0E1259B82DEDDC9%40MN2PR03MB4879.namprd03.prod.outlook.com
> 
> .
>

-- 
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