~ubuntu-branches/debian/sid/gmsh/sid

« back to all changes in this revision

Viewing changes to Plugin/CutGrid.cpp

  • Committer: Package Import Robot
  • Author(s): Anton Gladky
  • Date: 2015-08-26 22:45:27 UTC
  • mfrom: (1.3.25)
  • Revision ID: package-import@ubuntu.com-20150826224527-wiqnz3dwz6ap8tty
Tags: 2.10.1+dfsg1-1
* [4c450d1] Update d/watch.
* [90ca918] Imported Upstream version 2.10.1+dfsg1. (Closes: #793245)
* [af26665] Use any-arch instead of list of archs.
* [e94c6d2] Refresh patches.
* [6ab417f] Update d/copyright.

Show diffs side-by-side

added added

removed removed

Lines of Context:
44
44
  drawContext *ctx = (drawContext*)context;
45
45
 
46
46
  getPoint(0, 0, p);
47
 
  glRasterPos3d(p[0], p[1], p[2]);
48
 
  ctx->drawString("(X0, Y0, Z0)");
 
47
  ctx->drawString("(X0, Y0, Z0)", p[0], p[1], p[2]);
49
48
  if(getNbU() > 1){
50
49
    getPoint(getNbU() - 1, 0, p);
51
 
    glRasterPos3d(p[0], p[1], p[2]);
52
 
    ctx->drawString("(X1, Y1, Z1)");
 
50
    ctx->drawString("(X1, Y1, Z1)", p[0], p[1], p[2]);
53
51
  }
54
52
  if(getNbV() > 1){
55
53
    getPoint(0, getNbV() - 1, p);
56
 
    glRasterPos3d(p[0], p[1], p[2]);
57
 
    ctx->drawString("(X2, Y2, Z2)");
 
54
    ctx->drawString("(X2, Y2, Z2)", p[0], p[1], p[2]);
58
55
  }
59
56
 
60
57
  if(CutGridOptions_Number[11].def){
219
216
{
220
217
  double u = getNbU() > 1 ? (double)iU / (double)(getNbU() - 1.) : 0.;
221
218
  double v = getNbV() > 1 ? (double)iV / (double)(getNbV() - 1.) : 0.;
222
 
  X[0] = CutGridOptions_Number[0].def + 
 
219
  X[0] = CutGridOptions_Number[0].def +
223
220
    u  * (CutGridOptions_Number[3].def-CutGridOptions_Number[0].def) +
224
221
    v  * (CutGridOptions_Number[6].def-CutGridOptions_Number[0].def);
225
 
  X[1] = CutGridOptions_Number[1].def + 
 
222
  X[1] = CutGridOptions_Number[1].def +
226
223
    u  * (CutGridOptions_Number[4].def-CutGridOptions_Number[1].def) +
227
224
    v  * (CutGridOptions_Number[7].def-CutGridOptions_Number[1].def);
228
 
  X[2] = CutGridOptions_Number[2].def + 
 
225
  X[2] = CutGridOptions_Number[2].def +
229
226
    u  * (CutGridOptions_Number[5].def-CutGridOptions_Number[2].def) +
230
227
    v  * (CutGridOptions_Number[8].def-CutGridOptions_Number[2].def);
231
228
}
232
229
 
233
 
void GMSH_CutGridPlugin::addInView(int numsteps, int connect, int nbcomp, 
234
 
                                   double ***pnts, double ***vals, 
235
 
                                   std::vector<double> &P, int *nP, 
236
 
                                   std::vector<double> &L, int *nL, 
 
230
void GMSH_CutGridPlugin::addInView(int numsteps, int connect, int nbcomp,
 
231
                                   double ***pnts, double ***vals,
 
232
                                   std::vector<double> &P, int *nP,
 
233
                                   std::vector<double> &L, int *nL,
237
234
                                   std::vector<double> &Q, int *nQ)
238
235
{
239
236
  if(!connect || (getNbU() == 1 && getNbV() == 1)){ // generate points
240
 
    
 
237
 
241
238
    for(int i = 0; i < getNbU(); ++i){
242
239
      for(int j = 0; j < getNbV(); ++j){
243
240
        P.push_back(pnts[i][j][0]);
250
247
        }
251
248
      }
252
249
    }
253
 
    
 
250
 
254
251
  }
255
252
  else{ // generate lines or quads
256
 
    
 
253
 
257
254
    if(getNbU() == 1){
258
255
      for(int i = 0; i < getNbV()-1; ++i){
259
256
        L.push_back(pnts[0][i][0]); L.push_back(pnts[0][i+1][0]);
305
302
        }
306
303
      }
307
304
    }
308
 
  }     
 
305
  }
309
306
}
310
307
 
311
308
PView *GMSH_CutGridPlugin::GenerateView(PView *v1, int connect)
317
314
 
318
315
  PView *v2 = new PView();
319
316
  PViewDataList *data2 = getDataList(v2);
320
 
 
 
317
 
321
318
  OctreePost o(v1);
322
319
 
323
320
  int nbs = data1->getNumScalars();
337
334
      getPoint(i, j, pnts[i][j]);
338
335
    }
339
336
  }
340
 
  
 
337
 
341
338
  if(nbs){
342
339
    for(int i = 0; i < getNbU(); i++)
343
340
      for(int j = 0; j < getNbV(); j++)
344
341
        o.searchScalar(pnts[i][j][0], pnts[i][j][1], pnts[i][j][2], vals[i][j]);
345
 
    addInView(numsteps, connect, 1, pnts, vals, data2->SP, &data2->NbSP, 
 
342
    addInView(numsteps, connect, 1, pnts, vals, data2->SP, &data2->NbSP,
346
343
              data2->SL, &data2->NbSL, data2->SQ, &data2->NbSQ);
347
344
  }
348
345
 
358
355
    for(int i = 0; i < getNbU(); i++)
359
356
      for(int j = 0; j < getNbV(); j++)
360
357
        o.searchTensor(pnts[i][j][0], pnts[i][j][1], pnts[i][j][2], vals[i][j]);
361
 
    addInView(numsteps, connect, 9, pnts, vals, data2->TP, &data2->NbTP, 
 
358
    addInView(numsteps, connect, 9, pnts, vals, data2->TP, &data2->NbTP,
362
359
              data2->TL, &data2->NbTL, data2->TQ, &data2->NbTQ);
363
360
  }
364
361