is a general finite element solver using mixed elements to discretize de Rham-type complexes in one, two and three dimensions.
| Gauss_Pyramid.cpp (getdp-3.4.0-source.tgz) | : | Gauss_Pyramid.cpp (getdp-3.5.0-source.tgz) |
|---|
|
| // GetDP - Copyright (C) 1997-2021 P. Dular and C. Geuzaine, University of Liege | | // GetDP - Copyright (C) 1997-2022 P. Dular and C. Geuzaine, University of Liege |
| // | | // |
| // See the LICENSE.txt file for license information. Please report all | | // See the LICENSE.txt file for license information. Please report all |
| // issues on https://gitlab.onelab.info/getdp/getdp/issues. | | // issues on https://gitlab.onelab.info/getdp/getdp/issues. |
| | |
| #include "Gauss_Pyramid.h" | | #include "Gauss_Pyramid.h" |
| #include "Message.h" | | #include "Message.h" |
| | |
| /* Gauss integration over a pyramid */ | | /* Gauss integration over a pyramid */ |
| | |
|
| void Gauss_Pyramid(int Nbr_Points, int Num, | | void Gauss_Pyramid(int Nbr_Points, int Num, double *u, double *v, double *w, |
| double *u, double *v, double *w, double *wght) | | double *wght) |
| { | | { |
|
| switch (Nbr_Points) { | | switch(Nbr_Points) { |
| | | case 8: |
| case 8 : | | *u = upyr8[Num]; |
| *u = upyr8[Num]; *v = vpyr8[Num]; *w = wpyr8[Num]; *wght = ppyr8[Num]; | | *v = vpyr8[Num]; |
| break ; | | *w = wpyr8[Num]; |
| | *wght = ppyr8[Num]; |
| | break; |
| | |
|
| default : | | default: |
| Message::Error("Wrong number of Gauss points for Pyramid: " | | Message::Error("Wrong number of Gauss points for Pyramid: " |
| "valid choice: 8"); | | "valid choice: 8"); |
| break; | | break; |
| } | | } |
| } | | } |
| | | |
| End of changes. 4 change blocks. |
|---|
| 9 lines changed or deleted | | 11 lines changed or added |
|---|
"Fossies" - the Fresh Open Source Software Archive 