Geant4 Cross Reference |
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 // G4tgrEvaluator implementation 27 // 28 // Author: P.Arce, CIEMAT (November 2007) 29 // -------------------------------------------------------------------- 30 31 #include "G4tgrEvaluator.hh" 32 33 #include <cmath> 34 35 // -------------------------------------------------------------------- 36 G4tgrEvaluator::G4tgrEvaluator() 37 { 38 AddCommonFunctions(); 39 } 40 41 // -------------------------------------------------------------------- 42 G4tgrEvaluator::~G4tgrEvaluator() 43 { 44 } 45 46 // -------------------------------------------------------------------- 47 void G4tgrEvaluator::print_error(G4int estatus) const 48 { 49 switch(estatus) 50 { 51 case ERROR_SYNTAX_ERROR: 52 G4cerr << "G4tgrEvaluator: syntax error!" << G4endl; 53 return; 54 default: 55 G4Evaluator::print_error(); 56 return; 57 } 58 } 59 60 G4double fltsin(G4double arg) { return std::sin(arg); } 61 G4double fltcos(G4double arg) { return std::cos(arg); } 62 G4double flttan(G4double arg) { return std::tan(arg); } 63 G4double fltasin(G4double arg) { return std::asin(arg); } 64 G4double fltacos(G4double arg) { return std::acos(arg); } 65 G4double fltatan(G4double arg) { return std::atan(arg); } 66 G4double fltatan2(G4double arg1, G4double arg2) { return std::atan2(arg1, arg2); } 67 G4double fltsinh(G4double arg) { return std::sinh(arg); } 68 G4double fltcosh(G4double arg) { return std::cosh(arg); } 69 G4double flttanh(G4double arg) { return std::tanh(arg); } 70 // G4double fltasinh( G4double arg ){ return std::asinh(arg); } 71 // G4double fltacosh( G4double arg ){ return std::acosh(arg); } 72 // G4double fltatanh( G4double arg ){ return std::atanh(arg); } 73 G4double fltsqrt(G4double arg) { return std::sqrt(arg); } 74 G4double fltexp(G4double arg) { return std::exp(arg); } 75 G4double fltlog(G4double arg) { return std::log(arg); } 76 G4double fltlog10(G4double arg) { return std::log10(arg); } 77 G4double fltpow(G4double arg1, G4double arg2) { return std::pow(arg1, arg2); } 78 79 // -------------------------------------------------------------------- 80 void G4tgrEvaluator::AddCommonFunctions() 81 { 82 setFunction("sin", (*fltsin)); 83 setFunction("cos", (*fltcos)); 84 setFunction("tan", (*flttan)); 85 setFunction("asin", (*fltasin)); 86 setFunction("acos", (*fltacos)); 87 setFunction("atan", (*fltatan)); 88 setFunction("atan2", (*fltatan2)); 89 setFunction("sinh", (*fltsinh)); 90 setFunction("cosh", (*fltcosh)); 91 setFunction("tanh", (*flttanh)); 92 // setFunction("asinh", (*fltasinh)); 93 // setFunction("acosh", (*fltacosh)); 94 // setFunction("atanh", (*fltatanh)); 95 setFunction("sqrt", (*fltsqrt)); 96 setFunction("exp", (*fltexp)); 97 setFunction("log", (*fltlog)); 98 setFunction("log10", (*fltlog10)); 99 setFunction("pow", (*fltpow)); 100 } 101