| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 // G4VTwistSurface 26 // G4VTwistSurface
27 // 27 //
28 // Class description: 28 // Class description:
29 // 29 //
30 // Abstract base class for boundary surface of 30 // Abstract base class for boundary surface of G4VSolid.
31 31
32 // 01-Aug-2002 - Kotoyo Hoshina (hoshina@hepbu 32 // 01-Aug-2002 - Kotoyo Hoshina (hoshina@hepburn.s.chiba-u.ac.jp), created.
33 // 13-Nov-2003 - O.Link (Oliver.Link@cern.ch), 33 // 13-Nov-2003 - O.Link (Oliver.Link@cern.ch), Integration in Geant4
34 // from original version in Jupi 34 // from original version in Jupiter-2.5.02 application.
35 // ------------------------------------------- 35 // --------------------------------------------------------------------
36 #ifndef G4VTWISTSURFACE_HH 36 #ifndef G4VTWISTSURFACE_HH
37 #define G4VTWISTSURFACE_HH 37 #define G4VTWISTSURFACE_HH
38 38
39 #include <CLHEP/Units/SystemOfUnits.h> 39 #include <CLHEP/Units/SystemOfUnits.h>
40 40
41 #include "G4VSolid.hh" 41 #include "G4VSolid.hh"
42 #include "geomdefs.hh" 42 #include "geomdefs.hh"
43 43
44 #include "G4RotationMatrix.hh" 44 #include "G4RotationMatrix.hh"
45 45
46 #define G4VSURFACENXX 10 46 #define G4VSURFACENXX 10
47 47
48 class G4VTwistSurface 48 class G4VTwistSurface
49 { 49 {
50 public: << 50 public: // without description
51 51
52 enum EValidate { kDontValidate = 0, kValida 52 enum EValidate { kDontValidate = 0, kValidateWithTol = 1,
53 kValidateWithoutTol = 2, k 53 kValidateWithoutTol = 2, kUninitialized = 3 };
54 54
>> 55 public: // with description
>> 56
55 G4VTwistSurface (const G4String& name); 57 G4VTwistSurface (const G4String& name);
56 G4VTwistSurface (const G4String& name, 58 G4VTwistSurface (const G4String& name,
57 const G4RotationMatrix& ro << 59 const G4RotationMatrix& rot,
58 const G4ThreeVector& tl << 60 const G4ThreeVector& tlate,
59 G4int ha << 61 G4int handedness,
60 const EAxis ax << 62 const EAxis axis1,
61 const EAxis ax << 63 const EAxis axis2,
62 G4double ax << 64 G4double axis0min = -kInfinity,
63 G4double ax << 65 G4double axis1min = -kInfinity,
64 G4double ax << 66 G4double axis0max = kInfinity,
65 G4double ax << 67 G4double axis1max = kInfinity);
66 68
67 virtual ~G4VTwistSurface() = default; << 69 virtual ~G4VTwistSurface();
68 70
69 virtual G4int AmIOnLeftSide(const G4Thr 71 virtual G4int AmIOnLeftSide(const G4ThreeVector& me,
70 const G4Thr 72 const G4ThreeVector& vec,
71 G4boo 73 G4bool withTol = true);
72 74
73 virtual G4double DistanceToBoundary( 75 virtual G4double DistanceToBoundary( G4int areacode,
74 76 G4ThreeVector& xx,
75 const 77 const G4ThreeVector& p) ;
76 78
77 virtual G4double DistanceToIn(const G4Thre 79 virtual G4double DistanceToIn(const G4ThreeVector& gp,
78 const G4Thre 80 const G4ThreeVector& gv,
79 G4Thre 81 G4ThreeVector& gxxbest);
80 virtual G4double DistanceToOut(const G4Thr 82 virtual G4double DistanceToOut(const G4ThreeVector& gp,
81 const G4Thr 83 const G4ThreeVector& gv,
82 G4Thr 84 G4ThreeVector& gxxbest);
83 virtual G4double DistanceTo(const G4ThreeV 85 virtual G4double DistanceTo(const G4ThreeVector& gp,
84 G4ThreeV 86 G4ThreeVector& gxx);
85 87
86 virtual G4int DistanceToSurface(const G 88 virtual G4int DistanceToSurface(const G4ThreeVector& gp,
87 const G 89 const G4ThreeVector& gv,
88 G 90 G4ThreeVector gxx[],
89 G 91 G4double distance[],
90 G 92 G4int areacode[],
91 G 93 G4bool isvalid[],
92 EValida 94 EValidate validate=kValidateWithTol) = 0;
93 95
94 virtual G4int DistanceToSurface(const G 96 virtual G4int DistanceToSurface(const G4ThreeVector& gp,
95 G 97 G4ThreeVector gxx[],
96 G 98 G4double distance[],
97 G 99 G4int areacode[]) = 0;
98 100
99 void DebugPrint() const; 101 void DebugPrint() const;
100 102
>> 103 // get methods
>> 104
101 virtual G4ThreeVector GetNormal(const G4Thr 105 virtual G4ThreeVector GetNormal(const G4ThreeVector& xx,G4bool isGlobal) = 0;
102 106
103 virtual G4String GetName() const { ret 107 virtual G4String GetName() const { return fName; }
104 virtual void GetBoundaryParameters 108 virtual void GetBoundaryParameters(const G4int& areacode,
105 109 G4ThreeVector& d,
106 110 G4ThreeVector& x0,
107 111 G4int& boundarytype) const;
108 virtual G4ThreeVector GetBoundaryAtPZ(G4int 112 virtual G4ThreeVector GetBoundaryAtPZ(G4int areacode,
109 const 113 const G4ThreeVector& p) const;
110 114
111 inline G4double DistanceToPlaneWithV(const 115 inline G4double DistanceToPlaneWithV(const G4ThreeVector& p,
112 const 116 const G4ThreeVector& v,
113 const 117 const G4ThreeVector& x0,
114 const 118 const G4ThreeVector& n0,
115 119 G4ThreeVector& xx);
116 120
117 inline G4double DistanceToPlane(const G4Th 121 inline G4double DistanceToPlane(const G4ThreeVector& p,
118 const G4Th 122 const G4ThreeVector& x0,
119 const G4Th 123 const G4ThreeVector& n0,
120 G4Th 124 G4ThreeVector& xx);
121 125
122 inline G4double DistanceToPlane(const G4Th 126 inline G4double DistanceToPlane(const G4ThreeVector& p,
123 const G4Th 127 const G4ThreeVector& x0,
124 const G4Th 128 const G4ThreeVector& t1,
125 const G4Th 129 const G4ThreeVector& t2,
126 G4Th 130 G4ThreeVector& xx,
127 G4Th 131 G4ThreeVector& n);
128 132
129 inline G4double DistanceToLine (const G4Th 133 inline G4double DistanceToLine (const G4ThreeVector& p,
130 const G4Th 134 const G4ThreeVector& x0,
131 const G4Th 135 const G4ThreeVector& d,
132 G4Th 136 G4ThreeVector& xx);
133 137
134 inline G4bool IsAxis0 (G4int areacode) c 138 inline G4bool IsAxis0 (G4int areacode) const;
135 inline G4bool IsAxis1 (G4int areacode) c 139 inline G4bool IsAxis1 (G4int areacode) const;
136 inline G4bool IsOutside (G4int areacode) c 140 inline G4bool IsOutside (G4int areacode) const;
137 inline G4bool IsInside (G4int areacode, G 141 inline G4bool IsInside (G4int areacode, G4bool testbitmode = false) const;
138 inline G4bool IsBoundary (G4int areacode, G 142 inline G4bool IsBoundary (G4int areacode, G4bool testbitmode = false) const;
139 inline G4bool IsCorner (G4int areacode, G 143 inline G4bool IsCorner (G4int areacode, G4bool testbitmode = false) const;
140 inline G4bool IsValidNorm() const { return 144 inline G4bool IsValidNorm() const { return fIsValidNorm; }
141 G4bool IsSameBoundary (G4VTwistSurface* sur 145 G4bool IsSameBoundary (G4VTwistSurface* surface1, G4int areacode1,
142 G4VTwistSurface* sur 146 G4VTwistSurface* surface2, G4int areacode2 ) const;
143 inline G4int GetAxisType(G4int areacode, G 147 inline G4int GetAxisType(G4int areacode, G4int whichaxis) const;
144 148
145 inline G4ThreeVector ComputeGlobalPoint 149 inline G4ThreeVector ComputeGlobalPoint (const G4ThreeVector& lp) const;
146 inline G4ThreeVector ComputeLocalPoint 150 inline G4ThreeVector ComputeLocalPoint (const G4ThreeVector& gp) const;
147 inline G4ThreeVector ComputeGlobalDirection 151 inline G4ThreeVector ComputeGlobalDirection (const G4ThreeVector& lp) const;
148 inline G4ThreeVector ComputeLocalDirection 152 inline G4ThreeVector ComputeLocalDirection (const G4ThreeVector& gp) const;
149 153
150 // set methods 154 // set methods
151 155
152 inline void SetAxis(G4int i, const EAxis ax 156 inline void SetAxis(G4int i, const EAxis axis) { fAxis[i] = axis; }
153 inline void SetNeighbours(G4VTwistSurface* 157 inline void SetNeighbours(G4VTwistSurface* ax0min, G4VTwistSurface* ax1min,
154 G4VTwistSurface* 158 G4VTwistSurface* ax0max, G4VTwistSurface* ax1max);
155 159
156 virtual G4ThreeVector SurfacePoint(G4double 160 virtual G4ThreeVector SurfacePoint(G4double, G4double,
157 G4bool i 161 G4bool isGlobal = false ) = 0 ;
158 virtual G4double GetBoundaryMin(G4double) = 162 virtual G4double GetBoundaryMin(G4double) = 0 ;
159 virtual G4double GetBoundaryMax(G4double) = 163 virtual G4double GetBoundaryMax(G4double) = 0 ;
160 virtual G4double GetSurfaceArea() = 0 ; 164 virtual G4double GetSurfaceArea() = 0 ;
161 virtual void GetFacets(G4int m, G4int n, G4 165 virtual void GetFacets(G4int m, G4int n, G4double xyz[][3],
162 G4int faces[][4], G4 166 G4int faces[][4], G4int iside) = 0 ;
163 G4int GetNode( G4int i, G4int j, G4int m, G 167 G4int GetNode( G4int i, G4int j, G4int m, G4int n, G4int iside ) ;
164 G4int GetFace( G4int i, G4int j, G4int m, G 168 G4int GetFace( G4int i, G4int j, G4int m, G4int n, G4int iside ) ;
165 G4int GetEdgeVisibility( G4int i, G4int j, 169 G4int GetEdgeVisibility( G4int i, G4int j, G4int m, G4int n,
166 G4int number, G4in 170 G4int number, G4int orientation) ;
167 171
>> 172 public: // without description
>> 173
168 G4VTwistSurface(__void__&); 174 G4VTwistSurface(__void__&);
169 // Fake default constructor for usage res 175 // Fake default constructor for usage restricted to direct object
170 // persistency for clients requiring prea 176 // persistency for clients requiring preallocation of memory for
171 // persistifiable objects. 177 // persistifiable objects.
172 178
173 protected: << 179 protected: // with description
174 180
>> 181 // get methods
>> 182
175 inline G4VTwistSurface** GetNeighbours() { 183 inline G4VTwistSurface** GetNeighbours() { return fNeighbours; }
176 inline G4int GetNeighbours(G4int areacode, 184 inline G4int GetNeighbours(G4int areacode, G4VTwistSurface* surfaces[]);
177 inline G4ThreeVector GetCorner(G4int areaco 185 inline G4ThreeVector GetCorner(G4int areacode) const;
178 void GetBoundaryAxis(G4int areacode, EAxis << 186 void GetBoundaryAxis(G4int areacode, EAxis axis[]) const;
179 void GetBoundaryLimit(G4int areacode, G4dou << 187 void GetBoundaryLimit(G4int areacode, G4double limit[]) const;
180 virtual G4int GetAreaCode(const G4ThreeVect 188 virtual G4int GetAreaCode(const G4ThreeVector& xx, G4bool withtol=true) = 0;
>> 189
>> 190 // set methods
181 191
182 virtual void SetBoundary(const G4int& 192 virtual void SetBoundary(const G4int& axiscode,
183 const G4ThreeVecto 193 const G4ThreeVector& direction,
184 const G4ThreeVecto 194 const G4ThreeVector& x0,
185 const G4int& 195 const G4int& boundarytype);
186 // areacode must be one of them: 196 // areacode must be one of them:
187 // sAxis0 & sAxisMin, sAxis0 & sAxisMax, 197 // sAxis0 & sAxisMin, sAxis0 & sAxisMax,
188 // sAxis1 & sAxisMin, sAxis1 & sAxisMax. 198 // sAxis1 & sAxisMin, sAxis1 & sAxisMax.
189 // boundarytype represents the shape of l 199 // boundarytype represents the shape of locus
190 // from the start point to end point of b 200 // from the start point to end point of boundary.
191 // ex. 201 // ex.
192 // sAxisRho = linear line which start poi 202 // sAxisRho = linear line which start point is fixed at origin.
193 // sAxisPhi = part of circle which center 203 // sAxisPhi = part of circle which center placed at the origin.
194 204
195 void SetCorner(G4int areacode, G4double x, 205 void SetCorner(G4int areacode, G4double x, G4double y, G4double z);
196 206
197 private: 207 private:
198 208
199 virtual void SetBoundaries() = 0; 209 virtual void SetBoundaries() = 0;
200 virtual void SetCorners() = 0; 210 virtual void SetCorners() = 0;
201 211
202 // data members ----------------------------- 212 // data members ---------------------------------------------------------
203 213
204 public: 214 public:
205 215
206 static const G4int sOutside ; 216 static const G4int sOutside ;
207 static const G4int sInside ; 217 static const G4int sInside ;
208 static const G4int sBoundary; 218 static const G4int sBoundary;
209 static const G4int sCorner; 219 static const G4int sCorner;
210 static const G4int sC0Min1Min; 220 static const G4int sC0Min1Min;
211 static const G4int sC0Max1Min; 221 static const G4int sC0Max1Min;
212 static const G4int sC0Max1Max; 222 static const G4int sC0Max1Max;
213 static const G4int sC0Min1Max; 223 static const G4int sC0Min1Max;
214 static const G4int sAxisMin; 224 static const G4int sAxisMin;
215 static const G4int sAxisMax; 225 static const G4int sAxisMax;
216 static const G4int sAxisX; 226 static const G4int sAxisX;
217 static const G4int sAxisY; 227 static const G4int sAxisY;
218 static const G4int sAxisZ; 228 static const G4int sAxisZ;
219 static const G4int sAxisRho; 229 static const G4int sAxisRho;
220 static const G4int sAxisPhi; 230 static const G4int sAxisPhi;
221 static const G4int sAxis0; 231 static const G4int sAxis0;
222 static const G4int sAxis1; 232 static const G4int sAxis1;
223 static const G4int sSizeMask; 233 static const G4int sSizeMask;
224 static const G4int sAxisMask; 234 static const G4int sAxisMask;
225 static const G4int sAreaMask; 235 static const G4int sAreaMask;
226 236
227 protected: 237 protected:
228 238
229 class CurrentStatus 239 class CurrentStatus
230 { 240 {
231 public: << 241 public:
232 242
233 CurrentStatus(); 243 CurrentStatus();
234 virtual ~CurrentStatus(); 244 virtual ~CurrentStatus();
235 245
236 inline G4ThreeVector GetXX(G4int i) 246 inline G4ThreeVector GetXX(G4int i) const { return fXX[i]; }
237 inline G4double GetDistance(G4int i 247 inline G4double GetDistance(G4int i) const { return fDistance[i]; }
238 inline G4int GetAreacode(G4int i 248 inline G4int GetAreacode(G4int i) const { return fAreacode[i]; }
239 inline G4int GetNXX() 249 inline G4int GetNXX() const { return fNXX; }
240 inline G4bool IsDone() 250 inline G4bool IsDone() const { return fDone; }
241 inline G4bool IsValid(G4int i) 251 inline G4bool IsValid(G4int i) const { return fIsValid[i]; }
242 252
243 void SetCurrentStatus(G4int 253 void SetCurrentStatus(G4int i,
244 G4ThreeVector& 254 G4ThreeVector& xx,
245 G4double& 255 G4double& dist,
246 G4int& 256 G4int& areacode,
247 G4bool& 257 G4bool& isvalid,
248 G4int 258 G4int nxx,
249 EValidate 259 EValidate validate,
250 const G4ThreeVector* p, 260 const G4ThreeVector* p,
251 const G4ThreeVector* v = 261 const G4ThreeVector* v = nullptr);
252 262
253 void ResetfDone(EValidate validate, 263 void ResetfDone(EValidate validate,
254 const G4ThreeVector* p, 264 const G4ThreeVector* p,
255 const G4ThreeVector* v = nullp 265 const G4ThreeVector* v = nullptr);
256 266
257 267
258 void DebugPrint() const; 268 void DebugPrint() const;
259 269
260 private: << 270 private:
261 271
262 G4double fDistance[G4VSURFAC 272 G4double fDistance[G4VSURFACENXX];
263 G4ThreeVector fXX[G4VSURFACENXX]; 273 G4ThreeVector fXX[G4VSURFACENXX];
264 G4int fAreacode[G4VSURFAC 274 G4int fAreacode[G4VSURFACENXX];
265 G4bool fIsValid[G4VSURFACE 275 G4bool fIsValid[G4VSURFACENXX];
266 G4int fNXX; 276 G4int fNXX;
267 G4ThreeVector fLastp; 277 G4ThreeVector fLastp;
268 G4ThreeVector fLastv; 278 G4ThreeVector fLastv;
269 EValidate fLastValidate; 279 EValidate fLastValidate;
270 G4bool fDone; 280 G4bool fDone;
271 }; 281 };
272 282
273 class Boundary 283 class Boundary
274 { 284 {
275 public: << 285 public:
276 << 286 Boundary();
277 Boundary() = default; << 287 virtual ~Boundary();
278 virtual ~Boundary() = default; <<
279 288
280 void SetFields(const G4int& areac 289 void SetFields(const G4int& areacode,
281 const G4ThreeVector& d, 290 const G4ThreeVector& d,
282 const G4ThreeVector& x0, 291 const G4ThreeVector& x0,
283 const G4int& boun 292 const G4int& boundarytype);
284 293
285 G4bool IsEmpty() const; 294 G4bool IsEmpty() const;
286 295
287 G4bool GetBoundaryParameters(const G4int 296 G4bool GetBoundaryParameters(const G4int& areacode,
288 G4Thr 297 G4ThreeVector& d,
289 G4Thr 298 G4ThreeVector& x0,
290 G4int 299 G4int& boundarytype) const;
291 300
292 private: << 301 private:
293 << 302 G4int fBoundaryAcode;
294 G4int fBoundaryAcode{-1}; <<
295 G4ThreeVector fBoundaryDirection; 303 G4ThreeVector fBoundaryDirection;
296 G4ThreeVector fBoundaryX0; 304 G4ThreeVector fBoundaryX0;
297 G4int fBoundaryType{0}; << 305 G4int fBoundaryType;
298 }; 306 };
299 307
300 EAxis fAxis[2]; 308 EAxis fAxis[2];
301 G4double fAxisMin[2]; 309 G4double fAxisMin[2];
302 G4double fAxisMax[2]; 310 G4double fAxisMax[2];
303 CurrentStatus fCurStatWithV; 311 CurrentStatus fCurStatWithV;
304 CurrentStatus fCurStat; 312 CurrentStatus fCurStat;
305 G4RotationMatrix fRot; 313 G4RotationMatrix fRot;
306 G4ThreeVector fTrans; 314 G4ThreeVector fTrans;
307 G4int fHandedness; 315 G4int fHandedness;
308 class G4SurfCurNormal 316 class G4SurfCurNormal
309 { 317 {
310 public: 318 public:
311 <<
312 G4ThreeVector p; 319 G4ThreeVector p;
313 G4ThreeVector normal; 320 G4ThreeVector normal;
314 }; 321 };
315 G4SurfCurNormal fCurrentNormal; 322 G4SurfCurNormal fCurrentNormal;
316 G4bool fIsValidNorm; 323 G4bool fIsValidNorm;
317 G4double kCarTolerance; 324 G4double kCarTolerance;
318 325
319 private: 326 private:
320 327
321 G4VTwistSurface* fNeighbours[4]; // {0,1,2, 328 G4VTwistSurface* fNeighbours[4]; // {0,1,2,3} = sAxis0min, sAxis1min,
322 // 329 // sAxis0max, sAxis1max
323 G4ThreeVector fCorners[4]; // corners o 330 G4ThreeVector fCorners[4]; // corners of the surface in local coordinate
324 Boundary fBoundaries[4]; // boundarie 331 Boundary fBoundaries[4]; // boundaries of the surface.
325 G4String fName; 332 G4String fName;
326 333
327 class G4SurfSideQuery 334 class G4SurfSideQuery
328 { 335 {
329 public: 336 public:
330 <<
331 G4ThreeVector me; 337 G4ThreeVector me;
332 G4ThreeVector vec; 338 G4ThreeVector vec;
333 G4bool withTol; 339 G4bool withTol;
334 G4int amIOnLeftSide; 340 G4int amIOnLeftSide;
335 }; 341 };
336 G4SurfSideQuery fAmIOnLeftSide; 342 G4SurfSideQuery fAmIOnLeftSide;
337 }; 343 };
338 344
339 //============================================ 345 //========================================================
340 // inline functions 346 // inline functions
341 //============================================ 347 //========================================================
342 348
343 struct Intersection 349 struct Intersection
344 { 350 {
345 G4double phi ; // parameter phi 351 G4double phi ; // parameter phi
346 G4double u ; // parameter u 352 G4double u ; // parameter u
347 G4ThreeVector xx ; // intersection point i 353 G4ThreeVector xx ; // intersection point in cartesian
348 G4double distance ; // distance to intersec 354 G4double distance ; // distance to intersection
349 G4int areacode ; // the areacode of the 355 G4int areacode ; // the areacode of the intersection
350 G4bool isvalid ; // valid intersection ? 356 G4bool isvalid ; // valid intersection ??
351 357
352 }; 358 };
353 359
354 inline 360 inline
355 G4bool DistanceSort( const Intersection& a, co 361 G4bool DistanceSort( const Intersection& a, const Intersection& b)
356 { 362 {
357 return a.distance < b.distance ; 363 return a.distance < b.distance ;
358 } 364 }
359 365
360 inline 366 inline
361 G4bool EqualIntersection( const Intersection& 367 G4bool EqualIntersection( const Intersection& a, const Intersection& b)
362 { 368 {
363 return ( ( a.xx - b.xx ).mag() < 1E-9*CLHEP: 369 return ( ( a.xx - b.xx ).mag() < 1E-9*CLHEP::mm ) ;
364 } 370 }
365 371
366 #include "G4VTwistSurface.icc" 372 #include "G4VTwistSurface.icc"
367 373
368 #endif 374 #endif
369 375