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

fix

parent faf9c317
No related branches found
No related tags found
No related merge requests found
...@@ -20,6 +20,10 @@ void outputScalarField(std::list<BDS_Face*> t, const char *iii, int param, GFace ...@@ -20,6 +20,10 @@ void outputScalarField(std::list<BDS_Face*> t, const char *iii, int param, GFace
if (gf){ if (gf){
FILE* view_c = Fopen("param_c.pos","w"); FILE* view_c = Fopen("param_c.pos","w");
if(!view_c){
Msg::Error("Could not open file 'param_c.pos");
return;
}
fprintf(view_c,"View \"paramC\"{\n"); fprintf(view_c,"View \"paramC\"{\n");
std::set<MEdge,Less_Edge> all; std::set<MEdge,Less_Edge> all;
std::list<BDS_Face*>::iterator tit = t.begin(); std::list<BDS_Face*>::iterator tit = t.begin();
...@@ -37,7 +41,7 @@ void outputScalarField(std::list<BDS_Face*> t, const char *iii, int param, GFace ...@@ -37,7 +41,7 @@ void outputScalarField(std::list<BDS_Face*> t, const char *iii, int param, GFace
SPoint2 p = p1*(1.-t) + p2*t; SPoint2 p = p1*(1.-t) + p2*t;
GPoint pa = gf->point(p.x(),p.y()); GPoint pa = gf->point(p.x(),p.y());
GPoint pb = gf->point(prev.x(),prev.y()); GPoint pb = gf->point(prev.x(),prev.y());
fprintf(view_c,"SL(%g,%g,%g,%g,%g,%g){1,1,1};\n",pa.x(),pa.y(),pa.z(),pb.x(),pb.y(),pb.z()); fprintf(view_c,"SL(%g,%g,%g,%g,%g,%g){1,1,1};\n",pa.x(),pa.y(),pa.z(),pb.x(),pb.y(),pb.z());
prev = p; prev = p;
} }
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment