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 58: | Line 58: | ||
{{ :nuclide_chart.jpg?650 |}} | {{ :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 |}} | ||
+ |