Skip to content
Snippets Groups Projects
Commit 41a1b26d authored by Christophe Geuzaine's avatar Christophe Geuzaine
Browse files

corrige bug prism/hexa dans unv

parent cbaac534
Branches
Tags
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 11 additions and 9 deletions
Mesh/Print_Mesh.cpp
+ 11
− 9
View file @ 41a1b26d
/* $Id: Print_Mesh.cpp,v 1.3 2000-11-23 23:20:35 geuzaine Exp $ */ /* $Id: Print_Mesh.cpp,v 1.4 2000-11-24 00:45:53 geuzaine Exp $ */
#include "Gmsh.h" #include "Gmsh.h"
#include "Const.h" #include "Const.h"
...@@ -501,6 +501,8 @@ int process_3D_elements (FILE * funv, Mesh * m){ ...@@ -501,6 +501,8 @@ int process_3D_elements (FILE * funv, Mesh * m){
List_T *ListVolumes = Tree2List (m->Volumes); List_T *ListVolumes = Tree2List (m->Volumes);
List_T *Elements; List_T *Elements;
Simplex *sx; Simplex *sx;
Hexahedron *hx;
Prism *px;
Volume *v; Volume *v;
int nb = 0, i, j, nsup, n, ntot, k, geo, fetyp; int nb = 0, i, j, nsup, n, ntot, k, geo, fetyp;
...@@ -553,8 +555,8 @@ int process_3D_elements (FILE * funv, Mesh * m){ ...@@ -553,8 +555,8 @@ int process_3D_elements (FILE * funv, Mesh * m){
// PRISMS // PRISMS
Elements = Tree2List (v->Prisms); Elements = Tree2List (v->Prisms);
for (j = 0; j < List_Nbr (Elements); j++){ for (j = 0; j < List_Nbr (Elements); j++){
List_Read (Elements, j, &sx); List_Read (Elements, j, &px);
if (sx->VSUP){ if (px->VSUP){
fetyp = WEDGE; fetyp = WEDGE;
n = 6; n = 6;
nsup = 9; nsup = 9;
...@@ -570,13 +572,13 @@ int process_3D_elements (FILE * funv, Mesh * m){ ...@@ -570,13 +572,13 @@ int process_3D_elements (FILE * funv, Mesh * m){
ELEMENT_ID++, fetyp, geo, geo, 7, n + nsup); ELEMENT_ID++, fetyp, geo, geo, 7, n + nsup);
ntot = 0; ntot = 0;
for (k = 0; k < n; k++){ for (k = 0; k < n; k++){
fprintf (funv, "%10d", sx->V[k]->Num); fprintf (funv, "%10d", px->V[k]->Num);
if (ntot % 8 == 7) if (ntot % 8 == 7)
fprintf (funv, "\n"); fprintf (funv, "\n");
ntot++; ntot++;
} }
for (k = 0; k < nsup; k++){ for (k = 0; k < nsup; k++){
fprintf (funv, "%10d", sx->VSUP[k]->Num); fprintf (funv, "%10d", px->VSUP[k]->Num);
if (ntot % 8 == 7) if (ntot % 8 == 7)
fprintf (funv, "\n"); fprintf (funv, "\n");
ntot++; ntot++;
...@@ -590,8 +592,8 @@ int process_3D_elements (FILE * funv, Mesh * m){ ...@@ -590,8 +592,8 @@ int process_3D_elements (FILE * funv, Mesh * m){
// HEXAHEDRONS // HEXAHEDRONS
Elements = Tree2List (v->Hexahedra); Elements = Tree2List (v->Hexahedra);
for (j = 0; j < List_Nbr (Elements); j++){ for (j = 0; j < List_Nbr (Elements); j++){
List_Read (Elements, j, &sx); List_Read (Elements, j, &hx);
if (sx->VSUP){ if (hx->VSUP){
fetyp = BRICK; fetyp = BRICK;
n = 8; n = 8;
nsup = 12; nsup = 12;
...@@ -607,13 +609,13 @@ int process_3D_elements (FILE * funv, Mesh * m){ ...@@ -607,13 +609,13 @@ int process_3D_elements (FILE * funv, Mesh * m){
ELEMENT_ID++, fetyp, geo, geo, 7, n + nsup); ELEMENT_ID++, fetyp, geo, geo, 7, n + nsup);
ntot = 0; ntot = 0;
for (k = 0; k < n; k++){ for (k = 0; k < n; k++){
fprintf (funv, "%10d", sx->V[k]->Num); fprintf (funv, "%10d", hx->V[k]->Num);
if (ntot % 8 == 7) if (ntot % 8 == 7)
fprintf (funv, "\n"); fprintf (funv, "\n");
ntot++; ntot++;
} }
for (k = 0; k < nsup; k++){ for (k = 0; k < nsup; k++){
fprintf (funv, "%10d", sx->VSUP[k]->Num); fprintf (funv, "%10d", hx->VSUP[k]->Num);
if (ntot % 8 == 7) if (ntot % 8 == 7)
fprintf (funv, "\n"); fprintf (funv, "\n");
ntot++; ntot++;
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment