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
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
|
// SPDX-FileCopyrightText: Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
// SPDX-FileCopyrightText: Copyright 2008 Sandia Corporation
// SPDX-License-Identifier: LicenseRef-BSD-3-Clause-Sandia-USGov
#include <vtkArrayInterpolate.h>
#include <vtkDenseArray.h>
#include <vtkSmartPointer.h>
#include <iostream>
#include <stdexcept>
void test_expression(const bool expression, const std::string& message)
{
if (!expression)
throw std::runtime_error(message);
}
int TestArrayInterpolationDense(int vtkNotUsed(argc), char* vtkNotUsed(argv)[])
{
try
{
vtkSmartPointer<vtkDenseArray<double>> a = vtkSmartPointer<vtkDenseArray<double>>::New();
a->Resize(4);
a->SetValue(0, 0);
a->SetValue(1, 1);
a->SetValue(2, 2);
a->SetValue(3, 3);
vtkSmartPointer<vtkDenseArray<double>> b = vtkSmartPointer<vtkDenseArray<double>>::New();
b->Resize(vtkArrayExtents(2));
vtkInterpolate(a.GetPointer(),
vtkArrayExtentsList(
vtkArrayExtents(vtkArrayRange(0, 1)), vtkArrayExtents(vtkArrayRange(1, 2))),
vtkArrayWeights(0.5, 0.5), vtkArrayExtents(vtkArrayRange(0, 1)), b.GetPointer());
vtkInterpolate(a.GetPointer(),
vtkArrayExtentsList(
vtkArrayExtents(vtkArrayRange(2, 3)), vtkArrayExtents(vtkArrayRange(3, 4))),
vtkArrayWeights(0.5, 0.5), vtkArrayExtents(vtkArrayRange(1, 2)), b.GetPointer());
test_expression(b->GetValue(0) == 0.5, "expected 0.5");
test_expression(b->GetValue(1) == 2.5, "expected 2.5");
vtkSmartPointer<vtkDenseArray<double>> c = vtkSmartPointer<vtkDenseArray<double>>::New();
c->Resize(4, 2);
c->SetValue(0, 0, 0);
c->SetValue(0, 1, 1);
c->SetValue(1, 0, 2);
c->SetValue(1, 1, 3);
c->SetValue(2, 0, 4);
c->SetValue(2, 1, 5);
c->SetValue(3, 0, 6);
c->SetValue(3, 1, 7);
vtkSmartPointer<vtkDenseArray<double>> d = vtkSmartPointer<vtkDenseArray<double>>::New();
d->Resize(vtkArrayExtents(2, 2));
vtkInterpolate(c.GetPointer(),
vtkArrayExtentsList(vtkArrayExtents(vtkArrayRange(0, 1), vtkArrayRange(0, 2)),
vtkArrayExtents(vtkArrayRange(1, 2), vtkArrayRange(0, 2))),
vtkArrayWeights(0.5, 0.5), vtkArrayExtents(vtkArrayRange(0, 1), vtkArrayRange(0, 2)),
d.GetPointer());
vtkInterpolate(c.GetPointer(),
vtkArrayExtentsList(vtkArrayExtents(vtkArrayRange(2, 3), vtkArrayRange(0, 2)),
vtkArrayExtents(vtkArrayRange(3, 4), vtkArrayRange(0, 2))),
vtkArrayWeights(0.5, 0.5), vtkArrayExtents(vtkArrayRange(1, 2), vtkArrayRange(0, 2)),
d.GetPointer());
test_expression(d->GetValue(0, 0) == 1, "expected 1");
test_expression(d->GetValue(0, 1) == 2, "expected 2");
test_expression(d->GetValue(1, 0) == 5, "expected 5");
test_expression(d->GetValue(1, 1) == 6, "expected 6");
return 0;
}
catch (std::exception& e)
{
cerr << e.what() << endl;
return 1;
}
}
|
