feat: add read omega block

This commit is contained in:
Stig Rune Jensen
2025-02-12 16:23:44 +01:00
parent 02de20d50b
commit 9abd9d45bc

View File

@@ -165,6 +165,17 @@ let readOmega (ds: DataSet) t l =
Log.Error $"{e}"
Array.empty
let readOmegaBlock (ds: DataSet) t =
try
let n = ds.Dimensions["node"].Length
let l = ds.Dimensions["siglev"].Length
let u = ds[ "omega" ].GetData([| t; 0; 0 |], [| 1; l; n |]) :?> single [,,]
u[0, *, *]
with
| e ->
Log.Error $"{e}"
Array2D.zeroCreate 0 0
let readOmegas (ds: DataSet) t l es =
try
let n = ds.Dimensions["node"].Length