vtkPV3FoamMesh.C 14.4 KB
Newer Older
1
2
3
4
/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
5
    \\  /    A nd           | Copyright (C) 1991-2010 OpenCFD Ltd.
6
7
8
9
10
     \\/     M anipulation  |
-------------------------------------------------------------------------------
License
    This file is part of OpenFOAM.

11
12
13
14
    OpenFOAM is free software: you can redistribute it and/or modify it
    under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.
15
16
17
18
19
20
21

    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
    for more details.

    You should have received a copy of the GNU General Public License
22
    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44

\*---------------------------------------------------------------------------*/

#include "vtkPV3Foam.H"

// Foam includes
#include "cellSet.H"
#include "faceSet.H"
#include "pointSet.H"
#include "fvMeshSubset.H"
#include "vtkPV3FoamReader.h"

// VTK includes
#include "vtkDataArraySelection.h"
#include "vtkMultiBlockDataSet.h"
#include "vtkPolyData.h"
#include "vtkUnstructuredGrid.h"

// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //

void Foam::vtkPV3Foam::convertMeshVolume
(
45
46
    vtkMultiBlockDataSet* output,
    int& blockNo
47
48
)
{
49
50
    arrayRange& range = arrayRangeVolume_;
    range.block(blockNo);      // set output block
Mark Olesen's avatar
Mark Olesen committed
51
52
53
54
    label datasetNo = 0;       // restart at dataset 0
    const fvMesh& mesh = *meshPtr_;

    // resize for decomposed polyhedra
55
    regionPolyDecomp_.setSize(range.size());
Mark Olesen's avatar
Mark Olesen committed
56

57
58
    if (debug)
    {
59
        Info<< "<beg> Foam::vtkPV3Foam::convertMeshVolume" << endl;
Mark Olesen's avatar
Mark Olesen committed
60
        printMemory();
61
62
    }

Mark Olesen's avatar
Mark Olesen committed
63
    // Convert the internalMesh
64
    // this looks like more than one part, but it isn't
65
    for (int partId = range.start(); partId < range.end(); ++partId)
Mark Olesen's avatar
Mark Olesen committed
66
67
    {
        const word partName = "internalMesh";
68

Mark Olesen's avatar
Mark Olesen committed
69
        if (!partStatus_[partId])
70
        {
Mark Olesen's avatar
Mark Olesen committed
71
            continue;
72
73
        }

Mark Olesen's avatar
Mark Olesen committed
74
75
76
        vtkUnstructuredGrid* vtkmesh = volumeVTKMesh
        (
            mesh,
Mark Olesen's avatar
Mark Olesen committed
77
            regionPolyDecomp_[datasetNo]
Mark Olesen's avatar
Mark Olesen committed
78
        );
79

Mark Olesen's avatar
Mark Olesen committed
80
81
        if (vtkmesh)
        {
82
            AddToBlock(output, vtkmesh, range, datasetNo, partName);
Mark Olesen's avatar
Mark Olesen committed
83
            vtkmesh->Delete();
84

Mark Olesen's avatar
Mark Olesen committed
85
            partDataset_[partId] = datasetNo++;
Mark Olesen's avatar
Mark Olesen committed
86
        }
87
88
    }

Mark Olesen's avatar
Mark Olesen committed
89
    // anything added?
Mark Olesen's avatar
Mark Olesen committed
90
    if (datasetNo)
91
92
    {
        ++blockNo;
93
    }
94
95
96
97

    if (debug)
    {
        Info<< "<end> Foam::vtkPV3Foam::convertMeshVolume" << endl;
Mark Olesen's avatar
Mark Olesen committed
98
        printMemory();
99
    }
100
101
102
103
104
}


void Foam::vtkPV3Foam::convertMeshLagrangian
(
105
106
    vtkMultiBlockDataSet* output,
    int& blockNo
107
108
)
{
109
110
    arrayRange& range = arrayRangeLagrangian_;
    range.block(blockNo);      // set output block
Mark Olesen's avatar
Mark Olesen committed
111
112
113
    label datasetNo = 0;       // restart at dataset 0
    const fvMesh& mesh = *meshPtr_;

114
115
    if (debug)
    {
116
        Info<< "<beg> Foam::vtkPV3Foam::convertMeshLagrangian" << endl;
Mark Olesen's avatar
Mark Olesen committed
117
        printMemory();
118
119
    }

120
    for (int partId = range.start(); partId < range.end(); ++partId)
Mark Olesen's avatar
Mark Olesen committed
121
122
123
124
    {
        const word cloudName = getPartName(partId);

        if (!partStatus_[partId])
125
126
127
128
        {
            continue;
        }

Mark Olesen's avatar
Mark Olesen committed
129
        vtkPolyData* vtkmesh = lagrangianVTKMesh(mesh, cloudName);
Mark Olesen's avatar
Mark Olesen committed
130

Mark Olesen's avatar
Mark Olesen committed
131
132
        if (vtkmesh)
        {
133
            AddToBlock(output, vtkmesh, range, datasetNo, cloudName);
Mark Olesen's avatar
Mark Olesen committed
134
            vtkmesh->Delete();
135

Mark Olesen's avatar
Mark Olesen committed
136
            partDataset_[partId] = datasetNo++;
Mark Olesen's avatar
Mark Olesen committed
137
        }
138
    }
139

Mark Olesen's avatar
Mark Olesen committed
140
    // anything added?
Mark Olesen's avatar
Mark Olesen committed
141
    if (datasetNo)
142
143
144
145
    {
        ++blockNo;
    }

146
147
148
    if (debug)
    {
        Info<< "<end> Foam::vtkPV3Foam::convertMeshLagrangian" << endl;
Mark Olesen's avatar
Mark Olesen committed
149
        printMemory();
150
    }
151
152
153
154
155
}


void Foam::vtkPV3Foam::convertMeshPatches
(
156
157
    vtkMultiBlockDataSet* output,
    int& blockNo
158
159
)
{
160
161
    arrayRange& range = arrayRangePatches_;
    range.block(blockNo);      // set output block
Mark Olesen's avatar
Mark Olesen committed
162
163
164
165
    label datasetNo = 0;       // restart at dataset 0
    const fvMesh& mesh = *meshPtr_;
    const polyBoundaryMesh& patches = mesh.boundaryMesh();

166
167
    if (debug)
    {
168
        Info<< "<beg> Foam::vtkPV3Foam::convertMeshPatches" << endl;
Mark Olesen's avatar
Mark Olesen committed
169
        printMemory();
170
171
    }

172
    for (int partId = range.start(); partId < range.end(); ++partId)
173
    {
Mark Olesen's avatar
Mark Olesen committed
174
175
        const word patchName = getPartName(partId);
        const label  patchId = patches.findPatchID(patchName);
176

Mark Olesen's avatar
Mark Olesen committed
177
        if (!partStatus_[partId] || patchId < 0)
178
        {
Mark Olesen's avatar
Mark Olesen committed
179
180
            continue;
        }
181

Mark Olesen's avatar
Mark Olesen committed
182
183
184
185
186
        if (debug)
        {
            Info<< "Creating VTK mesh for patch[" << patchId <<"] "
                << patchName  << endl;
        }
187

Mark Olesen's avatar
Mark Olesen committed
188
        vtkPolyData* vtkmesh = patchVTKMesh(patches[patchId]);
189

Mark Olesen's avatar
Mark Olesen committed
190
191
        if (vtkmesh)
        {
192
            AddToBlock(output, vtkmesh, range, datasetNo, patchName);
Mark Olesen's avatar
Mark Olesen committed
193
            vtkmesh->Delete();
194

Mark Olesen's avatar
Mark Olesen committed
195
            partDataset_[partId] = datasetNo++;
196
197
        }
    }
198

Mark Olesen's avatar
Mark Olesen committed
199
    // anything added?
Mark Olesen's avatar
Mark Olesen committed
200
    if (datasetNo)
201
202
203
204
    {
        ++blockNo;
    }

205
206
207
    if (debug)
    {
        Info<< "<end> Foam::vtkPV3Foam::convertMeshPatches" << endl;
Mark Olesen's avatar
Mark Olesen committed
208
        printMemory();
209
    }
210
211
212
213
214
}


void Foam::vtkPV3Foam::convertMeshCellZones
(
215
216
    vtkMultiBlockDataSet* output,
    int& blockNo
217
218
)
{
219
220
    arrayRange& range = arrayRangeCellZones_;
    range.block(blockNo);      // set output block
Mark Olesen's avatar
Mark Olesen committed
221
222
223
224
    label datasetNo = 0;       // restart at dataset 0
    const fvMesh& mesh = *meshPtr_;

    // resize for decomposed polyhedra
225
    zonePolyDecomp_.setSize(range.size());
Mark Olesen's avatar
Mark Olesen committed
226

227
    if (range.empty())
Mark Olesen's avatar
Mark Olesen committed
228
229
230
231
    {
        return;
    }

232
233
    if (debug)
    {
234
        Info<< "<beg> Foam::vtkPV3Foam::convertMeshCellZones" << endl;
Mark Olesen's avatar
Mark Olesen committed
235
        printMemory();
236
237
    }

Mark Olesen's avatar
Mark Olesen committed
238
    const cellZoneMesh& zMesh = mesh.cellZones();
239
    for (int partId = range.start(); partId < range.end(); ++partId)
240
    {
Mark Olesen's avatar
Mark Olesen committed
241
242
        const word zoneName = getPartName(partId);
        const label  zoneId = zMesh.findZoneID(zoneName);
243

Mark Olesen's avatar
Mark Olesen committed
244
        if (!partStatus_[partId] || zoneId < 0)
245
        {
Mark Olesen's avatar
Mark Olesen committed
246
247
            continue;
        }
Mark Olesen's avatar
Mark Olesen committed
248

Mark Olesen's avatar
Mark Olesen committed
249
250
251
252
253
        if (debug)
        {
            Info<< "Creating VTK mesh for cellZone[" << zoneId << "] "
                << zoneName << endl;
        }
254

Mark Olesen's avatar
Mark Olesen committed
255
256
        fvMeshSubset subsetter(mesh);
        subsetter.setLargeCellSubset(zMesh[zoneId]);
257

Mark Olesen's avatar
Mark Olesen committed
258
259
260
261
262
        vtkUnstructuredGrid* vtkmesh = volumeVTKMesh
        (
            subsetter.subMesh(),
            zonePolyDecomp_[datasetNo]
        );
263

Mark Olesen's avatar
Mark Olesen committed
264
265
266
267
        if (vtkmesh)
        {
            // superCells + addPointCellLabels must contain global cell ids
            inplaceRenumber
268
            (
Mark Olesen's avatar
Mark Olesen committed
269
270
271
272
273
274
275
                subsetter.cellMap(),
                zonePolyDecomp_[datasetNo].superCells()
            );
            inplaceRenumber
            (
                subsetter.cellMap(),
                zonePolyDecomp_[datasetNo].addPointCellLabels()
276
277
            );

Mark Olesen's avatar
Mark Olesen committed
278
279
            // copy pointMap as well, otherwise pointFields fail
            zonePolyDecomp_[datasetNo].pointMap() = subsetter.pointMap();
Mark Olesen's avatar
Mark Olesen committed
280

281
            AddToBlock(output, vtkmesh, range, datasetNo, zoneName);
Mark Olesen's avatar
Mark Olesen committed
282
283
284
            vtkmesh->Delete();

            partDataset_[partId] = datasetNo++;
285
286
        }
    }
287

Mark Olesen's avatar
Mark Olesen committed
288
    // anything added?
Mark Olesen's avatar
Mark Olesen committed
289
    if (datasetNo)
290
291
292
293
    {
        ++blockNo;
    }

294
295
296
    if (debug)
    {
        Info<< "<end> Foam::vtkPV3Foam::convertMeshCellZones" << endl;
Mark Olesen's avatar
Mark Olesen committed
297
        printMemory();
298
    }
299
300
301
}


302
void Foam::vtkPV3Foam::convertMeshCellSets
303
(
304
305
    vtkMultiBlockDataSet* output,
    int& blockNo
306
307
)
{
308
309
    arrayRange& range = arrayRangeCellSets_;
    range.block(blockNo);      // set output block
Mark Olesen's avatar
Mark Olesen committed
310
311
312
313
    label datasetNo = 0;       // restart at dataset 0
    const fvMesh& mesh = *meshPtr_;

    // resize for decomposed polyhedra
314
    csetPolyDecomp_.setSize(range.size());
Mark Olesen's avatar
Mark Olesen committed
315

316
317
    if (debug)
    {
318
        Info<< "<beg> Foam::vtkPV3Foam::convertMeshCellSets" << endl;
Mark Olesen's avatar
Mark Olesen committed
319
        printMemory();
320
321
    }

322
    for (int partId = range.start(); partId < range.end(); ++partId)
323
    {
Mark Olesen's avatar
Mark Olesen committed
324
325
326
        const word partName = getPartName(partId);

        if (!partStatus_[partId])
327
        {
Mark Olesen's avatar
Mark Olesen committed
328
329
            continue;
        }
330

Mark Olesen's avatar
Mark Olesen committed
331
332
333
334
        if (debug)
        {
            Info<< "Creating VTK mesh for cellSet=" << partName << endl;
        }
335

Mark Olesen's avatar
Mark Olesen committed
336
337
338
        const cellSet cSet(mesh, partName);
        fvMeshSubset subsetter(mesh);
        subsetter.setLargeCellSubset(cSet);
339

Mark Olesen's avatar
Mark Olesen committed
340
341
342
343
344
        vtkUnstructuredGrid* vtkmesh = volumeVTKMesh
        (
            subsetter.subMesh(),
            csetPolyDecomp_[datasetNo]
        );
345

Mark Olesen's avatar
Mark Olesen committed
346
347
348
349
350
351
352
353
354
        if (vtkmesh)
        {
            // superCells + addPointCellLabels must contain global cell ids
            inplaceRenumber
            (
                subsetter.cellMap(),
                csetPolyDecomp_[datasetNo].superCells()
            );
            inplaceRenumber
355
            (
Mark Olesen's avatar
Mark Olesen committed
356
357
                subsetter.cellMap(),
                csetPolyDecomp_[datasetNo].addPointCellLabels()
358
359
            );

Mark Olesen's avatar
Mark Olesen committed
360
361
            // copy pointMap as well, otherwise pointFields fail
            csetPolyDecomp_[datasetNo].pointMap() = subsetter.pointMap();
Mark Olesen's avatar
Mark Olesen committed
362

363
            AddToBlock(output, vtkmesh, range, datasetNo, partName);
Mark Olesen's avatar
Mark Olesen committed
364
365
366
            vtkmesh->Delete();

            partDataset_[partId] = datasetNo++;
367
368
        }
    }
369

Mark Olesen's avatar
Mark Olesen committed
370
    // anything added?
Mark Olesen's avatar
Mark Olesen committed
371
    if (datasetNo)
372
373
374
375
    {
        ++blockNo;
    }

376
377
    if (debug)
    {
378
        Info<< "<end> Foam::vtkPV3Foam::convertMeshCellSets" << endl;
Mark Olesen's avatar
Mark Olesen committed
379
        printMemory();
380
    }
381
382
}

383

384
385
void Foam::vtkPV3Foam::convertMeshFaceZones
(
386
387
    vtkMultiBlockDataSet* output,
    int& blockNo
388
389
)
{
390
391
    arrayRange& range = arrayRangeFaceZones_;
    range.block(blockNo);      // set output block
Mark Olesen's avatar
Mark Olesen committed
392
393
394
    label datasetNo = 0;       // restart at dataset 0
    const fvMesh& mesh = *meshPtr_;

395
    if (range.empty())
Mark Olesen's avatar
Mark Olesen committed
396
397
398
399
    {
        return;
    }

400
401
    if (debug)
    {
402
        Info<< "<beg> Foam::vtkPV3Foam::convertMeshFaceZones" << endl;
Mark Olesen's avatar
Mark Olesen committed
403
        printMemory();
404
405
    }

Mark Olesen's avatar
Mark Olesen committed
406
    const faceZoneMesh& zMesh = mesh.faceZones();
407
    for (int partId = range.start(); partId < range.end(); ++partId)
408
    {
Mark Olesen's avatar
Mark Olesen committed
409
410
        const word zoneName = getPartName(partId);
        const label  zoneId = zMesh.findZoneID(zoneName);
411

Mark Olesen's avatar
Mark Olesen committed
412
        if (!partStatus_[partId] || zoneId < 0)
413
        {
Mark Olesen's avatar
Mark Olesen committed
414
415
            continue;
        }
416

Mark Olesen's avatar
Mark Olesen committed
417
418
419
420
421
        if (debug)
        {
            Info<< "Creating VTKmesh for faceZone[" << zoneId << "] "
                << zoneName << endl;
        }
422

Mark Olesen's avatar
Mark Olesen committed
423
424
425
        vtkPolyData* vtkmesh = faceZoneVTKMesh(mesh, zMesh[zoneId]);
        if (vtkmesh)
        {
426
            AddToBlock(output, vtkmesh, range, datasetNo, zoneName);
Mark Olesen's avatar
Mark Olesen committed
427
            vtkmesh->Delete();
Mark Olesen's avatar
Mark Olesen committed
428

Mark Olesen's avatar
Mark Olesen committed
429
            partDataset_[partId] = datasetNo++;
430
431
        }
    }
432

Mark Olesen's avatar
Mark Olesen committed
433
    // anything added?
Mark Olesen's avatar
Mark Olesen committed
434
    if (datasetNo)
435
436
437
438
    {
        ++blockNo;
    }

439
440
441
    if (debug)
    {
        Info<< "<end> Foam::vtkPV3Foam::convertMeshFaceZones" << endl;
Mark Olesen's avatar
Mark Olesen committed
442
        printMemory();
443
    }
444
445
446
}


447
void Foam::vtkPV3Foam::convertMeshFaceSets
448
(
449
450
    vtkMultiBlockDataSet* output,
    int& blockNo
451
452
)
{
453
454
    arrayRange& range = arrayRangeFaceSets_;
    range.block(blockNo);      // set output block
Mark Olesen's avatar
Mark Olesen committed
455
456
457
    label datasetNo = 0;       // restart at dataset 0
    const fvMesh& mesh = *meshPtr_;

458
459
    if (debug)
    {
460
        Info<< "<beg> Foam::vtkPV3Foam::convertMeshFaceSets" << endl;
Mark Olesen's avatar
Mark Olesen committed
461
        printMemory();
462
463
    }

464
    for (int partId = range.start(); partId < range.end(); ++partId)
465
    {
Mark Olesen's avatar
Mark Olesen committed
466
        const word partName = getPartName(partId);
467

Mark Olesen's avatar
Mark Olesen committed
468
469
470
471
        if (!partStatus_[partId])
        {
            continue;
        }
472

Mark Olesen's avatar
Mark Olesen committed
473
474
475
476
        if (debug)
        {
            Info<< "Creating VTK mesh for faceSet=" << partName << endl;
        }
477

Mark Olesen's avatar
Mark Olesen committed
478
        const faceSet fSet(mesh, partName);
479

Mark Olesen's avatar
Mark Olesen committed
480
481
482
        vtkPolyData* vtkmesh = faceSetVTKMesh(mesh, fSet);
        if (vtkmesh)
        {
483
            AddToBlock(output, vtkmesh, range, datasetNo, partName);
Mark Olesen's avatar
Mark Olesen committed
484
            vtkmesh->Delete();
Mark Olesen's avatar
Mark Olesen committed
485

Mark Olesen's avatar
Mark Olesen committed
486
            partDataset_[partId] = datasetNo++;
487
488
        }
    }
489

Mark Olesen's avatar
Mark Olesen committed
490
    // anything added?
Mark Olesen's avatar
Mark Olesen committed
491
    if (datasetNo)
492
493
494
495
    {
        ++blockNo;
    }

496
497
    if (debug)
    {
498
        Info<< "<end> Foam::vtkPV3Foam::convertMeshFaceSets" << endl;
Mark Olesen's avatar
Mark Olesen committed
499
        printMemory();
500
    }
501
502
503
504
505
}


void Foam::vtkPV3Foam::convertMeshPointZones
(
506
507
    vtkMultiBlockDataSet* output,
    int& blockNo
508
509
)
{
510
511
    arrayRange& range = arrayRangePointZones_;
    range.block(blockNo);      // set output block
Mark Olesen's avatar
Mark Olesen committed
512
    label datasetNo = 0;       // restart at dataset 0
513
    const fvMesh& mesh = *meshPtr_;
Mark Olesen's avatar
Mark Olesen committed
514

515
516
    if (debug)
    {
517
        Info<< "<beg> Foam::vtkPV3Foam::convertMeshPointZones" << endl;
Mark Olesen's avatar
Mark Olesen committed
518
        printMemory();
519
520
    }

521
    if (range.size())
522
    {
Mark Olesen's avatar
Mark Olesen committed
523
        const pointZoneMesh& zMesh = mesh.pointZones();
524
        for (int partId = range.start(); partId < range.end(); ++partId)
525
        {
Mark Olesen's avatar
Mark Olesen committed
526
            word zoneName = getPartName(partId);
Mark Olesen's avatar
Mark Olesen committed
527
            label zoneId = zMesh.findZoneID(zoneName);
528

Mark Olesen's avatar
Mark Olesen committed
529
            if (!partStatus_[partId] || zoneId < 0)
530
531
532
533
            {
                continue;
            }

Mark Olesen's avatar
Mark Olesen committed
534
535
            vtkPolyData* vtkmesh = pointZoneVTKMesh(mesh, zMesh[zoneId]);
            if (vtkmesh)
536
            {
537
                AddToBlock(output, vtkmesh, range, datasetNo, zoneName);
Mark Olesen's avatar
Mark Olesen committed
538
539
                vtkmesh->Delete();

Mark Olesen's avatar
Mark Olesen committed
540
                partDataset_[partId] = datasetNo++;
541
542
543
            }
        }
    }
544

Mark Olesen's avatar
Mark Olesen committed
545
    // anything added?
Mark Olesen's avatar
Mark Olesen committed
546
    if (datasetNo)
547
548
549
550
    {
        ++blockNo;
    }

551
552
553
    if (debug)
    {
        Info<< "<end> Foam::vtkPV3Foam::convertMeshPointZones" << endl;
Mark Olesen's avatar
Mark Olesen committed
554
        printMemory();
555
    }
556
557
558
559
}



560
void Foam::vtkPV3Foam::convertMeshPointSets
561
(
562
563
    vtkMultiBlockDataSet* output,
    int& blockNo
564
565
)
{
566
567
    arrayRange& range = arrayRangePointSets_;
    range.block(blockNo);      // set output block
Mark Olesen's avatar
Mark Olesen committed
568
569
570
    label datasetNo = 0;       // restart at dataset 0
    const fvMesh& mesh = *meshPtr_;

571
572
    if (debug)
    {
573
        Info<< "<beg> Foam::vtkPV3Foam::convertMeshPointSets" << endl;
Mark Olesen's avatar
Mark Olesen committed
574
        printMemory();
575
576
    }

577
    for (int partId = range.start(); partId < range.end(); ++partId)
578
    {
Mark Olesen's avatar
Mark Olesen committed
579
        word partName = getPartName(partId);
580

Mark Olesen's avatar
Mark Olesen committed
581
582
583
584
        if (!partStatus_[partId])
        {
            continue;
        }
585

Mark Olesen's avatar
Mark Olesen committed
586
587
588
589
        if (debug)
        {
            Info<< "Creating VTK mesh for pointSet=" << partName << endl;
        }
590

Mark Olesen's avatar
Mark Olesen committed
591
        const pointSet pSet(mesh, partName);
592

Mark Olesen's avatar
Mark Olesen committed
593
594
595
        vtkPolyData* vtkmesh = pointSetVTKMesh(mesh, pSet);
        if (vtkmesh)
        {
596
            AddToBlock(output, vtkmesh, range, datasetNo, partName);
Mark Olesen's avatar
Mark Olesen committed
597
            vtkmesh->Delete();
Mark Olesen's avatar
Mark Olesen committed
598

Mark Olesen's avatar
Mark Olesen committed
599
            partDataset_[partId] = datasetNo++;
600
601
        }
    }
602

Mark Olesen's avatar
Mark Olesen committed
603
    // anything added?
Mark Olesen's avatar
Mark Olesen committed
604
    if (datasetNo)
605
606
607
608
    {
        ++blockNo;
    }

609
610
    if (debug)
    {
611
        Info<< "<end> Foam::vtkPV3Foam::convertMeshPointSets" << endl;
Mark Olesen's avatar
Mark Olesen committed
612
        printMemory();
613
    }
614
615
}

616

617
// ************************************************************************* //