-
Notifications
You must be signed in to change notification settings - Fork 0
/
vtkConeSlice.cpp
77 lines (62 loc) · 1.95 KB
/
vtkConeSlice.cpp
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
/*=========================================================================
Program: Visualization Toolkit
Module: vtkConeSlice.cpp
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.
=========================================================================*/
// Copyright 2014-2016 Etienne Tang
#include "vtkConeSlice.h"
#include "vtkObjectFactory.h"
#include "vtkInformation.h"
#include "vtkInformationVector.h"
#include "vtkPointSet.h"
#include "vtkPoints.h"
#include "vtkDoubleArray.h"
#include "vtkPointData.h"
#include "vtkMath.h"
#include <cmath>
vtkStandardNewMacro(vtkConeXrFunction)
vtkStandardNewMacro(vtkConeSlice)
vtkConeXrFunction::vtkConeXrFunction()
{
this->Xr1[0] = 0.;
this->Xr1[1] = 0.;
this->Xr2[0] = 0.;
this->Xr2[1] = 1.;
}
double vtkConeXrFunction::EvaluateFunction(double pt[3])
{
double r = sqrt(pt[1]*pt[1] + pt[2]*pt[2]);
return (Xr1[1]-Xr2[1])*(pt[0]-Xr1[0]) + (Xr2[0]-Xr1[0])*(r - Xr1[1]);
}
void vtkConeXrFunction::EvaluateGradient(double pt[3], double g[3])
{
double r = sqrt(pt[1]*pt[1] + pt[2]*pt[2]);
g[0] = Xr1[1]-Xr2[1];
g[1] = (Xr2[0]-Xr1[0])*pt[1]/r;
g[2] = (Xr2[0]-Xr1[0])*pt[2]/r;
}
void vtkConeXrFunction::PrintSelf( ostream& os, vtkIndent indent )
{
this->Superclass::PrintSelf( os, indent );
}
void vtkConeSlice::PrintSelf(ostream& os, vtkIndent indent)
{
this->Superclass::PrintSelf(os, indent);
}
vtkConeSlice::vtkConeSlice():
vtkCutter()
{
this->ConeFunction = vtkConeXrFunction::New();
this->SetXr1(0., 0.);
this->SetXr2(0., 1.);
this->SetCutFunction(this->ConeFunction);
}
vtkConeSlice::~vtkConeSlice()
{
this->ConeFunction->Delete();
}