Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer << 3 // * DISCLAIMER * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th << 5 // * The following disclaimer summarizes all the specific disclaimers * 6 // * the Geant4 Collaboration. It is provided << 6 // * of contributors to this software. The specific disclaimers,which * 7 // * conditions of the Geant4 Software License << 7 // * govern, are listed with their locations in: * 8 // * LICENSE and available at http://cern.ch/ << 8 // * http://cern.ch/geant4/license * 9 // * include a list of copyright holders. << 10 // * 9 // * * 11 // * Neither the authors of this software syst 10 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 11 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 12 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 13 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file << 14 // * use. * 16 // * for the full disclaimer and the limitatio << 17 // * 15 // * * 18 // * This code implementation is the result << 16 // * This code implementation is the intellectual property of the * 19 // * technical work of the GEANT4 collaboratio << 17 // * GEANT4 collaboration. * 20 // * By using, copying, modifying or distri << 18 // * By copying, distributing or modifying the Program (or any work * 21 // * any work based on the software) you ag << 19 // * based on the Program) you indicate your acceptance of this * 22 // * use in resulting scientific publicati << 20 // * statement, and all its terms. * 23 // * acceptance of all terms of the Geant4 Sof << 24 // ******************************************* 21 // ******************************************************************** 25 // 22 // 26 // G4SimpleHeum implementation << 27 // 23 // 28 // Created: W.Wander <wwc@mit.edu>, 12/09/1997 << 24 // $Id: G4SimpleHeum.cc,v 1.5 2001/07/11 09:59:13 gunter Exp $ 29 // ------------------------------------------- << 25 // GEANT4 tag $Name: geant4-05-02 $ >> 26 // >> 27 // Simple Heum: >> 28 // x_1 = x_0 + h * >> 29 // 1/4 * dx(t0,x0) + >> 30 // 3/4 * dx(t0+2/3*h, x0+2/3*h*(dx(t0+h/3,x0+h/3*dx(t0,x0)))) >> 31 // >> 32 // third order solver >> 33 // >> 34 // W.Wander <wwc@mit.edu> 12/09/97 >> 35 // 6.11.98. V.Grichine new data member fNumberOfVariables >> 36 30 37 31 #include "G4SimpleHeum.hh" 38 #include "G4SimpleHeum.hh" 32 #include "G4ThreeVector.hh" 39 #include "G4ThreeVector.hh" 33 40 34 ////////////////////////////////////////////// 41 /////////////////////////////////////////////////////////////////////////// 35 // 42 // 36 // Constructor 43 // Constructor 37 // << 44 38 G4SimpleHeum::G4SimpleHeum(G4EquationOfMotion* << 45 G4SimpleHeum::G4SimpleHeum(G4Mag_EqRhs *EqRhs, G4int num_variables): 39 : G4MagErrorStepper(EqRhs, num_variables), << 46 G4MagErrorStepper(EqRhs, num_variables), 40 fNumberOfVariables(num_variables) << 47 fNumberOfVariables(num_variables) 41 { 48 { 42 dydxTemp = new G4double[fNumberOfVariables] 49 dydxTemp = new G4double[fNumberOfVariables] ; 43 dydxTemp2 = new G4double[fNumberOfVariables] 50 dydxTemp2 = new G4double[fNumberOfVariables] ; 44 yTemp = new G4double[fNumberOfVariables] 51 yTemp = new G4double[fNumberOfVariables] ; 45 yTemp2 = new G4double[fNumberOfVariables] 52 yTemp2 = new G4double[fNumberOfVariables] ; 46 } 53 } 47 54 >> 55 48 ////////////////////////////////////////////// 56 ////////////////////////////////////////////////////////////////////////// 49 // 57 // 50 // Destructor 58 // Destructor 51 // << 59 52 G4SimpleHeum::~G4SimpleHeum() 60 G4SimpleHeum::~G4SimpleHeum() 53 { 61 { 54 delete [] dydxTemp; << 62 delete[] dydxTemp; 55 delete [] dydxTemp2; << 63 delete[] dydxTemp2; 56 delete [] yTemp; << 64 delete[] yTemp; 57 delete [] yTemp2; << 65 delete[] yTemp2; 58 } 66 } 59 67 >> 68 60 ////////////////////////////////////////////// 69 ////////////////////////////////////////////////////////////////////// 61 // 70 // 62 // DumbStepper << 63 // 71 // >> 72 64 void 73 void 65 G4SimpleHeum::DumbStepper( const G4double yIn[ << 74 G4SimpleHeum::DumbStepper( const G4double yIn[], 66 const G4double dydx << 75 const G4double dydx[], 67 G4double h, << 76 G4double h, 68 G4double yOut << 77 G4double yOut[]) 69 { 78 { 70 for( G4int i = 0; i < fNumberOfVariables; ++ << 79 // const G4int nvar = 6 ; >> 80 >> 81 G4int i; >> 82 >> 83 for( i = 0; i < fNumberOfVariables; i++ ) 71 { 84 { 72 yTemp[i] = yIn[i] + (1.0/3.0) * h * dydx[ 85 yTemp[i] = yIn[i] + (1.0/3.0) * h * dydx[i] ; 73 } 86 } 74 87 75 RightHandSide(yTemp,dydxTemp); 88 RightHandSide(yTemp,dydxTemp); 76 89 77 for( G4int i = 0; i < fNumberOfVariables; ++ << 90 for( i = 0; i < fNumberOfVariables; i++ ) 78 { 91 { 79 yTemp2[i] = yIn[i] + (2.0/3.0) * h * dydxT 92 yTemp2[i] = yIn[i] + (2.0/3.0) * h * dydxTemp[i] ; 80 } 93 } 81 94 82 RightHandSide(yTemp2,dydxTemp2); 95 RightHandSide(yTemp2,dydxTemp2); 83 96 84 for( G4int i = 0; i < fNumberOfVariables; ++ << 97 for( i = 0; i < fNumberOfVariables; i++ ) 85 { 98 { 86 yOut[i] = yIn[i] + h * (0.25 * dydx[i] + 0 << 99 yOut[i] = yIn[i] + h * ( 0.25 * dydx[i] + >> 100 0.75 * dydxTemp2[i]); 87 } 101 } 88 102 89 if ( fNumberOfVariables == 12 ) { NormaliseP << 103 // NormaliseTangentVector( yOut ); >> 104 >> 105 return ; 90 } 106 } 91 107