[deal.II] Re: Calculation of local flux matrix in DG: looping over a cell's faces

2019-09-23 Thread vachan potluri
Found the face ordering here https://www.dealii.org/current/doxygen/deal.II/structGeometryInfo.html -- 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 subsc

Re: [deal.II] How to get the coordinates of dof with FE_Raviart_Thomas?

2019-09-23 Thread Praveen C
If you are going to use Cartesian meshes only, then you may want to check out https://dealii.org/developer/doxygen/deal.II/classFE__RaviartThomasNodal.html Best praveen > On 24-Sep-2019, at 10:42 AM, Phạm Ngọc Kiên

[deal.II] Calculation of local flux matrix in DG: looping over a cell's faces

2019-09-23 Thread vachan potluri
Hello all, I want to calculate local flux matrix of a cell (in 2D). The algorithm I thought of is the following: [image: CodeCogsEqn.png] Within loop over cells: start loop over faces: re-initialize fe_face_values for current face calculate the 1D mass matrix add

Re: [deal.II] How to get the coordinates of dof with FE_Raviart_Thomas?

2019-09-23 Thread Phạm Ngọc Kiên
Dear Prof. Wolfgang Bangerth, The function fe.has_support_points() returns false. This means that the FE_Raviart_Thomas does not have support points. I think generalized_support_points are used for this type of finite element. For FE_RaviartThomas

Re: [deal.II] How to get the coordinates of dof with FE_Raviart_Thomas?

2019-09-23 Thread Wolfgang Bangerth
Pham, > With FE_Raviart_Thomas in 3D, the dofs are defined normal to the faces of > each > cell. > For example, if we use FE_RaviartThomas > (0), > > then in a cell we have 6 dofs d

[deal.II] How to get the coordinates of dof with FE_Raviart_Thomas?

2019-09-23 Thread Phạm Ngọc Kiên
Hi colleagues, With FE_Raviart_Thomas in 3D, the dofs are defined normal to the faces of each cell. For example, if we use FE_RaviartThomas (0), then in a cell we have 6 dofs defined o

Re: [deal.II] Create a LinearAlgebra::distributed::Vector from a LinearAlgebra::distributed::Vector

2019-09-23 Thread 'Maxi Miller' via deal.II User Group
Yes, that worked, thanks! Am Montag, 23. September 2019 17:17:09 UTC+2 schrieb Martin Kronbichler: > > Dear Maxi, > > Did you try to call 'operator=' in the same line as you declare the > variable? In that case, you still go through the (missing) constructor. It > should work if you first declar

Re: [deal.II] Create a LinearAlgebra::distributed::Vector from a LinearAlgebra::distributed::Vector

2019-09-23 Thread Martin Kronbichler
Dear Maxi, Did you try to call 'operator=' in the same line as you declare the variable? In that case, you still go through the (missing) constructor. It should work if you first declare the float vector in one statement and then assign in the next. We have this function: https://dealii.org/

Re: [deal.II] Create a LinearAlgebra::distributed::Vector from a LinearAlgebra::distributed::Vector

2019-09-23 Thread 'Maxi Miller' via deal.II User Group
When doing that in my code, I get the error error: conversion from ‘Vector’ to non-scalar type ‘Vector< float,[...]>’ requested thus I assume I have to implement it myself... Am Freitag, 13. September 2019 17:53:37 UTC+2 schrieb Wolfgang Bangerth: > > On 9/13/19 8:39 AM, 'Maxi Miller' via deal.II