ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4BogackiShampine23.hh
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4BogackiShampine23.hh
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 // G4BogackiShampine23
27 //
28 // Class description:
29 //
30 // Bogacki-Shampine - 4 - 3(2) non-FSAL implementation
31 //
32 // An implementation of the embedded RK method from the paper
33 // [1] P. Bogacki and L. F. Shampine,
34 // "A 3(2) pair of Runge - Kutta formulas"
35 // Appl. Math. Lett., vol. 2, no. 4, pp. 321-325, Jan. 1989.
36 //
37 // This version does not utilise the FSAL property of the method,
38 // which would allow the reuse of the last derivative in the next step.
39 // (Alternative FSAL implementation created with revised interface)
40 
41 // Created: Somnath Banerjee, Google Summer of Code 2015, 20 May 2015
42 // Supervision: John Apostolakis, CERN
43 // --------------------------------------------------------------------
44 #ifndef G4BOGACKI_SHAMPINE23_HH
45 #define G4BOGACKI_SHAMPINE23_HH
46 
48 #include "G4FieldTrack.hh"
49 
51 {
52  public:
53 
55  G4int numberOfVariables = 6);
56 
57  virtual void Stepper(const G4double yInput[],
58  const G4double dydx[],
59  G4double hstep,
60  G4double yOutput[],
61  G4double yError[]) override;
62 
63  void Stepper(const G4double yInput[],
64  const G4double dydx[],
65  G4double hstep,
66  G4double yOutput[],
67  G4double yError[],
68  G4double dydxOutput[]);
69 
72 
73  virtual G4double DistChord() const override;
74  virtual G4int IntegratorOrder() const override { return 3; }
75 
76  private:
77 
78  void makeStep(const G4double yInput[],
79  const G4double dydx[],
80  const G4double hstep,
81  G4double yOutput[],
82  G4double* dydxOutput = nullptr,
83  G4double* yError = nullptr) const;
84 
89  G4double fhstep = -1.0;
90 };
91 
92 #endif