las_rs 0.2.1

High-performance parser and writer for LAS (Log ASCII Standard) well-log files — LAS 1.2 / 2.0 / 3.0. Pure-Rust core; the optional `python` feature adds the PyO3 bindings that back the `las-rs` PyPI wheel.
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
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118

"""Phase 5 extended: Pickle serialization for all core types."""

import pickle

import numpy as np
import pytest

import las_rs


# ---------------------------------------------------------------------------
# Helpers
# ---------------------------------------------------------------------------

_SIMPLE_LAS = """\
~VERSION INFORMATION
 VERS.   2.0 : LAS VERSION 2.0
 WRAP.    NO : ONE LINE PER DEPTH STEP
~WELL INFORMATION
 STRT.M  3000.0 : START DEPTH
 STOP.M  3004.0 : STOP DEPTH
 STEP.M     1.0 : STEP VALUE
 NULL. -999.25  : NULL VALUE
 COMP.  MESA OIL CORP  : COMPANY
 WELL.  FLATROCK-11 #2 : WELL NAME
~CURVE INFORMATION
 DEPT.M    : MEASURED DEPTH
 GR  .GAPI : GAMMA RAY
 RHOB.G/CC : BULK DENSITY
~ASCII LOG DATA
 3000.0   48.31  2.508
 3001.0   62.74  2.490
 3002.0   77.08  2.533
 3003.0   55.20  2.461
 3004.0   41.65  2.576
"""


# ---------------------------------------------------------------------------
# 1. test_pickle_lasfile
# ---------------------------------------------------------------------------

@pytest.mark.xfail(reason="not yet implemented")
def test_pickle_lasfile():
    """pickle.dumps / pickle.loads round-trips a LASFile, preserving the
    version string, well section, curve names, and numeric data."""
    original = las_rs.read(_SIMPLE_LAS)
    blob = pickle.dumps(original)
    restored = pickle.loads(blob)

    # Version preserved
    assert float(restored.version["VERS"].value) == pytest.approx(
        float(original.version["VERS"].value)
    )
    # Well name preserved
    assert restored.well["WELL"].value == original.well["WELL"].value
    # Curve count preserved
    assert len(restored.curves) == len(original.curves)
    # Numeric data preserved
    np.testing.assert_array_almost_equal(
        restored.curves["GR"].data,
        original.curves["GR"].data,
    )
    np.testing.assert_array_almost_equal(
        restored.curves["DEPT"].data,
        original.curves["DEPT"].data,
    )


# ---------------------------------------------------------------------------
# 2. test_pickle_curveitem
# ---------------------------------------------------------------------------

@pytest.mark.xfail(reason="not yet implemented")
def test_pickle_curveitem():
    """pickle.dumps / pickle.loads round-trips a CurveItem, preserving all
    header fields and the data array."""
    samples = np.array([10.0, 20.5, 31.2, 42.8, 55.1])
    original = las_rs.CurveItem(
        mnemonic="SP",
        unit="MV",
        value="52 530 01 00 10",
        descr="Spontaneous potential",
        data=samples,
    )
    blob = pickle.dumps(original)
    restored = pickle.loads(blob)

    assert restored.mnemonic == "SP"
    assert restored.unit == "MV"
    assert restored.value == "52 530 01 00 10"
    assert restored.descr == "Spontaneous potential"
    np.testing.assert_array_almost_equal(restored.data, samples)


# ---------------------------------------------------------------------------
# 3. test_pickle_sectionitems
# ---------------------------------------------------------------------------

@pytest.mark.xfail(reason="not yet implemented")
def test_pickle_sectionitems():
    """pickle.dumps / pickle.loads round-trips a SectionItems, preserving all
    items and the mnemonic_transforms flag."""
    original = las_rs.SectionItems(mnemonic_transforms=True)
    original.append(las_rs.HeaderItem(mnemonic="WELL",  value="PICKLE-WELL"))
    original.append(las_rs.HeaderItem(mnemonic="FLD",   value="PICKLE-FIELD"))
    original.append(las_rs.HeaderItem(mnemonic="PROV",  value="PICKLE-PROVINCE"))

    blob = pickle.dumps(original)
    restored = pickle.loads(blob)

    assert len(restored) == 3
    assert list(restored.keys()) == ["WELL", "FLD", "PROV"]
    assert restored["WELL"].value == "PICKLE-WELL"
    assert restored["FLD"].value  == "PICKLE-FIELD"
    assert restored["PROV"].value == "PICKLE-PROVINCE"
    # mnemonic_transforms flag must be preserved
    assert restored.mnemonic_transforms is True