Arcane  v3.14.10.0
Documentation développeur
Chargement...
Recherche...
Aucune correspondance
HyodaIceT.cc
1// -*- tab-width: 2; indent-tabs-mode: nil; coding: utf-8-with-signature -*-
2//-----------------------------------------------------------------------------
3// Copyright 2000-2022 CEA (www.cea.fr) IFPEN (www.ifpenergiesnouvelles.com)
4// See the top-level COPYRIGHT file for details.
5// SPDX-License-Identifier: Apache-2.0
6//-----------------------------------------------------------------------------
7/*****************************************************************************
8 * HyodaIceT.cc (C) 2000-2012 *
9 *****************************************************************************/
10
11#include "arcane/IMesh.h"
12#include "arcane/IApplication.h"
13#include "arcane/IParallelMng.h"
14#include "arcane/IVariableMng.h"
15#include "arcane/utils/ITraceMng.h"
16#include "arcane/IVariableAccessor.h"
17#include "arcane/FactoryService.h"
18#include "arcane/ServiceFinder2.h"
19#include "arcane/SharedVariable.h"
20#include "arcane/CommonVariables.h"
21#include "arcane/utils/ScopedPtr.h"
22#include "arcane/AbstractService.h"
23#include "arcane/VariableCollection.h"
24#include "arcane/utils/ArcaneGlobal.h"
25#include "arcane/utils/PlatformUtils.h"
26#include "arcane/utils/IOnlineDebuggerService.h"
27#include "arcane/ITransferValuesParallelOperation.h"
28
29#include "arcane/IVariableAccessor.h"
30#include <arcane/IParticleFamily.h>
31#include "arcane/datatype/ArrayVariant.h"
32
33#ifndef GL_MULTISAMPLE
34#define GL_MULTISAMPLE 0x809D
35#endif
36
37#include "arcane/hyoda/HyodaIceTGlobal.h"
38#include "arcane/hyoda/Hyoda.h"
39#include "arcane/hyoda/HyodaArc.h"
40#include "arcane/hyoda/HyodaTcp.h"
41#include "arcane/hyoda/HyodaIceT.h"
42#include "arcane/hyoda/HyodaMix.h"
43extern void ppmWrite(const char*, const unsigned char*, int width, int height);
44
45namespace
46{
47 static Arcane::HyodaIceT* HyodaIceTCallbackHandle = NULL;
48 void drawGLCallback(void){
49 if (!HyodaIceTCallbackHandle) return;
50 HyodaIceTCallbackHandle->drawGL();
51 }
52}
53
54
55/*---------------------------------------------------------------------------*/
56/*---------------------------------------------------------------------------*/
57
58ARCANE_BEGIN_NAMESPACE
59
60/*---------------------------------------------------------------------------*/
61/*---------------------------------------------------------------------------*/
62
63
64/******************************************************************************
65 * HyodaIceT
66 *****************************************************************************/
67HyodaIceT::HyodaIceT(Hyoda *hd,
68 ISubDomain *sd,
69 ITraceMng *tm,
70 unsigned int w,
71 unsigned int h,
72 HyodaTcp *_tcp): TraceAccessor(tm),
73 hyoda(hd),
74 m_sub_domain(sd),
75 icetContext(NULL),
76 icetCommunicator(NULL),
77 osMesaContext(NULL),
78 rank(sd->parallelMng()->commRank()),
79 num_proc(sd->parallelMng()->commSize()),
80 m_screen_width(w),
81 m_screen_height(h),
82 m_image_buffer(NULL),
83 scale(1.),
84 rot_x(0.1),
85 rot_y(0.),
86 rot_z(0.),
87 variable_index(0),
88 variable_plugin(0),
89 m_tcp(_tcp),
90 m_hyoda_mix(new HyodaMix(hd,sd,tm)),
91 m_variable(NULL)
92{
93 debug()<<"\33[33m[HyodaIceT::HyodaIceT] New "
94 << ", screen_width="<<m_screen_width
95 << ", screen_height="<<m_screen_height
96 << "\33[m";
97
98 m_pov_sxyzip[0]=
99 m_pov_sxyzip[1]=
100 m_pov_sxyzip[2]=
101 m_pov_sxyzip[3]=
102 m_pov_sxyzip[4]=
103 m_pov_sxyzip[5]=0.0;
104
105 #warning Duplicated code here
106 // On cherche la premier IVariable à afficher
107 VariableCollection variables = m_sub_domain->variableMng()->usedVariables();
108 for(VariableCollection::Enumerator ivar(variables); ++ivar; ){
109 IVariable* var = *ivar;
110 // Pas de références, pas de variable
111 if (var->nbReference()==0) {debug() << "\33[33m[HyodaIceT::HyodaIceT] No nbReference\33[m"; continue;}
112 // Pas sur le bon support, pas de variable
113 if (var->itemKind() != IK_Node &&
114 var->itemKind() != IK_Cell &&
115 var->itemKind() != IK_Face &&
116 var->itemKind() != IK_Particle) continue;
117 // Pas réclamée en tant que PostProcess'able, pas de variable
118 if (var->itemKind() != IK_Particle && (!var->hasTag("PostProcessing")))
119 {debug() << "\33[33m[HyodaIceT::HyodaIceT] "<<var->name()<<" not PostProcessing\33[m"; continue;}
120 // Pas de type non supportés
121 //if (var->dataType()>=DT_String) continue;
122 if (var->dataType()!=DT_Real) continue;
123 m_variable=var;
124 debug() << "\33[33m[HyodaIceT::HyodaIceT] m_variable is "<<var->name()<<"\33[m";
125 break;
126 }
127 if (!m_variable)
128 ARCANE_FATAL("\33[33m[HyodaIceT::HyodaIceT] No m_variable to draw!\33[m");
129 debug()<<"\33[33m[HyodaIceT::HyodaIceT] Focusing on variable: "<<m_variable->name()<<"\33[m";
130
131 /* specify Z, stencil, accum sizes */
132 osMesaContext = OSMesaCreateContextExt(OSMESA_RGBA, 8, 0, 0, NULL);
133 if (!osMesaContext) fatal()<<"OSMesaCreateContext failed!\n";
134 /* Allocate the image buffer */
135 m_image_buffer=malloc(renderSize());
136 if (!m_image_buffer) fatal()<<"Alloc image failed!\n";
137 /* Bind the image to the context and make it current */
138 if (!OSMesaMakeCurrent(osMesaContext,
139 m_image_buffer,
140 GL_UNSIGNED_BYTE,
141 m_screen_width,
142 m_screen_height ))
143 fatal()<<"OSMesaMakeCurrent failed!\n";
144 {
145 int z, s, a;
146 glGetIntegerv(GL_DEPTH_BITS, &z);
147 glGetIntegerv(GL_STENCIL_BITS, &s);
148 glGetIntegerv(GL_ACCUM_RED_BITS, &a);
149 debug()<<"\33[33m[HyodaIceT] Depth="<<z<<", Stencil="<<s<<", Accum="<<a<<"\33[m";
150 }
151 // IceT system init
152 #warning IceT needs MPI_Init
153 IceTBitField diag_level = ICET_DIAG_ALL_NODES | ICET_DIAG_WARNINGS;
154 if (sd->parallelMng()->isParallel()){
155 debug() << "\33[33m[HyodaIceT] isParallel icetCreateMPICommunicator @"
156 << sd->parallelMng()->getMPICommunicator()<< "\33[m";
157 icetCommunicator = icetCreateMPICommunicator(*(MPI_Comm*)sd->parallelMng()->getMPICommunicator());
158 }else{
159 debug()<<"\33[33m[HyodaIceT] MPI_COMM_SELF icetCreateMPICommunicator\33[m";
160 icetCommunicator = icetCreateMPICommunicator(MPI_COMM_SELF);
161 }
162 icetContext = icetCreateContext(icetCommunicator);
163
164 // iceT GL
165 icetGLInitialize();
166 icetGLDrawCallback(drawGLCallback);
167 HyodaIceTCallbackHandle=this;
168 icetResetTiles();
169 icetAddTile(0, 0, m_screen_width, m_screen_height, 0);
170 icetStrategy(ICET_STRATEGY_SEQUENTIAL);
171 icetSingleImageStrategy(ICET_SINGLE_IMAGE_STRATEGY_AUTOMATIC);
172 //icetDiagnostics(diag_level);
173 icetSetColorFormat(ICET_IMAGE_COLOR_RGBA_UBYTE);
174 icetSetDepthFormat(ICET_IMAGE_DEPTH_FLOAT);
175 icetCompositeMode(ICET_COMPOSITE_MODE_Z_BUFFER);
176 //icetBoundingBoxd(-1.0f, +1.0f, -1.0f, +1.0f, -4.0f, 16.0f);
177 setLeftRightBottomTop();
178 Real length=lrbtnf[1]-lrbtnf[0];
179 Real height=lrbtnf[3]-lrbtnf[2];
180 icetBoundingBoxd(lrbtnf[0]-length, lrbtnf[1]+length, lrbtnf[2]-height, lrbtnf[3]+height, lrbtnf[4], lrbtnf[5]);
181 icetDisable(ICET_FLOATING_VIEWPORT);
182 icetEnable(ICET_ORDERED_COMPOSITE);
183 initGL();
184 checkIceTError();
185 debug()<<"\33[33m[HyodaIceT] checkOglError\33[m";
186 checkOglError();
187 traceMng()->flush();
188 }
189
190
191
192/******************************************************************************
193 * Initialisation de la base ortho qui servira à glOrtho: left, right, bottom, top
194 *****************************************************************************/
195void HyodaIceT::setLeftRightBottomTop(void){
196 Real scale=1.0;
197 const VariableItemReal3 nodes_coords=m_sub_domain->defaultMesh()->nodesCoordinates();
198 debug()<<"[HyodaIceT::setLeftRightBottomTop] mesh()->dimension()="<<m_sub_domain->mesh()->dimension();
199 lrbtnf[0]=lrbtnf[1]=lrbtnf[2]=lrbtnf[3]=0.0;
200 lrbtnf[4]=-4.0; // On commence avec un Z-min par défaut pour le 2D
201 lrbtnf[5]=+16.0; // Idem avec le Z-max, pour le 2D par défaut
202 //if (m_sub_domain->mesh()->dimension()==3) scale=1.0;
203 ENUMERATE_NODE(node,m_sub_domain->defaultMesh()->allCells().outerFaceGroup().nodeGroup()){
204 //ENUMERATE_NODE(node, m_sub_domain->defaultMesh()->ownNodes()){
205 lrbtnf[0]=math::min(lrbtnf[0],nodes_coords[node].x); // left
206 lrbtnf[1]=math::max(lrbtnf[1],nodes_coords[node].x); // right
207 lrbtnf[2]=math::min(lrbtnf[2],nodes_coords[node].y); // bottom
208 lrbtnf[3]=math::max(lrbtnf[3],nodes_coords[node].y); // top
209 lrbtnf[4]=math::min(lrbtnf[4],nodes_coords[node].z); // near
210 lrbtnf[5]=math::max(lrbtnf[5],nodes_coords[node].z); // far
211 }
212 UniqueArray<Real> rdcMin(3);
213 UniqueArray<Real> rdcMax(3);
214 rdcMin[0]=scale*lrbtnf[0]; // left
215 rdcMin[1]=scale*lrbtnf[2]; // bottom
216 rdcMin[2]=scale*lrbtnf[4]; // near
217 m_sub_domain->parallelMng()->reduce(Parallel::ReduceMin,rdcMin.view());
218 rdcMax[0]=scale*lrbtnf[1]; // right
219 rdcMax[1]=scale*lrbtnf[3]; // top
220 rdcMax[2]=scale*lrbtnf[5]; // far
221 m_sub_domain->parallelMng()->reduce(Parallel::ReduceMax,rdcMax.view());
222 lrbtnf[0]=rdcMin[0]; // left
223 lrbtnf[1]=rdcMax[0]; // right
224 lrbtnf[2]=rdcMin[1]; // bottom
225 lrbtnf[3]=rdcMax[1]; // top
226 lrbtnf[4]=rdcMin[2]; // near
227 lrbtnf[5]=rdcMax[2]; // far
228 debug()<<"\33[33m[HyodaIceT::setLeftRightBottomTop] left =" << lrbtnf[0]<<"\33[m";
229 debug()<<"\33[33m[HyodaIceT::setLeftRightBottomTop] right =" << lrbtnf[1]<<"\33[m";
230 debug()<<"\33[33m[HyodaIceT::setLeftRightBottomTop] bottom =" << lrbtnf[2]<<"\33[m";
231 debug()<<"\33[33m[HyodaIceT::setLeftRightBottomTop] top =" << lrbtnf[3]<<"\33[m";
232 debug()<<"\33[33m[HyodaIceT::setLeftRightBottomTop] near =" << lrbtnf[4]<<"\33[m";
233 debug()<<"\33[33m[HyodaIceT::setLeftRightBottomTop] far =" << lrbtnf[5]<<"\33[m";
234}
235
236
237/******************************************************************************
238 * sxyz
239 *****************************************************************************/
240void HyodaIceT::sxyzip(double *pov){
241 scale=pov[0];
242 rot_x=pov[1];
243 rot_y=pov[2];
244 rot_z=pov[3];
245 debug()<<"\33[33m[HyodaIceT::sxyzi] scale="<<scale<<", rot_x="<<rot_x<<", rot_y="<<rot_y<<", rot_z="<<rot_z<<"\33[m";
246 variable_index=(int)pov[4];
247 variable_plugin=(int)pov[5];
248 if (variable_index<0){
249 debug()<<"\33[33m[HyodaIceT::sxyzi] Fetched negative variable_index!\33[m";
250 variable_index=0;
251 }
252 if (variable_plugin<0 || variable_plugin>2){
253 debug()<<"\33[33m[HyodaIceT::sxyzi] Fetched uncompatible variable_plugin!\33[m";
254 variable_plugin=0;
255 }
256
257 // Scrutation de la variable ARCANE ciblée pour être affichée
258 int i=0;
259 VariableCollection variables = m_sub_domain->variableMng()->usedVariables();
260 debug()<<"\33[33m[HyodaIceT::sxyzi] Looking for variable #"<<variable_index<<"\33[m";
261 for(VariableCollection::Enumerator ivar(variables); ++ivar;){
262 m_variable = *ivar;
263 if (m_variable->nbReference()==0) continue;
264 if (m_variable->itemKind()!=IK_Node &&
265 m_variable->itemKind()!=IK_Cell &&
266 m_variable->itemKind()!=IK_Face &&
267 m_variable->itemKind()!=IK_Particle) continue;
268 if (m_variable->itemKind()!=IK_Particle && (!m_variable->hasTag("PostProcessing"))) continue;
269 //if (variable->dataType()>=DT_String) continue;
270 if (m_variable->dataType()!=DT_Real) continue;
271 if (i==variable_index) break;
272 if (i>variable_index)
273 fatal()<<"\33[33m[HyodaIceT::sxyzi] Could not find variable to draw (variable_index="<<variable_index<<", i="<<i<<")!\33[m\n";
274 i+=1;
275 }
276 debug()<<"\33[33m[HyodaIceT::sxyzi] " << m_variable->name()<<"\33[m";
277}
278
279
280/******************************************************************************
281 * imgSize
282 *****************************************************************************/
283int HyodaIceT::renderSize(void){
284 return m_screen_width*m_screen_height*4*sizeof(GLubyte);
285}
286
287
288/******************************************************************************
289 * ~HyodaIceT
290 *****************************************************************************/
291HyodaIceT::~HyodaIceT(){
292 free(m_image_buffer);
293 checkOglError();
294 checkIceTError();
295 debug()<<"\33[33m~HyodaIceT\33[m";
296 icetDestroyContext(icetContext);
297 OSMesaDestroyContext(osMesaContext);
298}
299
300
301
302/******************************************************************************
303 * render
304 *****************************************************************************/
305void HyodaIceT::render(void){
306 debug()<<"\33[33m[HyodaIceT] render\33[m";
307 icetGLDrawCallback(drawGLCallback);
308 // Tout le monde participe au rendering
309 // Mais on blinde cette section contre les erreurs FPU qu'iceT a tendance à jeter
310 //platform::enableFloatingException(false);
311 platform::enableFloatingException(false);
312 IceTImage image=icetGLDrawFrame();
313 platform::enableFloatingException(true);
314 checkOglError();
315 checkIceTError();
316
317 m_sub_domain->parallelMng()->barrier();
318
319 // Seul le CPU 0 possède l'image résultante pour l'instant
320 if (rank==0){
321 // Now fetch image's address
322 unsigned char* imgAddress=icetImageGetColorub(image);
323
324 // Si la variable d'environement ARCANE_HYODA_MESH_PPM est settée, on dump l'image dans un fichier
325 if (!platform::getEnvironmentVariable("ARCANE_HYODA_MESH_PPM").null()){
326 char filename[3883];
327 if (29!=snprintf(filename,3883,"/tmp/HyODA_Msh_%04d.ppm", m_sub_domain->commonVariables().globalIteration()))
328 return;
329 debug()<<"\33[33m[HyodaIceT::render] \33[33mARCANE_HYODA_MESH_PPM\33[m "
330 << icetImageGetWidth(image) << "x" << icetImageGetHeight(image)
331 << " image " << filename;
332 ppmWrite(filename, imgAddress,(int)m_screen_width, (int)m_screen_height);
333 }
334
335 debug()<<"\33[33m[HyodaIceT::render] sending header packet for image!\33[m";
336 char icePacketHeader[8];
337 // On écrit le 'MeshIceTHeader' QHyodaTcpSwitch
338 *(unsigned int*)&icePacketHeader[0]=0xcbce69bcul;
339 // On pousse la taille du paquet
340 *(unsigned int*)&icePacketHeader[4]=8;
341 m_tcp->send(icePacketHeader,8);
342 m_tcp->waitForAcknowledgment();
343 debug()<<"\33[33m[HyodaIceT::render] sending image packet!, renderSize="<<renderSize()<<"\33[m";
344 m_tcp->send((void*)imgAddress, renderSize());
345 m_tcp->waitForAcknowledgment();
346 debug()<<"\33[33m[HyodaIceT::render] Waiting for POV answer!\33[m";
347 m_tcp->recvPov(m_pov_sxyzip);
348 }
349
350 m_sub_domain->parallelMng()->barrier();
351
352 // On dépiote la réponse
353 UniqueArray<Real> pov_to_broadcast(0);
354 pov_to_broadcast.add(m_pov_sxyzip[0]);
355 pov_to_broadcast.add(m_pov_sxyzip[1]);
356 pov_to_broadcast.add(m_pov_sxyzip[2]);
357 pov_to_broadcast.add(m_pov_sxyzip[3]);
358 pov_to_broadcast.add(m_pov_sxyzip[4]);
359 pov_to_broadcast.add(m_pov_sxyzip[5]);
360 debug()<<"\33[33m[Hyoda::hook] broadcasting pov...\33[m";
361 m_sub_domain->parallelMng()->broadcast(pov_to_broadcast.view(),0);
362 m_pov_sxyzip[0]=pov_to_broadcast[0];
363 m_pov_sxyzip[1]=pov_to_broadcast[1];
364 m_pov_sxyzip[2]=pov_to_broadcast[2];
365 m_pov_sxyzip[3]=pov_to_broadcast[3];
366 m_pov_sxyzip[4]=pov_to_broadcast[4];
367 m_pov_sxyzip[5]=pov_to_broadcast[5];
368 sxyzip(m_pov_sxyzip);
369 debug()<<"\33[33m[Hyoda::hook] pov done\33[m";
370 debug()<<"\33[33m[HyodaIceT::render] sxyzip:"
371 << "scale=" << m_pov_sxyzip[0]
372 << ", x=" << m_pov_sxyzip[1]
373 << ", y=" << m_pov_sxyzip[2]
374 << ", z=" << m_pov_sxyzip[3]
375 << ", i=" << m_pov_sxyzip[4]
376 << ", p=" << m_pov_sxyzip[5]
377 << "\33[m";
378}
379
380
381/******************************************************************************
382 * GL init
383 *****************************************************************************/
384void HyodaIceT::initGL(void){
385 debug()<<"\33[33m[HyodaIceT] initGL in\33[m";
386 checkOglError();
387 glViewport(0, 0, m_screen_width,m_screen_height);
388 glEnable(GL_DITHER);
389 glEnable(GL_DEPTH_TEST);
390 glEnable(GL_MULTISAMPLE);
391 glEnable(GL_POINT_SMOOTH);
392 //glEnable(GL_POLYGON_SMOOTH);
393 //glEnable(GL_LINE_SMOOTH);
394 glPointSize(1.0);
395 glLineWidth(1.0);
396 glDepthFunc(GL_LESS);
397 glShadeModel(GL_SMOOTH);
398 glShadeModel(GL_FLAT);
399 glLineStipple(1, 0x0101);
400 glPolygonMode(GL_FRONT_AND_BACK, GL_FILL);
401 glClearDepth(1.0);
402 glColor3f(0.0f, 0.0f, 0.0f); // Black points
403 glClearColor(0.0, 0.0, 0.0, 0.0); // Black background
404 glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
405 debug()<<"\33[33m[HyodaIceT] initGL out\33[m";
406 checkOglError();
407}
408
409
410/******************************************************************************
411 * drawGLCallback
412 *****************************************************************************/
413void HyodaIceT::drawGL(void){
414 debug()<<"\33[33m[HyodaIceT] drawGL\33[m";
415 glPushMatrix();
416 glShadeModel(GL_SMOOTH);
417 glMatrixMode(GL_PROJECTION);
418 glLoadIdentity();
419 Real length=lrbtnf[1]-lrbtnf[0];
420 Real height=lrbtnf[3]-lrbtnf[2];
421 //glOrtho(lrbtnf[0], lrbtnf[1], lrbtnf[2], lrbtnf[3], lrbtnf[4], lrbtnf[5]);
422 glOrtho(lrbtnf[0]-length/4., lrbtnf[1]+length/4., lrbtnf[2]-height/4., lrbtnf[3]+height/4., lrbtnf[4], lrbtnf[5]);
423 //glFrustum(lrbtnf[0], lrbtnf[1], lrbtnf[2], lrbtnf[3], lrbtnf[4], lrbtnf[5]);
424 //glTranslatef(0.0f, 0.0f, -4.0f);
425 // et on applique les transformation venant de l'IHM
426 glRotatef(rot_x, 1.0f, 0.0f, 0.0f);
427 glRotatef(rot_y, 0.0f, 1.0f, 0.0f);
428 glRotatef(rot_z, 0.0f, 0.0f, 1.0f);
429 glScalef(scale, scale, scale);
430 glMatrixMode(GL_MODELVIEW);
431 glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
432 drawArcaneMesh();
433 glPopMatrix();
434 glFlush();
435 glFinish();
436 checkOglError();
437}
438
439
440/******************************************************************************
441 * drawArcaneMesh
442 *****************************************************************************/
443void HyodaIceT::drawArcaneMesh(void){
444 //const register VariableItemReal3 nodes_coords=m_sub_domain->defaultMesh()->nodesCoordinates();
445 // S'il n'y a pas de variable à afficher, il n'y a rien à scruter
446 /*if (m_tcp->getVariablesNumber()==0 || m_tcp->getVariablesNumber()==-1) {
447 debug()<<"[HyodaIceT::drawArcaneMesh] No variables have been set as PostProcessing!";
448 drawArcPolygonsWithoutVariable();
449 return;
450 }*/
451 debug()<<"\33[33m[HyodaIceT] drawArcaneMesh\33[m";
452 drawArcPolygons();
453 checkOglError();
454}
455
456
457// ****************************************************************************
458// * iceColorMinMax
459// ****************************************************************************
460void HyodaIceT::iceColorMinMax(Real &min, Real &max){
461 String item_family_name(m_variable->itemFamilyName()==NULL?
462 "particles":m_variable->itemFamilyName());
463 VariableItemReal var(VariableBuildInfo(m_sub_domain->defaultMesh(),
464 m_variable->name(),item_family_name),
465 m_variable->itemKind());
466 ENUMERATE_ITEM(item, m_variable->itemFamily()->allItems().own()){
467 Real val=var[item];
468 min = math::min(min,val);
469 max = math::max(max,val);
470 }
471 min=m_sub_domain->parallelMng()->reduce(Parallel::ReduceMin, min);
472 max=m_sub_domain->parallelMng()->reduce(Parallel::ReduceMax, max);
473 debug()<<"\33[33m[iceColorMinMax] min="<<min<<", max="<<max<<"\33[m";
474}
475
476
477// *****************************************************************************
478// * drawArcPolygons
479// *****************************************************************************
480void HyodaIceT::drawArcPolygons(void){
481 Real3 rgb;
482 Real min=+HUGE_VAL;
483 Real max=-HUGE_VAL;
484 String item_family_name(m_variable->itemFamilyName()==NULL?
485 "particles":m_variable->itemFamilyName());
486 VariableNodeReal3 nCoords = m_sub_domain->defaultMesh()->nodesCoordinates();
487 VariableItemReal var(VariableBuildInfo(m_sub_domain->defaultMesh(),
488 m_variable->name(),item_family_name),
489 m_variable->itemKind());
490
491 // On calcule le min max des variables globales
492 iceColorMinMax(min,max);
493 if (m_sub_domain->defaultMesh()->dimension()==2){
494 debug()<<"\33[33m[drawArcPolygons] 2D\33[m";
495 // On regarde si un plugins est réclamé depuis l'interface
496 if (variable_plugin!=0){
497 debug()<<"\33[33m[drawArcPolygons] variable_plugin="<<variable_plugin<<"\33[m";
498 m_hyoda_mix->xLine2Cell(variable_plugin, m_variable, min, max);
499 //debug()<<"[HyodaIceT] drawArcPolygons xLine2Cell";
500 checkOglError();
501 return;
502 }
503
504 debug()<<"\33[33m[drawArcPolygons] variable_plugin="<<variable_plugin<<"\33[m";
505 if ((m_variable->itemKind()==IK_Cell)||(m_variable->itemKind()==IK_Node)){
506 debug()<<"\33[33m[drawArcPolygons] m_variable->name()="<<m_variable->name()<<"\33[m";
507 debug()<<"\33[33m[drawArcPolygons] m_variable->itemKind()="<<m_variable->itemKind()<<"\33[m";
508 ENUMERATE_CELL(cell, m_sub_domain->defaultMesh()->ownCells()){
509 glBegin(GL_POLYGON);
510 if (m_variable->itemKind()==IK_Cell) setColor(min,max,var[cell],rgb);
511 ENUMERATE_NODE(node, cell->nodes()){
512 if (m_variable->itemKind()==IK_Node) setColor(min,max,var[node],rgb);
513 glColor3d(rgb[0], rgb[1], rgb[2]);
514 glVertex2d(nCoords[node].x, nCoords[node].y);
515 }
516 glEnd();
517 }
518 //debug()<<"[HyodaIceT] drawArcPolygons CELL NODE";
519 checkOglError();
520 return;
521 }
522
523 if (m_variable->itemKind()==IK_Face){
524 debug()<<"\33[33m[drawArcPolygons] IK_Face!\33[m";
525 ENUMERATE_FACE(face, m_sub_domain->defaultMesh()->ownFaces()){
526 glBegin(GL_LINES);
527 setColor(min,max,var[face],rgb);
528 glColor3d(rgb[0], rgb[1], rgb[2]);
529 ENUMERATE_NODE(node, face->nodes()){
530 glVertex2d(nCoords[node].x, nCoords[node].y);
531 }
532 glEnd();
533 }
534 //debug()<<"[HyodaIceT] drawArcPolygons FACE NODE";
535 checkOglError();
536 return;
537 }
538 } // dimension()==2
539
540 if (m_sub_domain->defaultMesh()->dimension()==3){
541 debug()<<"\33[33m[drawArcPolygons] 3D\33[m";
542 if (m_variable->itemKind()==IK_Cell){
543 ENUMERATE_FACE(face, m_sub_domain->defaultMesh()->outerFaces()){
544 if (!face->isOwn()) continue;
545 glBegin(GL_POLYGON);
546 setColor(min,max,var[face->cell(0)],rgb);
547 glColor3d(rgb[0], rgb[1], rgb[2]);
548 ENUMERATE_NODE(node, face->nodes()){
549 glVertex3d(nCoords[node].x,
550 nCoords[node].y,
551 nCoords[node].z);
552 }
553 glEnd();
554 }
555 }
556
557 if (m_variable->itemKind()==IK_Particle){
558 VariableParticleReal3 m_particle_r(VariableBuildInfo(m_sub_domain->defaultMesh(),"r","particles"));
559 ENUMERATE_PARTICLE(particle, m_sub_domain->defaultMesh()->findItemFamily("particles",true)->allItems()){
560 if (!particle->cell().isOwn()) continue;
561 glBegin(GL_POINTS);
562 setColor(min,max,var[particle],rgb);
563 glColor3d(rgb[0], rgb[1], rgb[2]);
564 glVertex3d(m_particle_r[particle].x,m_particle_r[particle].y,m_particle_r[particle].z);
565 glEnd();
566 }
567 }
568
569 // Permet de dessiner les lignes des arêtes
570 /*glBegin(GL_LINE_STRIP);
571 glColor3d(0.0, 0.0, 0.0);
572 ENUMERATE_NODE(node, m_sub_domain->defaultMesh()->outerFaces()){
573 glVertex3d(nCoords[node].x,
574 nCoords[node].y,
575 nCoords[node].z);
576 }
577 glEnd();*/
578 } // dimension()==3
579}
580
581
582// ****************************************************************************
583// * drawArcPolygons
584// ****************************************************************************
585void HyodaIceT::drawArcPolygonsWithoutVariable(void){
586 VariableNodeReal3 nCoords = m_sub_domain->defaultMesh()->nodesCoordinates();
587 ENUMERATE_CELL(cell, m_sub_domain->defaultMesh()->ownCells()){
588 glBegin(GL_POLYGON);
589 glColor3d(0.0, 0.0, 1.0);
590 ENUMERATE_NODE(node, cell->nodes()){
591 glVertex3d(nCoords[node].x,
592 nCoords[node].y,
593 nCoords[node].z);
594 }
595 glEnd();
596 }
597}
598
599
600/******************************************************************************
601 * drawArcPoints
602 *****************************************************************************/
603void HyodaIceT::drawArcPoints(const VariableItemReal3 &nodes_coords,
604 VariableItemReal &var,
605 double min,double max){
606 Real3 rgb;
607 // On trace les points
608 glBegin(GL_POINTS);
609 /*ENUMERATE_NODE(node, m_sub_domain->defaultMesh()->ownNodes()){
610 setColor(min,max,var[node],rgb);
611 glColor3d(rgb[0], rgb[1], rgb[2]);
612 glVertex3d(nodes_coords[node].x,
613 nodes_coords[node].y,
614 nodes_coords[node].z);
615 }*/
616 ENUMERATE_CELL(cell, m_sub_domain->defaultMesh()->ownCells()){
617 #warning Only cell Variables for drawArcPoints
618 setColor(min,max,var[cell],rgb);
619 ENUMERATE_NODE(node, cell->nodes()){
620 glColor3d(rgb[0], rgb[1], rgb[2]);
621 glVertex3d(nodes_coords[node].x,
622 nodes_coords[node].y,
623 nodes_coords[node].z);
624 }
625 }
626 glEnd();
627}
628
629
630/******************************************************************************
631 * drawArcLines
632 *****************************************************************************/
633void HyodaIceT::drawArcLines(const VariableItemReal3 &nodes_coords,
634 VariableItemReal&){
635// On trace les lignes
636 glBegin(GL_LINES);
637 //glColor3f(1.0, 1.0, 1.0);
638 glColor3f(0.0f, 0.0f, 0.0f);
639 ENUMERATE_FACE(face, m_sub_domain->defaultMesh()->allFaces()){
640 for(Integer i=0,mx=face->nodes().size()-1;i<mx;++i){
641 glVertex3d(nodes_coords[face->node(i)].x,
642 nodes_coords[face->node(i)].y,
643 nodes_coords[face->node(i)].z);
644 glVertex3d(nodes_coords[face->node(i+1)].x,
645 nodes_coords[face->node(i+1)].y,
646 nodes_coords[face->node(i+1)].z);
647 }
648 glVertex3d(nodes_coords[face->node(0)].x,
649 nodes_coords[face->node(0)].y,
650 nodes_coords[face->node(0)].z);
651 glVertex3d(nodes_coords[face->node(face->nodes().size()-1)].x,
652 nodes_coords[face->node(face->nodes().size()-1)].y,
653 nodes_coords[face->node(face->nodes().size()-1)].z);
654 }
655 glEnd();
656}
657
658
659/******************************************************************************
660 * color
661 *****************************************************************************/
662void HyodaIceT::
663setColor(double min, double max, double v, Real3 &rgb)
664{
665 ARCANE_ASSERT(min<=max,("setColor min<=max"));
666 if (min==max){
667 min=0.0;
668 max=1.0;
669 }
670 double mid=max-min;
671 ARCANE_ASSERT(mid!=0.,("setColor mid!=0."));
672 rgb.x=rgb.y=rgb.z=1.0;
673 if (v<min) v=min;
674 if (v>max) v=max;
675 if (v<(min+0.25*mid)){
676 rgb.x=0.;
677 rgb.y=4.*(v-min)/mid;
678 }else if (v<(min+0.5*mid)){
679 rgb.x=0.;
680 rgb.z=1.+4*(min+0.25*mid-v)/mid;
681 }else if (v<(min+0.75*mid)){
682 rgb.x=4*(v-min-0.5*mid)/mid;
683 rgb.z=0.;
684 }else{
685 rgb.y=1.+4*(min+0.75*mid-v)/mid;
686 rgb.z=0.;
687 }
688}
689
690
691/******************************************************************************
692 * checkIceTError
693 *****************************************************************************/
694void HyodaIceT::checkIceTError(void){
695#define CASE_ERROR(ename) case ename: debug() << "\33[33m## IceT status = " << #ename << "\33[m"; break;
696 switch (icetGetError()) {
697 case(ICET_NO_ERROR): break;
698 CASE_ERROR(ICET_SANITY_CHECK_FAIL);
699 CASE_ERROR(ICET_INVALID_ENUM);
700 CASE_ERROR(ICET_BAD_CAST);
701 CASE_ERROR(ICET_OUT_OF_MEMORY);
702 CASE_ERROR(ICET_INVALID_OPERATION);
703 CASE_ERROR(ICET_INVALID_VALUE);
704 default:debug()<<"\33[33m## UNKNOWN ICET ERROR CODE!!!!!\33[m";
705 }
706#undef CASE_ERROR
707}
708
709
710/******************************************************************************
711 * checkIceTError
712 *****************************************************************************/
713void HyodaIceT::checkOglError(void){
714#define CASE_ERROR(ename) case ename: debug() << "\33[33m## OpenGL status = " << #ename << "\33[m"; break;
715 switch (glGetError()) {
716 case(GL_NO_ERROR):break;
717 CASE_ERROR(GL_INVALID_ENUM);
718 CASE_ERROR(GL_INVALID_VALUE);
719 CASE_ERROR(GL_INVALID_OPERATION);
720 CASE_ERROR(GL_STACK_OVERFLOW);
721 CASE_ERROR(GL_STACK_UNDERFLOW);
722 CASE_ERROR(GL_OUT_OF_MEMORY);
723#ifdef GL_TABLE_TOO_LARGE
724 CASE_ERROR(GL_TABLE_TOO_LARGE);
725#endif
726 default:debug()<<"\33[33m## UNKNOWN GL ERROR CODE!!!!!\33[m";
727 }
728#undef CASE_ERROR
729}
730
731
732/*---------------------------------------------------------------------------*/
733/*---------------------------------------------------------------------------*/
734
735ARCANE_END_NAMESPACE
736
737/*---------------------------------------------------------------------------*/
738/*---------------------------------------------------------------------------*/
739
ArcaneGlobal.h
Fichier de configuration d'Arcane.
ARCANE_FATAL
#define ARCANE_FATAL(...)
Macro envoyant une exception FatalErrorException.
Definition ArcaneGlobal.h:809
ENUMERATE_FACE
#define ENUMERATE_FACE(name, group)
Enumérateur générique d'un groupe de faces.
Definition ItemEnumerator.h:424
ENUMERATE_PARTICLE
#define ENUMERATE_PARTICLE(name, group)
Enumérateur générique d'un groupe de particules.
Definition ItemEnumerator.h:430
ENUMERATE_CELL
#define ENUMERATE_CELL(name, group)
Enumérateur générique d'un groupe de mailles.
Definition ItemEnumerator.h:427
ENUMERATE_ITEM
#define ENUMERATE_ITEM(name, group)
Enumérateur générique d'un groupe de noeuds.
Definition ItemEnumerator.h:413
ENUMERATE_NODE
#define ENUMERATE_NODE(name, group)
Enumérateur générique d'un groupe de noeuds.
Definition ItemEnumerator.h:418
Arcane::LimaWrapper
Lecteur des fichiers de maillage via la bibliothèque LIMA.
Definition Lima.cc:120
Arcane::math::max
T max(const T &a, const T &b, const T &c)
Retourne le maximum de trois éléments.
Definition MathUtils.h:392
Arcane::math::min
ARCCORE_HOST_DEVICE Real2 min(Real2 a, Real2 b)
Retourne le minimum de deux Real2.
Definition MathUtils.h:336
Arcane::VariableParticleReal3
MeshVariableScalarRefT< Particle, Real3 > VariableParticleReal3
Grandeur particulaire de type coordonnées.
Definition VariableTypedef.h:383
Arcane::VariableItemReal
ItemVariableScalarRefT< Real > VariableItemReal
Grandeur de type réel.
Definition VariableTypedef.h:266
Arcane::VariableNodeReal3
MeshVariableScalarRefT< Node, Real3 > VariableNodeReal3
Grandeur au noeud de type coordonnées.
Definition VariableTypedef.h:359
Arcane::VariableItemReal3
ItemVariableScalarRefT< Real3 > VariableItemReal3
Grandeur de type coordonn?es 3D.
Definition VariableTypedef.h:353