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 16:34] schutrumpg |
nu_rnet [2014/06/05 17:21] schutrumpg |
||
---|---|---|---|
Line 14: | Line 14: | ||
{{:core-collapse_figure.png?650|}} | {{:core-collapse_figure.png?650|}} | ||
- | + | Fig: Schematic representation of the evolutionary stages from stellar core collapse. (Ref. H.Th. Janka, arXiv:astro-ph/0612072) | |
- | Fig. Schematic representation of the evolutionary stages from stellar core collapse. (Ref. H.Th. Janka, arXiv:astro-ph/0612072) | + | |
=====Method and aims:===== | =====Method and aims:===== | ||
Line 60: | Line 58: | ||
This network includes 4510 species and 53763 reaction rates. This setup requires a very fast matrix solver. | 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 |}} | ||
+ |