OpenCores
URL https://opencores.org/ocsvn/openrisc/openrisc/trunk

Subversion Repositories openrisc

[/] [openrisc/] [trunk/] [gnu-dev/] [or1k-gcc/] [gcc/] [testsuite/] [gfortran.fortran-torture/] [compile/] [pr24136.f] - Blame information for rev 695

Details | Compare with Previous | View Log

Line No. Rev Author Line
1 695 jeremybenn
      subroutine  electra(ro,t,ye,ee,pe,se
2
     a                   ,eer,eet,per,pet,ser,set,keyps)
3
      implicit real*8 (a-h,o-z)
4
      common  /nunu/  nu,dnudr,dnudb,eta,detadnu,nup
5
      data facen,facpr,facs,rg /2.037300d+24,1.358200d+24,1.686304d-10
6
     1,8.314339d+07/
7
      data a1,a2,a3,a4 /2.059815d-03,-7.027778d-03
8
     1,4.219747d-02,-1.132427d+00/
9
      beta=facs*t
10
      b32=b12*beta
11
      u=(f62/f52)**2
12
      dudnu=2.0d0*u*(df62/f62-df52/f52)
13
      x=beta*u
14
      f=1.0d0+x*(2.5d0+x*(2.0d0+0.5d0*x))
15
      df=2.5d0+x*(4.0d0+1.5d0*x)
16
      dfdb=u*df
17
      fi32=f32+(f-1.0d0)*f52/u
18
      dfidnu=dfidu*dudnu+df32+(f-1.0d0)*df52/u
19
      dfidb=dfdb*f52/u
20
      dfidbet=dfidb+dfidnu*dnudb
21
      gs=sqrt(g)
22
      dg=0.75d0*gs
23
      dgdb=u*dg
24
      dgdu=beta*dg
25
      gi32=f32+(g-1.0d0)*f52/u
26
      dgidu=(u*dgdu-g+1.0d0)*f52/us
27
      dgidnu=dgidu*dudnu+df32+(g-1.0d0)*df52/u
28
      dgidb=dgdb*f52/u
29
      dgidbet=dgidb+dgidnu*dnudb
30
      dgidroe=dgidnu*dnudr
31
      em=facen*b52*fi32
32
      demdbet=facen*b32*(2.5d0*fi32+beta*dfidbet)
33
      dpmdbet=facpr*b32*(2.5d0*gi32+beta*dgidbet)
34
      demdroe=facen*b52*dfidroe
35
      dpmdroe=facpr*b52*dgidroe
36
      call  divine(nup,fp12,dfp12,s12)
37
      s42=2.0d0
38
      call  divine(nup,fp42,dfp42,s42)
39
      eer=(ye*(demdroe+depdroe)-(em+ep)/ro)/ro
40
      eet=facs*(demdbet+depdbet)/ro
41
      per=ye*(dpmdroe+dppdroe)
42
      pet=facs*(dpmdbet+dppdbet)
43
      end

powered by: WebSVN 2.1.0

© copyright 1999-2024 OpenCores.org, equivalent to Oliscience, all rights reserved. OpenCores®, registered trademark.