include/Geant4/G4VFSALIntegrationStepper.icc

0001 //
0002 // ********************************************************************
0003 // * License and Disclaimer                                           *
0004 // *                                                                  *
0005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
0006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
0007 // * conditions of the Geant4 Software License,  included in the file *
0008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
0009 // * include a list of copyright holders.                             *
0010 // *                                                                  *
0011 // * Neither the authors of this software system, nor their employing *
0012 // * institutes,nor the agencies providing financial support for this *
0013 // * work  make  any representation or  warranty, express or implied, *
0014 // * regarding  this  software system or assume any liability for its *
0015 // * use.  Please see the license in the file  LICENSE  and URL above *
0016 // * for the full disclaimer and the limitation of liability.         *
0017 // *                                                                  *
0018 // * This  code  implementation is the result of  the  scientific and *
0019 // * technical work of the GEANT4 collaboration.                      *
0020 // * By using,  copying,  modifying or  distributing the software (or *
0021 // * any work based  on the software)  you  agree  to acknowledge its *
0022 // * use  in  resulting  scientific  publications,  and indicate your *
0023 // * acceptance of all terms of the Geant4 Software license.          *
0024 // ********************************************************************
0025 //
0026 // G4VFSALIntegrationStepper inline methods implementation
0027 //
0028 // Author: Somnath Banerjee, Google Summer of Code 2015
0029 // Supervision: John Apostolakis, CERN
0030 // --------------------------------------------------------------------
0031
0032 inline
0033 G4EquationOfMotion* G4VFSALIntegrationStepper::GetEquationOfMotion()
0034 {
0035   return fEquation_Rhs;
0036 }
0037
0038 inline void
0039  G4VFSALIntegrationStepper::SetEquationOfMotion(G4EquationOfMotion* newEquation)
0040 {
0041   if( newEquation != nullptr )
0042   {
0043     fEquation_Rhs = newEquation;
0044   }
0045 }
0046
0047 inline
0048 G4int G4VFSALIntegrationStepper::GetNumberOfVariables() const
0049 {
0050   return fNoIntegrationVariables;
0051 }
0052
0053 inline
0054 G4int G4VFSALIntegrationStepper::GetNumberOfStateVariables() const
0055 {
0056   return fNoStateVariables;
0057 }
0058
0059 inline
0060 void G4VFSALIntegrationStepper::NormaliseTangentVector( G4double vec[6] )
0061 {
0062   G4double drds2 = vec[3]*vec[3]+vec[4]*vec[4]+vec[5]*vec[5];
0063
0064   if( std::fabs(drds2 - 1.0) > 1.e-14 )
0065   {
0066     G4double normx = 1.0 / std::sqrt(drds2);
0067     for(auto i=3; i<6; ++i) { vec[i] *= normx; }
0068   }
0069 }
0070
0071 inline
0072 void G4VFSALIntegrationStepper::NormalisePolarizationVector( G4double vec[12] )
0073 {
0074   G4double drds2 = vec[9]*vec[9]+vec[10]*vec[10]+vec[11]*vec[11];
0075
0076   if( drds2 > 0. )
0077   {
0078     if( std::fabs(drds2 - 1.0) > 1.e-14 )
0079     {
0080       G4double normx = 1.0 / std::sqrt(drds2);
0081       for(auto i=9; i<12; ++i)  { vec[i] *= normx; }
0082     }
0083   }
0084 }