iniconstant.F

Go to the documentation of this file.

Copyright>        OpenRadioss
Copyright>        Copyright (C) 1986-2025 Altair Engineering Inc.
Copyright>
Copyright>        This program is free software: you can redistribute it and/or modify
Copyright>        it under the terms of the GNU Affero General Public License as published by
Copyright>        the Free Software Foundation, either version 3 of the License, or
Copyright>        (at your option) any later version.
Copyright>
Copyright>        This program is distributed in the hope that it will be useful,
Copyright>        but WITHOUT ANY WARRANTY; without even the implied warranty of
Copyright>        MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
Copyright>        GNU Affero General Public License for more details.
Copyright>
Copyright>        You should have received a copy of the GNU Affero General Public License
Copyright>        along with this program.  If not, see <https://www.gnu.org/licenses/>.
Copyright>
Copyright>
Copyright>        Commercial Alternative: Altair Radioss Software
Copyright>
Copyright>        As an alternative to this open-source version, Altair also offers Altair Radioss
Copyright>        software under a commercial license.  Contact Altair to discuss further if the
Copyright>        commercial version may interest you: https://www.altair.com/radioss/.
!||====================================================================
!||    iniconstant   ../engine/source/system/iniconstant.F
!||--- called by ------------------------------------------------------
!||    radioss2      ../engine/source/engine/radioss2.F
!||====================================================================
      SUBROUTINE iniconstant
#include "implicit_f.inc"   
C--------------------------------------
C      INTEGERS
C--------------------------------------
       two          = 2.
       three         = 3.
       four        = 4.
       five          = 5.
       six           = 6.
       seven          = 7.
       eight          = 8.
       nine          = 9.
       ten           = 10.
       eleven          = 11.
       twelve         = 12.
       fourteen      = 14.
       fifteen        = 15.
       sixteen         = 16.
       eighteen         = 18.
       twenty         = 20.
       twenty4        = 24.
       twenty5        = 25.
       twenty7        = 27.
       thirty        = 30.
       thirty2       = 32.
       thirty6       = 36.
       fourty      = 40.
       fourty5     = 45.
       fourty8     = 48.
       fifty     = 50.
       cinquante4     = 54.
       sixty4     = 64.
       eighty   = 80.
       eighty16 = 96.
       eighty19 = 99.
       hundred          = 100.
       hundred28        = 128.
       hundred50        = 150.
       hundred80        = 180.
       three100      = 300.
       twohundred25    = 225.
       twohundred95    = 295.
       twohundred98    = 298.
       twohundred93    = 293.
       threehundred24   = 324.
       fivehundred12    = 512.
       two1000      = 2000. 
       three1000     = 3000.  
C--------------------------------------
C      10^n
C--------------------------------------
       ep02  = hundred     
       ep03  = ep02 * ten    
       ep04  = ep03 * ten   
       ep05  = ep04 * ten    
       ep06  = ep05 * ten     
       ep07  = ep06 * ten              
       ep08  = ep07 * ten              
       ep09  = ep08 * ten              
       ep10  = ep09 * ten              
       ep11  = ep10 * ten           
       ep12  = ep11 * ten           
       ep13  = ep12 * ten           
       ep14  = ep13 * ten           
       ep15  = ep14 * ten           
       ep16  = ep15 * ten             
       ep17  = ep16 * ten             
       ep18  = ep17 * ten             
       ep19  = ep18 * ten             
       ep20  = ep19 * ten             
       ep29  = ep20 * ep09           
       ep30  = ep20 * ep10            
       ep31  = ep30 * ten        
       ep35  = ep30 * ep05             
       ep36  = ep30 * ep06             
       ep37  = ep30 * ep07               
C--------------------------------------
C      10^-n
C--------------------------------------
       em01 = one / ten  
       em02 = one / ep02 
       em03 = one / ep03 
       em04 = one / ep04 
       em05 = one / ep05 
       em06 = one / ep06 
       em07 = one / ep07 
       em08 = one / ep08 
       em09 = one / ep09 
       em10 = one / ep10 
       em11 = one / ep11 
       em12 = one / ep12 
       em13 = one / ep13 
       em14 = one / ep14 
       em15 = one / ep15 
       em16 = one / ep16 
       em17 = one / ep17 
       em18 = one / ep18 
       em19 = one / ep19 
       em20 = one / ep20 
       em30 = one / ep30
       em35 = one / ep35 
C--------------------------------------
C      m10^-n
C--------------------------------------
       zep1      = em01 
       zep2      = two   / ten 
       zep3      = three  / ten 
       zep4      = four / ten 
       zep5      = five   / ten 
       zep6      = one - zep4 
       zep7      = one - zep3 
       zep8      = one - zep2 
       zep9      = one - zep1 
C
       zep01     =        em02
       zep02     = two   / ep02
       zep03     = three  / ep02
       zep04     = four / ep02
       zep05     = five   / ep02
       zep06     = zep1 - zep04
       zep07     = zep1 - zep03
       zep08     = zep1 - zep02
       zep09     = zep1 - zep01   
C
       twoem3   = two   / ep03
       threeem3  = three  / ep03
       fourem3 = four / ep03
       fiveem3   = five   / ep03
       sixem3    = zep01 - fourem3
       sevenem3   = zep01 - threeem3 
       eightem3   = zep01 - twoem3
       nineem3   = zep01 - em03      
C
       twoem4   = two   / ep04
       threeem4  = three  / ep04
       fourem4 = four / ep04
       fiveem4   = five   / ep04
       sixem4    = em03 - fourem4 
       sevenem4   = em03 - threeem4 
       eightem4   = em03 - twoem4
       nineem4   = em03 - em04             
C
       twoem5   = two   / ep05
       threeem5  = three  / ep05
       fourem5 = four / ep05
       fiveem5   = five   / ep05
       sixem5    = em04 - fourem5
       sevenem5   = em04 - threeem5 
       eightem5   = em04 - twoem5
       nineem5   = em04 - em05             
C
       twoem6   = two   / ep06
       threeem6  = three  / ep06
       fourem6 = four / ep06
       fiveem6   = five   / ep06
       sixem6    = em05 - fourem6
       sevenem6   = em05 - threeem6
       eightem6   = em05 - twoem6
       nineem6   = em05 - em06        
C--------------------------------------
C      m10^n
C--------------------------------------
       twoep30   = two  *ep30
       threeep30  = three *ep30
       fourep30 = four*ep30  
       zep9ep30   = nine  *ep29
       onep1ep30   = eleven  * ep29
       onep2ep30   = twelve * ep29 
C--------------------------------------
C      Fractions
C--------------------------------------
       half    = one/two         
       third     = one/three        
       fourth     = one/four       
       one_fifth    = one/five         
       one_over_6     = one/six          
       one_over_8    = one/eight         
       one_over_9    = one/nine         
       one_over_12   = one/twelve        
       one_over_16   = one/sixteen        
       one_over_64     = one/sixty4    
       one_over_512    = one/fivehundred12 
       two_third    = two/three      
       three_half     = three/two      
       three_over_4    = three/four    
       trhee_over_14 = three/fourteen  
       tseize    = three/sixteen
       four_over_3    = four/three    
       four_over_5     = four/five     
       qneuf     = four/nine
       five_over_3    = five/three      
       five_over_4   = five/four     
       five_over_6      = five/six        
       six_over_5    = six/five        
       seven_half    = seven/two       
       seven_over_5     = seven/five       
       seven_over_9    = seven/nine       
       nine_over_2  = nine/two       
       nseize    = nine/sixteen
       nine_over_20    = nine/twenty      
C--------------------------------------
C      0.xxxx
C--------------------------------------
       zep11     = zep1 + zep01
       zep12     = zep1 + zep02
       zep13     = zep1 + zep03
       zep14     = zep1 + zep04
       zep15     = zep1 + zep05
       zep16     = zep1 + zep06
       zep17     = zep1 + zep07
       zep18     = zep1 + zep08
       zep19     = zep1 + zep09
C
       zep21     = zep2  + zep01
       zep23     = zep2  + zep03
       zep25     = zep2  + zep05
       zep26     = zep2  + zep06
       zep213    = zep21 + threeem3
       zep29     = zep2  + zep09
C
       zep33     = zep3  + zep03
       zep34     = zep3  + zep04
       zep35     = zep3  + zep05
       zep37     = zep3  + zep07
       zep38     = zep3  + zep08
       zep333    = zep33 + threeem3
       zep373    = zep37 + threeem3         
C
       zep41     = zep4  + zep01
       zep43     = zep4  + zep03
       zep45     = zep4  + zep05
       zep46     = zep4  + zep06
       zep47     = zep4  + zep07
       zep499    = zep4  + zep09 + nineem3
       zep4187   = zep41 + eightem3 + sevenem4
C
       zep62     = one   - zep38
       zep65     = one   - zep35
       zep66     = one   - zep34 
       zep67     = one   - zep33
c
       zep71     = one   - zep29
       zep74     = one   - zep26         
       zep707    = one   - zep3 + sevenem3
C
       zep81     = one   - zep19
       zep82     = one   - zep18
       zep84     = one   - zep16         
       zep85     = one   - zep15
       zep87     = one   - zep13
       zep89     = one   - zep11
C
       zep92     = one   - zep08
       zep93     = one   - zep07       
       zep94     = one   - zep06
       zep95     = one   - zep05 
       zep96     = one   - zep04
       zep97     = one   - zep03
       zep99     = one   - zep01
       zep995    = one   - fiveem3
       zep997    = one   - threeem3
       zep999    = one   - em03
       zep9999   = one   - em04
       zep99999  = one   - em05
       zep931    = zep93  + em03        
c
       zep015    = zep15 / ten
       zep015625 = one / sixty4
       zep025    = one / fourty
       zep037    = zep03 + sevenem3                                          
       zep072169 = zep07 + twoem3 + em04 + sixem5 + nineem6
       zep075    = three / fourty
C
       onep01    = one + zep01
       onep001   = one + em03
       onep0001  = one + em04        
       onep00001 = one + em05
       onep02  = one + zep02
       onep03  = one + zep03
       onep04  = one + zep04
       onep05  = one + zep05
       onep06  = one + zep06
       onep08  = one + zep08
       onep09  = one + zep09
C
       onep1   = one    + zep1
       onep11  = onep1  + zep01 
       onep12  = onep1  + zep02
       onep13  = onep1  + zep03
       onep14  = onep1  + zep04
       onep15  = onep1  + zep05 
       onep16  = onep1  + zep06
       onep17  = onep1  + zep07
       onep18  = onep1  + zep08
       onep19  = onep1  + zep09
       onep136 = onep13  + sixem3
       onep138 = onep13  + eightem3
C
       onep2   = one + zep2 
       onep21  = onep2  + zep01
       onep23  = onep2  + zep03
       onep24  = onep2  + zep04
       onep25  = five / four      
       onep26  = onep2  + zep06
       onep27  = onep2  + zep07
       onep28  = onep2  + zep08        
       onep251 = onep25  + em03
       onep252 = onep25  + twoem3
C
       onep3   = one + zep3
       onep31  = onep3  + zep01  
       onep32  = onep3  + zep02
       onep33  = onep3  + zep03
       onep34  = onep3  + zep04
       onep35  = onep3  + zep05
       onep36  = onep3  + zep06
       onep37  = onep3  + zep07
       onep38  = onep3  + zep08
       onep39  = onep3  + zep09
       onep333 = onep33  + threeem3
       onep3348= onep33  + fourem3 + eightem4      
       onep369 = onep36  + nineem3
C
       onep4   = one + zep4
       onep41  = onep4  + zep01
       onep42  = onep4  + zep02
       onep44  = onep4  + zep04         
       onep45  = onep4  + zep05
       onep46  = onep4  + zep06
       onep48  = onep4  + zep08
       onep49  = onep4  + zep09
       onep414 = onep41 + fourem3
       onep484 = onep48 + fourem3
C
       onep5   = three_half
       onep51  = onep5 + zep01
       onep52  = onep5 + zep02
       onep57  = onep5 + zep07
       onep59  = onep5 + zep09
C
       onep6   = one + zep6
       onep62  = onep6 + zep02
       onep65  = onep6 + zep05
       onep66  = onep6 + zep06
       onep67  = onep6 + zep07
       onep68  = onep6 + zep08 
C
       onep7   = one + zep7
       onep72  = onep7 + zep02
       onep73  = onep7 + zep03        
       onep75  = onep7 + zep05                
       onep78  = onep7 + zep08
       onep732 = onep73 + twoem3
C
       onep8   = one + zep8
       onep82  = onep8 + zep02
       onep84  = onep8 + zep04     
       onep88  = onep8 + zep08
       onep833  = onep8 + zep03 + threeem3
       onep8333 = one   + five_over_6
       onep8688 = onep8 + zep06 + eightem3 + eightem4
C
       onep9   = one + zep9
       onep92  = onep9 + zep02
       onep96  = onep9 + zep06
       onep98  = onep9 + zep08
CC  
       twop03 = two + zep03
       twop05 = two + zep05
       twop06 = two + zep06        
       twop07 = two + zep07
       twop08 = two + zep08
       twop1  = two + zep1    
       twop2  = two + zep2
       twop30 = two + zep3
       twop4  = two + zep4
       twop5  = five / two
       twop60 = two + zep6
       twop8  = two + zep8
C
       twop11 = twop1 + zep01
       twop12 = twop1 + zep02
       twop13 = twop1 + zep03
       twop16 = twop1 + zep06    
       twop19 = twop1 + zep09
       twop23 = twop2 + zep03
       twop24 = twop2 + zep04
       twop25 = twop2 + zep05
       twop33 = twop30 + zep03
       twop38 = twop30 + zep08
       twop45 = twop4 + zep05
       twop47 = twop4 + zep07
       twop52 = twop5 + zep02
       twop83 = twop8 + zep03        
       twop84 = twop8 + zep04        
       twop85 = twop8 + zep05
C
       twop185 = twop1 + zep08 + fiveem3  
       twop261 = twop2 + zep06 + em03    
       twop645 = twop60 + zep04 + fiveem3
C
       threep02  = three + zep02
       threep06  = three + zep06
       threep14  = three + zep14
       threep15  = three + zep15
       threep25  = three + zep2 + zep05
       threep27  = three + zep2 + zep07
       threep44  = three + zep4 + zep04
       threep46  = three + zep4 + zep06
       threep52  = three + zep5 + zep02
       threep72  = three + zep7 + zep02
       threep76  = three + zep7 + zep06
       threep8   = three + zep8
       threep94  = three + zep9 + zep04
       threep266 = three + zep2 + zep06 + sixem3
       threep464 = three + zep4 + zep06 + fourem3
       threep4641= threep464 + em04
C
       fourp186 = four + zep1 + zep08 + sixem3
       fourp40  = four + zep4
       fourp76  = four + zep7 + zep06
C
       fivep5  = five + zep5
       fivep56 = five + zep5 + zep06
       fivep76 = five + zep7 + zep06
C
       sip2832 = six + zep2 + zep08 + threeem3 + twoem4 
C       
       ninep24  = nine + zep2 + zep04
       ninep38  = nine + zep3 + zep08
       ninep793 = nine + zep7 + zep09 + threeem3
 
       elevenp5     = eleven     + zep5
       onzep225   = eleven  + zep2 + zep02+fiveem3
       twelvep6    = twelve    + zep6
       fourteenp2 = fourteen + zep2
       sixteenp5    = sixteen    + zep5 
       o88p9844   = 88. + zep9 + zep08 + fourem3 + fourem4
C
       onep3em5     = onep3 * em05
       onep001em32  = 1001. * em35
       ninep637em5 = 9637. * em08
       eightp314em5 = 8314. * em08        
C
       onep1ep30 = onep1 * ep30
       onep2ep30 = onep2 * ep30 
C
       zep999ep31 = zep999 * ep31
C--------------------------------------
       zep54     = one   - zep46    
       zep1381966= zep13 + eightem3 + em04 + nineem5 + sixem6 
     +           + six  *em07  
       zep5854102= zep5 + zep08 + fiveem3 + fourem4 + em05
     +           + two * em07  
C--------------------------------------
        onep001em32  = 1001. * em35
        ninep637em5 = 9637. * em08
        onep3em5     = onep3  * em05
        eightp314em5 = 8314.* em08
        nineem2 = zep09        
C--------------------------------------
C      Fractions
C--------------------------------------
       zep426667 = hundred28 / three100
       zep222    = ten / fourty5
       zep00666666667 = two / three100
       zep01666666667 = five / three100
       zep444    = hundred / twohundred25
       onep222    = one + zep222
       onep8333   = one + twenty5 / thirty
       twop444  = two + zep444                  
       twop6666666667= eight / three 
       twop6667 = twop6666666667
       fivep333 = sixteen / three    
C--------------------------------------
C      Physics 
C--------------------------------------
       r_igc    = 8.314472
C--------------------------------------
C      Math 
C--------------------------------------
       sqr2    = sqrt(two)
       sqr3    = sqrt(three)
       sqr5    = sqrt(five)
       sqr6    = sqrt(six)
       sqr7    = sqrt(seven)
       sqr3_2  = sqrt(three_half)
       sqr4_3  = two * sqrt(third)
       sqr16_3 = two * sqr4_3
       pi      = atan2(zero,-one)
       invpi   = one / pi
c ------------------------------
        em3 = em03        
        em4 = em04 
        em5 = em05 
        em6 = em06 
        em7 = em07 
        em8 = em08 
        em9 = em09
        ep9 = ep09 
C ------------------------       
        fiveem2  = zep05
        twoem2  = zep02
        nineem2   = zep09
        zep80 = zep8              
        unp20 = onep2                                             
        unp30 = onep3
        unp70 = onep7
        unp90 = onep9
C-----------------------------   
       RETURN
    
       END