Re: [deal.II] make hanging node constraint for locally refined mesh

2020-05-26 Thread Phạm Ngọc Kiên
Dear Prof. Wolfgang Bangerth,

Thank you very much.

I will change the codes as your guidance.

Yours sincerely,
Kien

Vào Th 4, 27 thg 5, 2020 vào lúc 12:58 Wolfgang Bangerth <
bange...@colostate.edu> đã viết:

> On 5/26/20 7:30 PM, Phạm Ngọc Kiên wrote:
> >
> > Thus, the conforming FE_Nedelec and FE_NedelecSZ can only work with the
> > conforming mesh. The non-conforming mesh, which has hanging nodes, does
> not
> > work with these elements.
>
> The FE_NedelecSZ element doesnt, but FE_Nedelec does work with hanging
> nodes.
>
> Best
>   W.
>
>
> --
> 
> Wolfgang Bangerth  email: bange...@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/44ffc21e-c3e2-4f48-a4ec-1537f37e8b42%40colostate.edu
> .
>

-- 
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/CAAo%2BSZcyBae4j_-0ahdsQSrVCbTCrrVOHpixTGiraB8_zONUjg%40mail.gmail.com.


Re: [deal.II] make hanging node constraint for locally refined mesh

2020-05-26 Thread Wolfgang Bangerth

On 5/26/20 7:30 PM, Phạm Ngọc Kiên wrote:


Thus, the conforming FE_Nedelec and FE_NedelecSZ can only work with the 
conforming mesh. The non-conforming mesh, which has hanging nodes, does not 
work with these elements.


The FE_NedelecSZ element doesnt, but FE_Nedelec does work with hanging nodes.

Best
 W.


--

Wolfgang Bangerth  email: bange...@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/44ffc21e-c3e2-4f48-a4ec-1537f37e8b42%40colostate.edu.


Re: [deal.II] table.h errors

2020-05-26 Thread Jammi Ashok
Actually I cloned deal-ii from github. Day before yesterday.

-- 
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/CANsJfM33BOUthrRfM3krHH4edyGmGXHs39MB8uedAdMDKX3pdg%40mail.gmail.com.


Re: [deal.II] table.h errors

2020-05-26 Thread Jammi Ashok
Visual studio - 2017 and deal-ii-9.1.1.

On Wed, May 27, 2020 at 2:13 AM Wolfgang Bangerth 
wrote:

> On 5/26/20 12:52 AM, Jammi Ashok wrote:
> > I am getting error in the following lines in "table.h".
> >
> > // Make the underlying iterator class a friend.
> > 958 friend class LinearIndexIterator
> > <
> https://nam01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwww.dealii.org%2Fdeveloper%2Fdoxygen%2Fdeal.II%2FclassLinearIndexIterator.html=02%7C01%7CWolfgang.Bangerth%40colostate.edu%7C07d8ed42c7934fc9907508d801414cd0%7Cafb58802ff7a4bb1ab21367ff2ecfc8b%7C0%7C0%7C637260727265355249=gOS6uSo%2BhH9AoOdYF5sMvholtVS7I%2FDlH6%2BoKnznpQY%3D=0
> ><
> > 959 Iterator
> > <
> https://nam01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwww.dealii.org%2Fdeveloper%2Fdoxygen%2Fdeal.II%2FclassMatrixTableIterators_1_1Iterator.html=02%7C01%7CWolfgang.Bangerth%40colostate.edu%7C07d8ed42c7934fc9907508d801414cd0%7Cafb58802ff7a4bb1ab21367ff2ecfc8b%7C0%7C0%7C637260727265365241=GpBlt90%2FL2VnfXI8UEtlnfnmrTrr7aP%2FZU6CCYZw4go%3D=0>
> > Constness, storage_order>,
> > 960
> > <
> https://nam01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwww.dealii.org%2Fdeveloper%2Fdoxygen%2Fdeal.II%2FclassMatrixTableIterators_1_1AccessorBase.html%23aea52115385a3e007cab612d2977f2eb9=02%7C01%7CWolfgang.Bangerth%40colostate.edu%7C07d8ed42c7934fc9907508d801414cd0%7Cafb58802ff7a4bb1ab21367ff2ecfc8b%7C0%7C0%7C637260727265365241=hZh%2BJuMaUdhDAkxEw3WzslQIvrH91FZL7nCxtY%2FDPZk%3D=0>
>
> > Accessor
> > <
> https://nam01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwww.dealii.org%2Fdeveloper%2Fdoxygen%2Fdeal.II%2FclassMatrixTableIterators_1_1Accessor.html=02%7C01%7CWolfgang.Bangerth%40colostate.edu%7C07d8ed42c7934fc9907508d801414cd0%7Cafb58802ff7a4bb1ab21367ff2ecfc8b%7C0%7C0%7C637260727265375236=QxqvRa6XbHCI9wrGydH7mN4KmPgESE6ZYjIYJyJT5A4%3D=0>
> > Constness, storage_order>>;
> > 961  };
> >
> >
> > I am using windows x64 machine.
> > THe error says "C2990 'dealii::LinearIndexIterator': non-class template
> has
> > already been declared as a class template"
> > Also,
> > "syntax error '<'.
> > "unexpected tokens preceding ';'
>
> Jammi,
> what version of the Microsoft compiler are you using? And what version of
> deal.II?
>
> The most recent deal.II versions all work just fine with the more recent
> Visual Studio compilers, so I am surprised that you are getting this error.
>
> Best
>   W.
>
>
> --
> 
> Wolfgang Bangerth  email: bange...@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/51d86946-09f8-b778-22ea-d93f204542d1%40colostate.edu
> .
>

-- 
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/CANsJfM0DFEn9x-Tmuf%2BsfJqzOYh8Snwz6gCsA1D8SDiu_t04TQ%40mail.gmail.com.


Re: [deal.II] make hanging node constraint for locally refined mesh

2020-05-26 Thread Phạm Ngọc Kiên
Dear Prof. Wolfgang Bangerth,

I used FE_NedelecSZ element.

Thus, the conforming FE_Nedelec and FE_NedelecSZ can only work with the
conforming mesh. The non-conforming mesh, which has hanging nodes, does not
work with these elements.

Thank you very much.

Yours sincerely,
Kien

Vào Th 4, 27 thg 5, 2020 vào lúc 05:42 Wolfgang Bangerth <
bange...@colostate.edu> đã viết:

> On 5/26/20 1:34 AM, Phạm Ngọc Kiên wrote:
> > /*  const dealii::FullMatrix& dealii::FiniteElement<,
> >  >::constraints(const dealii::internal::SubfaceCase&)
> const
> > [with int dim = 3; int spacedim = 3]
> > The violated condition was:
> >  (this->dofs_per_face == 0) || (interface_constraints.m() != 0)
> > Additional information:
> >  The finite element for which you try to obtain hanging node
> constraints
> > does not appear to implement them.*/
> > /*
> > */
> > Does this error comes from the AffineConstraints when I try to implement
> the
> > DoFTools::make_hanging_node_constraints?
> >
> > Could you please help me to overcome this error?
>
> Well, the error message says pretty clearly what is happening: The element
> you
> use does not provide hanging node information. What element are you using?
>
> Best
>   W.
>
> --
> 
> Wolfgang Bangerth  email: bange...@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/fc63875b-df35-eb80-02cd-c503ca81d82e%40colostate.edu
> .
>

-- 
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/CAAo%2BSZfXpVTE1NW1gMpbsSsng08KXFaJfHqO%2B%2BUVDiMoqFjGXA%40mail.gmail.com.


[deal.II] Re: Can dealii handle different type of element in a grid

2020-05-26 Thread 徐辉
Thanks sir, I will have a look at the tutorial you listed. Thank you for 
you all's time.

-- 
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/683c7555-daa9-4cca-a4d2-411ab0126bf3%40googlegroups.com.


Re: [deal.II] table.h errors

2020-05-26 Thread Wolfgang Bangerth

On 5/26/20 12:52 AM, Jammi Ashok wrote:

I am getting error in the following lines in "table.h".

// Make the underlying iterator class a friend.
958 friend class LinearIndexIterator 
<
959 Iterator 
Constness, storage_order>,
960 
 
Accessor 
Constness, storage_order>>;

961  };


I am using windows x64 machine.
THe error says "C2990 'dealii::LinearIndexIterator': non-class template has 
already been declared as a class template"

Also,
"syntax error '<'.
"unexpected tokens preceding ';'


Jammi,
what version of the Microsoft compiler are you using? And what version of 
deal.II?

The most recent deal.II versions all work just fine with the more recent 
Visual Studio compilers, so I am surprised that you are getting this error.


Best
 W.


--

Wolfgang Bangerth  email: bange...@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/51d86946-09f8-b778-22ea-d93f204542d1%40colostate.edu.


Re: [deal.II] make hanging node constraint for locally refined mesh

2020-05-26 Thread Wolfgang Bangerth

On 5/26/20 1:34 AM, Phạm Ngọc Kiên wrote:
/*  const dealii::FullMatrix& dealii::FiniteElement<, 
 >::constraints(const dealii::internal::SubfaceCase&) const 
[with int dim = 3; int spacedim = 3]

The violated condition was:
     (this->dofs_per_face == 0) || (interface_constraints.m() != 0)
Additional information:
     The finite element for which you try to obtain hanging node constraints 
does not appear to implement them.*/

/*
*/
Does this error comes from the AffineConstraints when I try to implement the 
DoFTools::make_hanging_node_constraints?


Could you please help me to overcome this error?


Well, the error message says pretty clearly what is happening: The element you 
use does not provide hanging node information. What element are you using?


Best
 W.

--

Wolfgang Bangerth  email: bange...@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/fc63875b-df35-eb80-02cd-c503ca81d82e%40colostate.edu.


Re: [deal.II] vertex and dofcells

2020-05-26 Thread Jean-Paul Pelteret
Dear Alberto,

Quick answer: There’s a entry in the FAQ that may answer your question:
https://github.com/dealii/dealii/wiki/Frequently-Asked-Questions#how-do-i-get-the-degree-of-freedom-indices-at-vertices
 


As mentioned at the bottom of that entry, is a small point on how to go from 
DoFs -> vertices using DoFTools::map_dofs_to_support_points() 
.

Does this help you?
Best,
Jean-Paul

> On 26 May 2020, at 18:35, Alberto Salvadori  
> wrote:
> 
> Dear community,
> 
> this looks a simple question but I just cannot find a simple answer, I 
> apologize.
> 
> The code I aim at writing is similar to:
> 
> unsigned vv = 0 ;
> 
> 
> for (unsigned int ii=0; ii 
> {
> 
>   const unsigned int ii_group = 
> manifold_fe.system_to_base_index(ii).first.first;
> 
>   
>   if ( ii_group == manifold_dofs.cS_dof )
> 
>   {
> 
> Point this_node = cell->vertex(vv) ;
> 
> 
> [ ... omitted code ... ]
> 
> 
> 
> cell_matrix( ii,ii ) = 1 ;
> 
> cell_rhs( ii ) = some_function_of( this_node );
> 
> 
> 
> ...
> 
> ++vv;
> 
>   }
> 
>  }
> 
> 
> 
> In a nutshell, I want to retrieve information that I stored at each node of 
> the cell, while looping on the dofcells. 
> I am not sure if the order of the dofs and of the vertex coincide, suspect 
> they do not.
> Can you please address me to a function that connects the dofcells to the 
> associated vertex,or suggest a better approach?
> 
> Thank you so much
> 
> Alberto
> 
> 
> Informativa sulla Privacy: http://www.unibs.it/node/8155 
> 
> 
> -- 
> 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/3a275ee6-9416-4547-8a81-188c29579e2a%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/F28279CC-8C0A-4CFC-84AC-4D97BE725C2E%40gmail.com.


[deal.II] vertex and dofcells

2020-05-26 Thread Alberto Salvadori
Dear community,

this looks a simple question but I just cannot find a simple answer, I 
apologize.

The code I aim at writing is similar to:

*unsigned* vv = 0 ;



*for* (*unsigned* *int* ii=0; ii this_node = cell->vertex(vv) ;



[ ... omitted code ... ]


cell_matrix( ii,ii ) = 1 ;

cell_rhs( ii ) = some_function_of( this_node );


...

++vv;

  }

 }



In a nutshell, I want to retrieve information that I stored at each node of 
the cell, while looping on the dofcells. 
I am not sure if the order of the dofs and of the vertex coincide, suspect 
they do not.
Can you please address me to a function that connects the dofcells to the 
associated vertex,or suggest a better approach?

Thank you so much

Alberto

-- 


Informativa sulla Privacy: http://www.unibs.it/node/8155 


-- 
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/3a275ee6-9416-4547-8a81-188c29579e2a%40googlegroups.com.


Re: [deal.II] Can dealii handle different type of element in a grid

2020-05-26 Thread Wolfgang Bangerth

On 5/26/20 8:13 AM, 徐辉 wrote:
Or another example in 2D, I have some some one order quadrilateral element and 
some two order quadrilateral element.


Or, if you mean the polynomial degree, you want to use the hp framework -- see 
step-27 and step-46.


Best
 W.

--

Wolfgang Bangerth  email: bange...@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/c9046fd7-dd46-d29d-092f-39901b3a3728%40colostate.edu.


[deal.II] Re: Can dealii handle different type of element in a grid

2020-05-26 Thread Bruno Turcksin
Hi,

deal.II only supports quads in 2D and hex in 3D.So if you have triangles 
you will need to split them. For the quads, you probably want to load order 
one quads and then, use a manifold to get the higher order. step 49 
(https://dealii.org/current/doxygen/deal.II/step_49.html)  shows how to 
manipulate meshes.

Best,

Bruno

On Tuesday, May 26, 2020 at 10:13:35 AM UTC-4, 徐辉 wrote:
>
> I want using different type of element in a grid in dealii.
>
> I already have a grid,  I want solve problem based on this grid.  The grid 
> contain different type of element. For example in 2D, I have some one order 
> triangle element and some one order quadrilateral element. Or another 
> example in 2D, I have some some one order quadrilateral element and some 
> two order quadrilateral element.
>
> I search dealii's tutorial do not found any topic talk about using 
> different type element in a grid. And I found in dealii FAQ first item that 
> seems dealii not support triangle or tetrahedra.
>
> Does dealii can only support one kind of element to assemble a discrete 
> system? All element in grid must be same? Is there any reference can help 
> me?
>
> Thanks for your time.
>

-- 
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/175d5891-8b6f-4665-b895-b3a32101f8fa%40googlegroups.com.


[deal.II] Can dealii handle different type of element in a grid

2020-05-26 Thread 徐辉
I want using different type of element in a grid in dealii.

I already have a grid,  I want solve problem based on this grid.  The grid 
contain different type of element. For example in 2D, I have some one order 
triangle element and some one order quadrilateral element. Or another 
example in 2D, I have some some one order quadrilateral element and some 
two order quadrilateral element.

I search dealii's tutorial do not found any topic talk about using 
different type element in a grid. And I found in dealii FAQ first item that 
seems dealii not support triangle or tetrahedra.

Does dealii can only support one kind of element to assemble a discrete 
system? All element in grid must be same? Is there any reference can help 
me?

Thanks for your time.

-- 
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/9be66e02-fa6e-4d47-b012-a80c76a930d3%40googlegroups.com.


Re: [deal.II] Blog post about how we develop deal.II

2020-05-26 Thread Thomas WICK

Hi Wolfgang,



Am 26.05.20 um 04:52 schrieb Wolfgang Bangerth:


Well, not really about how *we* develop deal.II -- but at least about 
my approach in general to developing software:
https://bssw.io/blog_posts/making-myself-better-what-craftspeople-can-teach-us-about-software 

I think that what I describe is how a lot of good programmers think 
about their craft. 


I like that.

Best Thomas

--
---

Prof. Dr. Thomas Wick
Leibniz Universität Hannover (LUH)
Institut für Angewandte Mathematik (IfAM)
Arbeitsgruppe Wissenschaftliches Rechnen (GWR)

Welfengarten 1
30167 Hannover, Germany

Tel.:   +49 511 762 3360
Email:  thomas.w...@ifam.uni-hannover.de
www:http://www.ifam.uni-hannover.de/wick
www:http://www.thomaswick.org

---
--




If you're still learning how to be a programmer, this might be of 
interest to you.


See many of you tomorrow for the deal.II workshop!
Best
 Wolfgang



--
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/0667ff0e-736b-078f-a6c4-03f95237c4ea%40gmail.com.


[deal.II] make hanging node constraint for locally refined mesh

2020-05-26 Thread Phạm Ngọc Kiên
Dear all,
When I try to create a locally refined mesh by the following codes:











* GridGenerator::hyper_cube(triangulation, 0, 1);Point 
center (0.5, 0.5, 0.5);triangulation.refine_global(5);
for (auto cell : triangulation.active_cell_iterators()) {
if (cell->is_locally_owned()) {if 
(cell->point_inside(center)){
cell->set_refine_flag();}}
}triangulation.execute_coarsening_and_refinement();*

The library return this error 




*  const dealii::FullMatrix& dealii::FiniteElement<, 
 >::constraints(const dealii::internal::SubfaceCase&) const 
[with int dim = 3; int spacedim = 3]The violated condition was: 
(this->dofs_per_face == 0) || (interface_constraints.m() != 0)Additional 
information: The finite element for which you try to obtain hanging 
node constraints does not appear to implement them.*

Does this error comes from the AffineConstraints when I try to implement 
the DoFTools::make_hanging_node_constraints?

Could you please help me to overcome this error?

Thank you very much.

Best regards,
Kien

-- 
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/fdcc4fd8-761d-4cc4-a954-8631e0cf117d%40googlegroups.com.


[deal.II] table.h errors

2020-05-26 Thread Jammi Ashok
Dear sir,
I am getting error in the following lines in "table.h".

// Make the underlying iterator class a friend.
958  friend class LinearIndexIterator 

<
959  Iterator 
,
960 

  
Accessor 
>;
961  };


I am using windows x64 machine.
THe error says "C2990 'dealii::LinearIndexIterator': non-class template has 
already been declared as a class template"
Also,
"syntax error '<'.
"unexpected tokens preceding ';'


Kindly help me out!!!

Ashok

-- 
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/d882ab60-8f4f-4bbb-8258-818c4d88e790%40googlegroups.com.