OPAL (Object Oriented Parallel Accelerator Library) 2022.1
OPAL
CoordinateTransform.h
Go to the documentation of this file.
1/*
2 * Copyright (c) 2018, Martin Duy Tat
3 * All rights reserved.
4 * Redistribution and use in source and binary forms, with or without
5 * modification, are permitted provided that the following conditions are met:
6 * 1. Redistributions of source code must retain the above copyright notice,
7 * this list of conditions and the following disclaimer.
8 * 2. Redistributions in binary form must reproduce the above copyright notice,
9 * this list of conditions and the following disclaimer in the documentation
10 * and/or other materials provided with the distribution.
11 * 3. Neither the name of STFC nor the names of its contributors may be used to
12 * endorse or promote products derived from this software without specific
13 * prior written permission.
14 *
15 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
16 * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
17 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
18 * ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
19 * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
20 * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
21 * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
22 * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
23 * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
24 * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
25 * POSSIBILITY OF SUCH DAMAGE.
26 */
27
28#ifndef COORDINATE_TRANSFORM_H
29#define COORDINATE_TRANSFORM_H
30
48#include <vector>
49
50namespace coordinatetransform {
51
52struct myParams {
53 double s_0;
54 double lambdaleft;
56 double rho;
57};
58
60public:
73 CoordinateTransform(const double &xlab,
74 const double &ylab,
75 const double &zlab,
76 const double &s_0,
77 const double &lambdaleft,
78 const double &lambdaright,
79 const double &rho);
87 std::vector<double> getTransformation() const;
91 std::vector<double> calcReferenceTrajectory(const double &s) const;
95 std::vector<double> getUnitTangentVector(const double &s) const;
96private:
101 void calcSCoordinate(const double &xlab, const double &ylab);
106 void calcXCoordinate(const double &xlab, const double &ylab);
114 void transformFromEntranceCoordinates(std::vector<double> &coordinates,
115 const double &boundingBoxLength);
116 double s_0_m;
119 double rho_m;
120 double x_m;
121 double z_m;
122 double s_m;
123 static const double error;
124 static const int workspaceSize;
125 static const int algorithm;
126};
127
129double getUnitTangentVectorX(double s, void *p);
131double getUnitTangentVectorY(double s, void *p);
132
133}
134#endif
double getUnitTangentVectorX(double s, void *p)
double getUnitTangentVectorY(double s, void *p)
std::vector< double > getUnitTangentVector(const double &s) const
void transformFromEntranceCoordinates(std::vector< double > &coordinates, const double &boundingBoxLength)
void calcXCoordinate(const double &xlab, const double &ylab)
std::vector< double > getTransformation() const
std::vector< double > calcReferenceTrajectory(const double &s) const
CoordinateTransform & operator=(const CoordinateTransform &transform)
void calcSCoordinate(const double &xlab, const double &ylab)