A package for ANT-MOC data manipulation.
Project description
ANT-MOC Data
A package for ANT-MOC data manipulation.
Prerequisites
- Python >= 3.6
- baseopt
- numpy
- h5py
Install
$ pip install antmocdata
License
MIT
ANT-MOC Log
Package antmocdata.log
contains tools for exploring ANT-MOC logs.
Examples
Please check the directory examples/log
for live examples.
Each of the sample scripts accepts command line arguments.
python ./examples/log/extract-records.py --help
Log file
A log file of ANT-MOC contains many data fields. If a group of log files are provided to a LogDB
object, these files will be serialized and stored as LogFile
objects.
from antmocdata.log import Options, LogDB
options = Options()
# ...
# setup options eigher through CLI or direct assignments
# ...
logdb = LogDB()
logdb.setup(options)
Extractor
A log extractor is used to make LogDB
queries and save the results to a .csv
file. For example, to get log files with specific fields, one could write down
options["output"].value = "antmoc-records.csv"
options["specs"].value = ["Azims", "Polars"]
extractor = TinyExtractor(logdb)
extractor.extract()
This would list the Azims
and Polars
fields of all the log files and save the results to antmoc-records.csv
.
In this case, Azims
and Polars
are called FieldSpec
. A spec can be used to filter out results. For example,
# List Azims and Polars fields of all the logs
options["specs"].value = ["Azims", "Polars"]
# List Azims and Polars fields of all the logs, and
# only show logs with Azims=64
options["specs"].value = ["Azims=64", "Polars"]
# List Azims and Polars fields of all the logs, and
# only show logs with Azims=64 and Polars>2
options["specs"].value = ["Azims=64", "Polars>2"]
extractor.extract()
A FieldSpec
consists of three parts: field name, binary operator, and value. Perl regex is supported for the field name and the operator. For example, ".*Time"
could match all the fields with an ending Time
string.
The binary operator could be ==
, <
, <=
, >
, or >=
. Operator ==
is for string comparison and perl regex is supported in this case. Inequality symbols are for string or numerical comparisons.
Be careful if you want to use inequality symbols on string fields. Values from these fields are compared through string comparison. Predefined field types are located in
antmocdata.log.default_fields.json
.
Fields
Predefined log fields are located in antmocdata.log.default_fields.json
. If these fields are outdated due to ANT-MOC updates, please update this json file or load a new one in your scripts.
from antmocdata.log import fields
fields.load("path/to/your/fields.json")
There is also an add
method for appending single field to the field dictionary.
from antmocdata.log import Field
fields.add(Field(name="NewField1", patterns=["NewField1.*"]))
# or adding the field directly to LogFields
from antmocdata.log import LogFields
LogFields().add(Field(name="NewField2", patterns=["NewField2.*"]))
Save log database
A LogDB
object can be dumped as json files to a specific directory.
options["savedb"].value = "antmoc-logdb/"
# ...
# setup the LogDB object
# ...
logdb.save(options("savedb"))
ANT-MOC MGXS
Package antmocdata.mgxs
contains tools for checking, manipulating, and generating MGXS files for ANT-MOC.
Examples
Please check the directory examples/mgxs
for live examples.
Each of the sample scripts accepts command line arguments.
python ./examples/mgxs/h5/fix-materials.py --help
python ./examples/mgxs/xml/fix-materials-in-xml.py --help
HDF5 data layout
There are two layouts of material data in an H5 file. Materials are treated as data groups in both of the layouts.
In addition to materials, the H5 file must contain a top-level attribute named '# groups' for the number of energy groups.
Layout: named
This is the default cross-section data layout for ANT-MOC.
Cross-section arrays are stored in H5 datasets individually. In this layout, a group for a material consists of several datasets:
absorption
fission
nu-fission
transport
, ortotal
chi
scatter matrix
, ornu-scatter matrix
, orconsistent scatter matrix
For example, a simple cross-section file could have the following hierarchy
attribute "# groups"
group /
group /material/MOX-8.7%
dataset /material/MOX-8.7%/chi
dataset /material/MOX-8.7%/fission
dataset /material/MOX-8.7%/nu-fission
dataset /material/MOX-8.7%/scatter matrix
dataset /material/MOX-8.7%/total
group /material/UO2
dataset /material/UO2/chi
dataset /material/UO2/fission
dataset /material/UO2/nu-fission
dataset /material/UO2/scatter matrix
dataset /material/UO2/total
The scatter matrix is usually a flattened n-by-n matrix, where n is the number of energy groups. Elements in the scatter dataset respect the source-major order, which is much like the row-major order.
For example, a scatter matrix with 2 energy groups has 4 elements, which are stored as
1->1
1->2
2->1
2->2
The number before symbol
->
is the source group, and the number after the symbol is the destination group.
Layout: compact/compressed
TODO
Common modules
material
: classMaterial
representing cross-sections. A material object could be written to an HDF5 file as a dataset.materialxml
: representation for the XML material definition, which is used to handlematerials.xml
.manip
: data manipulation utilities.options
: representation of command line options.
Type A
Package antmocdata.mgxs.type_a
defines a generator which accepts two files to create an mgxs input for antmoc:
infilecross
: cross-sections in plain text.materials.xml
: material definitions in XML, including nuclear densities.
Modules
material
: definition ofMaterialTypeA
, which is a sub-class ofMaterial
.nuclides
: representations of nuclides and nuclide sets, which are basically defined in a plain text file called "infilecross".infilecross
: functions for parsing an "infilecross" file. The file must be well-formed.generate
: functions for mgxs generation.options
: definition ofOptionsTypeA
, which is a sub-class ofOptions
.
File formats
materials.xml
This is an XML file consisting of material definitions.
<?xml version="1.0" encoding="utf-8"?>
<MATERIAL>
<material name="1" set="1" density="0.0E+00" temperature="523.15K" label="Some material">
<nuclide id="1102301" radio="1.6098e-2"/>
<nuclide id="601201" radio="7.3771e-5"/>
</material>
</MATERIAL>
material
: definition of a material object.material.name
: material name which will be written into the H5 file (string).material.set
: nuclide set ID for MGXS calculations (int, defaults toname
).material.label
: a short description.nuclide
: nuclide information for MGXS calculations.nuclide.id
: nuclide ID containing its atomic number and mass (int).nuclide.radio
: density used by MGXS calculations (float).
infilecross
This is a plain text file consisting of nuclide set definitions.
$SOMESTRING 1 7SETs 1SET
2 30 6 5 0 0
0.57564402E+00 0.29341474E+00 0.12731624E+00 ...
...
2511102302 0 3 6SODIUM-23
...
4311402802 0 3 6SILICON-28
...
$SOMESTRING 1 7SETs SET2
...
Type B
TODO
Modules
File formats
Project details
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
Built Distribution
Hashes for antmocdata-0.1.2-py3-none-any.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | 1e5f244a61e796bd0f165beb7608dc6692293b8d53b0340a780c234738beba2c |
|
MD5 | 3298870a27db65535b8da148fbf63119 |
|
BLAKE2b-256 | 105be27ed96828c036a26bdce198d693414ca2a32db3d2c1092cfe5df7026c75 |