This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision Next revision Both sides next revision | ||
nu_rnet [2014/06/05 15:45] schutrumpg |
nu_rnet [2014/06/05 17:17] schutrumpg |
||
---|---|---|---|
Line 56: | Line 56: | ||
Last but not least we need a reaction library. This was taken from the JINA Reaclib database. For the r-process only nucleides have to be considered which are stable or even more neutron rich. The proton rich nucleides can be neglected. In the chart of nuclei the considered rates look like this: | Last but not least we need a reaction library. This was taken from the JINA Reaclib database. For the r-process only nucleides have to be considered which are stable or even more neutron rich. The proton rich nucleides can be neglected. In the chart of nuclei the considered rates look like this: | ||
- | {{ :nuclide_chart.jpg?700 |}} | + | {{ :nuclide_chart.jpg?650 |}} |
The black squares mark the stable isotopes. | The black squares mark the stable isotopes. | ||
+ | |||
+ | This network includes 4510 species and 53763 reaction rates. This setup requires a very fast matrix solver. | ||
+ | |||
+ | ===Matrix Solver=== | ||
+ | |||
+ | The standart matrix solver included in the xnet programm "LAPACK" is too slow to solve a network like this. Faster solutions are sparse matrix solver. Two interfaces are included in xnet: For MA48 and Pardiso. Instructions to include these matrix solvers in the code can be found in the documentation of the code. We experienced problems for both with Linux: Pardiso should be much faster than MA48 needs a full compiled LAPACK library. This is included in the ifort compiler. Finally we succeeded in compiling Pardiso version with ifort. The program can be found here. | ||
+ | |||
+ | ====Results==== | ||
+ | ==$Y_e=0.2$== | ||
+ | {{ :ev20.png?400 |}} | ||
+ | ==$Y_e=0.25$== | ||
+ | {{ :ev25.png?400 |}} | ||
+ | ==$Y_e=0.4$== | ||
+ | {{ :ev40.png?400 |}} | ||
+ |