Skip to content

Commit

Permalink
Fix absolute to relative in notebook_sampling_ify_u235th_cea_c1.ipynb (
Browse files Browse the repository at this point in the history
…#309)

I have divided the absolute covariance matrix elementwise by the correct independent fission yields to convert it to a relative covariance matrix to draw samples from.
  • Loading branch information
Lars-Engelen authored May 22, 2024
1 parent c2fb17e commit c1fc6cb
Showing 1 changed file with 3 additions and 1 deletion.
4 changes: 3 additions & 1 deletion notebooks/notebook_sampling_ify_u235th_cea_c1.ipynb
Original file line number Diff line number Diff line change
Expand Up @@ -46,9 +46,11 @@
"\n",
"idx = nfpy.data.query(f\"E==0.0253 & MT==454\").index\n",
"ifyth = nfpy.data.loc[idx]\n",
"IFY_np = np.array([ifyth['FY'].to_numpy()])\n",
"\n",
"cov = pd.read_csv(join(path, file_cov), sep=r\"\\s+\", header=None)\n",
"cov.index = cov.columns = pd.Index(ifyth.ZAP)\n",
"cov = np.divide(cov.data.to_numpy(), IFY_np.T @ IFY_np) # convert to relative terms\n",
"cov = pd.DataFrame(cov, index = pd.Index(ifyth.ZAP), columns = pd.Index(ifyth.ZAP))\n",
"cov = sandy.CategoryCov(cov)"
]
},
Expand Down

0 comments on commit c1fc6cb

Please sign in to comment.