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] - Rev 742
Go to most recent revision | Compare with Previous | Blame | View Log
subroutine electra(ro,t,ye,ee,pe,se a ,eer,eet,per,pet,ser,set,keyps) implicit real*8 (a-h,o-z) common /nunu/ nu,dnudr,dnudb,eta,detadnu,nup data facen,facpr,facs,rg /2.037300d+24,1.358200d+24,1.686304d-10 1,8.314339d+07/ data a1,a2,a3,a4 /2.059815d-03,-7.027778d-03 1,4.219747d-02,-1.132427d+00/ beta=facs*t b32=b12*beta u=(f62/f52)**2 dudnu=2.0d0*u*(df62/f62-df52/f52) x=beta*u f=1.0d0+x*(2.5d0+x*(2.0d0+0.5d0*x)) df=2.5d0+x*(4.0d0+1.5d0*x) dfdb=u*df fi32=f32+(f-1.0d0)*f52/u dfidnu=dfidu*dudnu+df32+(f-1.0d0)*df52/u dfidb=dfdb*f52/u dfidbet=dfidb+dfidnu*dnudb gs=sqrt(g) dg=0.75d0*gs dgdb=u*dg dgdu=beta*dg gi32=f32+(g-1.0d0)*f52/u dgidu=(u*dgdu-g+1.0d0)*f52/us dgidnu=dgidu*dudnu+df32+(g-1.0d0)*df52/u dgidb=dgdb*f52/u dgidbet=dgidb+dgidnu*dnudb dgidroe=dgidnu*dnudr em=facen*b52*fi32 demdbet=facen*b32*(2.5d0*fi32+beta*dfidbet) dpmdbet=facpr*b32*(2.5d0*gi32+beta*dgidbet) demdroe=facen*b52*dfidroe dpmdroe=facpr*b52*dgidroe call divine(nup,fp12,dfp12,s12) s42=2.0d0 call divine(nup,fp42,dfp42,s42) eer=(ye*(demdroe+depdroe)-(em+ep)/ro)/ro eet=facs*(demdbet+depdbet)/ro per=ye*(dpmdroe+dppdroe) pet=facs*(dpmdbet+dppdbet) end
Go to most recent revision | Compare with Previous | Blame | View Log