materials_for_mc 0.2.0

Material builder for Monte Carlo simulations
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84

import pytest
import gc
import sys
from materials_for_mc import Material, Materials, Config

def test_building_up_materials():
    mat1 = Material()
    mat1.add_nuclide("H", 1.0)
    
    mat2 = Material()
    mat2.add_nuclide("Fe", 0.5)
    
    # mats = Materials([mat1, mat2])
    mats= Materials()
    mats.append(mat1)
    mats.append(mat2)
    assert len(mats) == 2
    assert mats[0].nuclides == [("H", 1.0)]
    assert mats[1].nuclides == [("Fe", 0.5)]


def test_building_up_materials():
    mat1 = Material()
    mat1.add_nuclide("H", 1.0)
    
    mat2 = Material()
    mat2.add_nuclide("Fe", 0.5)
    
    # mats = Materials([mat1, mat2])
    mats= Materials([mat1, mat2])
    assert len(mats) == 2
    assert mats[0].nuclides == [("H", 1.0)]
    assert mats[1].nuclides == [("Fe", 0.5)]

def test_materials_data_reading():
    from materials_for_mc import Config
    
    # Set the cross-section paths in the global Config
    Config.set_cross_sections({"Li6": "tests/Li6.json", "Li7": "tests/Li7.json"})
    
    material1 = Material()
    material1.add_nuclide("Li6", 1.0)

    material2 = Material()
    material2.add_nuclide("Li7", 1.0)
    
    # Create a Materials collection
    materials = Materials([material1, material2])
    
    # Accessing data will automatically load nuclides from the global Config
    assert len(materials) == 2


def test_materials_be9_selective_temperature_union_no_extra():
    from materials_for_mc import Config, clear_nuclide_cache
    # Ensure clean cache so previous tests that may have loaded multiple temps don't leak
    clear_nuclide_cache()
    # Two materials at 294 K; first uses Be9 (has 294 & 300), second uses Fe56
    m1 = Material()
    m1.temperature = "294"
    m1.add_nuclide("Be9", 1.0)
    m2 = Material()
    m2.temperature = "294"
    m2.add_nuclide("Fe56", 1.0)
    Config.set_cross_sections({"Be9": "tests/Be9.json", "Fe56": "tests/Fe56.json"})
    mats = Materials([m1, m2])  # eager union should only request 294 for Be9
    # Explicitly trigger eager union load via read_nuclides_from_json
    mats.read_nuclides_from_json({"Be9": "tests/Be9.json", "Fe56": "tests/Fe56.json"})

def test_materials_be9_selective_temperature_union():
    from materials_for_mc import Config, clear_nuclide_cache
    # Ensure clean cache so previous tests that may have loaded multiple temps don't leak
    clear_nuclide_cache()
    # Two materials, both at 294 K, each referencing Be9
    m1 = Material()
    m1.temperature = "294"
    m1.add_nuclide("Be9", 1.0)
    m2 = Material()
    m2.temperature = "300"
    m2.add_nuclide("Be9", 2.0)
    Config.set_cross_sections({"Be9": "tests/Be9.json"})
    mats = Materials([m1, m2])
    # Explicitly trigger eager union load across temperatures 294 & 300
    mats.read_nuclides_from_json({"Be9": "tests/Be9.json"})