Package com.nr.sf

Examples of com.nr.sf.Fermi.val()

63
64
65
66
67
68
69
70
71
72
73
    // Test cases with theta=0
    // Fermi(k,eta,0)=exp(eta)*LerchPhi(exp(eta),k+1,1)
    Fermi dirac = new Fermi();

    for (i=0;i<N;i++) {
      zz[i]=dirac.val(k[i],eta[i],0.0);
      expect[i]=exp(gammln(k[i]+1.0))*exp(eta[i])*lerchphi[i];
//      System.out.printf(setw(15) << expect[i] << setw(15) << zz[i];
//      System.out.printf(setw(15) << (zz[i]/expect[i]-1.0));
    }

View Full Code Here
83
84
85
86
87
88
89
90
91
92
93
    // Limiting cases for large eta and k=1/2
    maxerr=0.0;
    for (i=0;i<N;i++) {
      arg=1000.0;
      theta=0.1*(i+2);
      zz[i]=dirac.val(0.5,arg,theta);
      u=1.0+arg*theta;
      y=sqrt(SQR(u)-1.0);
      x=log((y+sqrt(SQR(y)+4.0))/2.0);
      expect[i]=(y*u-x)/pow(2.0*theta,1.5);
      err=abs(zz[i]/expect[i]-1.0);
View Full Code Here
TOP
Copyright © 2018 www.massapi.com. All rights reserved.
All source code are property of their respective owners. Java is a trademark of Sun Microsystems, Inc and owned by ORACLE Inc. Contact coftware#gmail.com.