BDSIM
BDSIM is a Geant4 extension toolkit for simulation of particle transport in accelerator beamlines.
Loading...
Searching...
No Matches
BDSInterpolator4DLinearMag.cc
1/*
2Beam Delivery Simulation (BDSIM) Copyright (C) Royal Holloway,
3University of London 2001 - 2023.
4
5This file is part of BDSIM.
6
7BDSIM is free software: you can redistribute it and/or modify
8it under the terms of the GNU General Public License as published
9by the Free Software Foundation version 3 of the License.
10
11BDSIM is distributed in the hope that it will be useful, but
12WITHOUT ANY WARRANTY; without even the implied warranty of
13MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14GNU General Public License for more details.
15
16You should have received a copy of the GNU General Public License
17along with BDSIM. If not, see <http://www.gnu.org/licenses/>.
18*/
19#include "BDSArray4DCoords.hh"
20#include "BDSFieldValue.hh"
21#include "BDSInterpolator4DLinearMag.hh"
22#include "BDSInterpolatorRoutines.hh"
23
24#include "G4Types.hh"
25
27 BDSInterpolator4D(arrayIn)
28{;}
29
30BDSInterpolator4DLinearMag::~BDSInterpolator4DLinearMag()
31{;}
32
34 G4double y,
35 G4double z,
36 G4double t) const
37{
38 BDSFieldValue localData[2][2][2][2];
39 G4double xFrac, yFrac, zFrac, tFrac;
40 array->ExtractSection2x2x2x2(x, y, z, t, localData, xFrac, yFrac, zFrac, tFrac);
41 return BDS::Linear4DMag(localData, xFrac, yFrac, zFrac, tFrac);
42}
Overlay of 4D array that provides uniform only spatial coordinate mapping.
virtual void ExtractSection2x2x2x2(G4double x, G4double y, G4double z, G4double t, BDSFieldValue(&localData)[2][2][2][2], G4double &xFrac, G4double &yFrac, G4double &zFrac, G4double &tFrac) const
Extract 2x2x2x2 points lying around coordinate x.
virtual BDSFieldValue GetInterpolatedValueT(G4double x, G4double y, G4double z, G4double t) const
BDSInterpolator4DLinearMag()=delete
Private default constructor to force use of provided one.
Interface for all 4D interpolators.
BDSArray4DCoords * array
The field data.
T Linear4DMag(const T p[2][2][2][2], G4double x, G4double y, G4double z, G4double t)
Linear interpolation in 4 dimensions including magnitude interpolation.