WPEC 2018 / SG43 Paris, France May

Transcription

1 WPEC 2018 / SG43 Paris, France May LLNL-PRES This work was performed under the auspices of the U.S. Department of Energy by under contract DE-AC52-07NA Lawrence Livermore National Security, LLC

2 For brevity I have removed most const and & from C++ examples For example const_iterator find( std::string const &a_label ) const ; is displayed as iterator find( std::string a_label ); 2

3 This talk will mainly focus on GIDI FUDGE already has a lot of stuff 3

4 GNDS has lots of nodes that store a list of objects of a particular type Examples: <reactions> node stores a list of <reaction> nodes <crosssection> node stores a list of 1d functions <distribution> node stores a list of P(µ,E E) s FUDGE and GIDI use a Suite class to store these: next slide The <styles> node has a special instance of the Suite class GNDS has links that need to be constructured/followed Example from a <resonancereactions> node <link href= /reactionsuite/reactions/reaction[@label='u239 + photon ]"/> FUDGE and GIDI use an Ancestry class inherited by all other classes: more later 4

5 The Suite class acts like the python OrderedDict from collections import OrderedDict That is, items can be accessed via their label (which must be unique) or their order in the Suited (i.e., index) For example the distribution node <angulartwobody label="eval >... </angulartwobody> <angulartwobody label="montecarlo >... </angulartwobody> <multigroup3d label="multigroup_0 >... </multigroup3d> <multigroup3d label="multigroup_1 >... </multigroup3d> distribution[ MultiGroup_0 ] or distribution[2] are the same Iteration returns items in the order that they were entered An important other method is add 5

6 Has three main members moniker the node name of the node (e.g., reaction in <reaction>) ancestor a pointer to the objects parent attribute used to construct an objects xlink. For example, each reaction needs its label added to its xlink to uniquely refer to it <reaction label= n + (U238_e14 -> U238) > href="/reactionsuite/reactions/reaction[@label='n + (U238_e14 -> U238)']" Has two main members 6

7 reactionsuite styles reactions sums reaction reaction reaction crosssection dist crosssection dist crosssection dist 7

8 std::string moniker( ); Ancestry *ancestor( ); std::string attribute( ); void setancestor( Ancestry *a_ancestor ); Ancestry *root( ); bool ischild( Ancestry *a_child ); bool isparent( Ancestry *a_parent ); bool isroot( ); Ancestry *findinancestry( std::string a_href ); Ancestry *findinancestry2( std::size_t a_index, std::vector<std::string> a_segments ); virtual Ancestry *findinancestry3( std::string a_item ); virtual std::string xlinkitemkey( ); std::string toxlink( ); // Relative and absolute 8

9 void Suite::add( Form *a_form ) { int i1 = 0; for( Suite::iterator iter = begin( ); iter!= end( ); ++iter, ++i1 ) { if( (*iter)->label( ) == a_form->label( ) ) { m_forms[i1] = a_form; return; } } m_map[a_form->label( )] = (int) m_forms.size( ); m_forms.push_back( a_form ); } a_form->setancestor( this ); // Oops, missing a_form->setancestor( this ); in replacement. 9

10 The styles node is a Suite instance but adds the additional method findlabelinlineage std::string *findlabelinlineage( GIDI::Suite a_suite, std::string a_label ); Used to find nearest derived data. For example, consider elastic neutron and proton distributions with multi-group upscattering <angulartwobody label="eval >... </angulartwobody> <multigroup3d label="multigroup_0 >... </multigroup3d> <multigroup3d label="upscat_0 >... </multigroup3d> <angulartwobody label="eval >... </angulartwobody> <multigroup3d label="multigroup_0 >... </multigroup3d> findlabelinlineage( distribution, Upscat_0 ) 10

11 In FUDGE when a child node is added its parent, its parent via setancestror takes ownership There can be an issue if another object also contains the child Example, let crosssection1 and crosssection2 be to crosssection nodes and let data1 be a cross section representation. Now consider the following crosssection1.add( data1 ) crosssection2.add( data1 ) # data1.ancestor = crosssection1 # data1.ancestor = crosssection2 Changing data1 changes both. Changing via crosssection1 also changes both Correct way to do this in FUDGE crosssection1.add( data1 ) crosssection2.add( data1.copy( ) ) 11

12 See Caleb s talk GIDI only does reading Reading XML is slow even with pugixml MCGIDI is passed a GIDI instance FUDGE does reading and writing 12

13 One vs many XYs1d classes GIDI only has one XYs1d class FUDGE has one for every physical property (e.g., cross section, multiplicity, P(u), P(E )) For example crosssection.py as a XYs1d multiplitiy.py as a XYs1d They all inherit from a base XYs1d The cross section suite only allows the XYs1d from crosssection.py to be added to it. Each XYs1d may add/over write additional methods For example, the cross section XYs1d adds the heat method which does not make sense for multiplity 13

14 All reside in the namespace Functions Base - included by all functions axes, dimension, interpolation, etc. members domainmin, domainmax virtual methods Base1d included by all 1d functions As evaluate( double x1 ) method Base2d included by all 2d functions As evaluate( double x2, double x1 ) method Base3d included by all 3d functions As evaluate( double x3, double x2, double x1 ) method 14

15 All reside in the namespace Probabilities Base - included by all probabilities Inherits from Functions::Base Add Xs member Base1d included by all 1d functions As methods evaluate( double x1 ) and sample( double rngvalue, Sampling::Input &a_input ) Note, rngvalue is x0 (y) value and x1 is sampled. Base2d included by all 2d functions As methods evaluate and sample( double x2, double rngvalue, Sampling::Input &a_input ) Base3d included by all 3d functions As evaluate and sample methods 15

16 reaction.crosssection( ) returns a reaction s cross section suite. For example for crosssection in reaction.crosssection reaction.crosssection( Settings::MG settings ) returns a multi-group cross section for reactions 16

17 Styles::Suite styles( ); Suite reactions( ) Suite orphanproducts( ); void productindices( std::set<int> Indices, bool transportablesonly ); int maximumlegendreorder( Settings::MG settings, int productindex ); Styles::TemperatureInfos temperatures( ); bool hasfission( ); 17

18 MG derived class Used to get multi-group data Vector crosssection( Settings::MG a_settings ) Defines style to get data from, multi-group collapsing info, want delayed neutrons included, etc. 18

19 std::vector<double> groupboundaries( Settings::MG settings, int productindex ); Vector inversespeed( Settings::MG settings ); Vector crosssection( Settings::MG settings ); Vector multiplicity( Settings::MG settings, int productindex ); Vector Q( Settings::MG settings, bool final ); Matrix productmatrix( Settings::MG settings, int productindex, int order ); Matrix fissionmatrix( Settings::MG settings, int order ); Vector transportcorrection( Settings::MG settings, int order, transportcorrectiontype transportcorrectiontype ); Vector availableenergy( Settings::MG settings ); Vector averageenergy( Settings::MG settings, int productindex ); Vector depositionenergy( Settings::MG settings ); Vector availablemomentum( Settings::MG settings ); Vector averagemomentum( Settings::MG settings, int productindex ); Vector depositionmomentum( Settings::MG settings ); 19

20 Vector crosssection( Settings::MG settings ); Vector multiplicity( Settings::MG settings, int productindex ); Vector Q( Settings::MG settings, bool final ); Matrix productmatrix( Settings::MG settings, int productindex, int order ); Matrix fissionmatrix( Settings::MG settings, int order ); Vector availableenergy( Settings::MG settings ); Vector averageenergy( Settings::MG settings, int productindex ); Vector depositionenergy( Settings::MG settings ); Vector availablemomentum( Settings::MG settings ); Vector averagemomentum( Settings::MG settings, int productindex ); Vector depositionmomentum( Settings::MG settings ); All of these methods are called by their Protare equivalent methods. The Product class has many of these methods. 20