1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
|
/*=========================================================================
Program: Visualization Toolkit
Module: TestAMRReadWrite.cxx
Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
All rights reserved.
See Copyright.txt or http://www.kitware.com/Copyright.htm for details.
This software is distributed WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE. See the above copyright notice for more information.
=========================================================================*/
// .NAME Test of vtkSimplePointsReader and vtkSimplePointsWriter
// .SECTION Description
//
#include "vtkSmartPointer.h"
#include "vtkSetGet.h"
#include "vtkOverlappingAMR.h"
#include "vtkAMRInformation.h"
#include "vtkUniformGrid.h"
#include "vtkNew.h"
#include "vtkAMREnzoReader.h"
#include "vtkTestUtilities.h"
#include "vtkCompositeDataWriter.h"
#include "vtkCompositeDataReader.h"
namespace
{
vtkSmartPointer<vtkOverlappingAMR> CreateTestAMR(int argc, char *argv[] )
{
char* fname = vtkTestUtilities::ExpandDataFileName(argc, argv,"Data/AMR/Enzo/DD0010/moving7_0010.hierarchy");
vtkNew<vtkAMREnzoReader> reader;
reader->SetFileName(fname);
delete [] fname;
reader->SetMaxLevel(8);
reader->SetCellArrayStatus( "TotalEnergy",1);
reader->Update();
vtkSmartPointer<vtkOverlappingAMR> ret = vtkOverlappingAMR::SafeDownCast(reader->GetOutputDataObject(0));
return ret;
}
}
int TestAMRReadWrite( int argc, char *argv[] )
{
vtkSmartPointer<vtkOverlappingAMR> amr = CreateTestAMR(argc, argv);
vtkNew<vtkCompositeDataWriter> writer;
writer->SetInputData(amr);
writer->SetFileName("testamr");
writer->Update();
vtkNew<vtkCompositeDataReader> reader;
reader->SetFileName("testamr");
reader->Update();
vtkSmartPointer<vtkOverlappingAMR> amr1 = vtkOverlappingAMR::SafeDownCast(reader->GetOutputDataObject(0));
int errors(0);
errors+= !( *amr1->GetAMRInfo() == *amr->GetAMRInfo());
return errors;
}
|
