sparse-ir-capi 0.8.4

C API for SparseIR Rust implementation
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

# ---
# jupyter:
#   jupytext:
#     cell_metadata_filter: -all
#     custom_cell_magics: kql
#     text_representation:
#       extension: .jl
#       format_name: percent
#       format_version: '1.3'
#       jupytext_version: 1.11.2
#   kernelspec:
#     display_name: Julia 1.12.1
#     language: julia
#     name: julia-1.12
# ---

# %%
using Libdl: dlext
# Load the shared library"
const libpath = joinpath(@__DIR__, "../../target/debug/libsparse_ir_capi.$(dlext)")

# %%
mutable struct spir_kernel end

"""
Create a Logistic kernel
"""
function kernel_logistic_new(lambda::Float64)
    status = Ref{Int32}(0)
    kernel = ccall(
        (:spir_logistic_kernel_new, libpath),
        Ptr{spir_kernel},
        (Float64, Ref{Int32}),
        lambda, status
    )

    if kernel == C_NULL
        error("Failed to create kernel: status = $(status[])")
    end

    return kernel
end

kernel = kernel_logistic_new(10.0)

# %%
epsilon_coarse = 1e-6
epsilon_fine = 1e-10

ngauss_coarse = Ref{Cint}(0)
ngauss_fine = Ref{Cint}(0)

ccall(
    (:spir_kernel_get_sve_hints_ngauss, libpath),
    Cint,
    (Ptr{spir_kernel}, Cdouble, Ptr{Cint}),
    kernel, epsilon_coarse, ngauss_coarse
)
@assert ngauss_coarse[] == 10
ccall(
    (:spir_kernel_get_sve_hints_ngauss, libpath),
    Cint,
    (Ptr{spir_kernel}, Cdouble, Ptr{Cint}),
    kernel, epsilon_fine, ngauss_fine
)
@assert ngauss_fine[] == 16

# %%
epsilon = 1e-8

nsvals = Ref{Cint}(0)
ccall(
    (:spir_kernel_get_sve_hints_nsvals, libpath),
    Cint,
    (Ptr{spir_kernel}, Cdouble, Ptr{Cint}),
    kernel, epsilon, nsvals
)

@assert nsvals[] > 0