sundials-rs 0.1.0

Safe Rust bindings to the SUNDIALS ODE/DAE solver library (CVODE, CVODES, IDA, IDAS)
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

use std::ptr::NonNull;
use sundials_rs_sys as sys;
use crate::error::SundialsError;

/// Wraps a dense `SUNMatrix`.
pub struct DenseMatrix {
    pub(crate) ptr: NonNull<sys::_generic_SUNMatrix>,
}

impl DenseMatrix {
    /// Allocate an `m × n` dense matrix.
    pub(crate) fn new(m: usize, n: usize, ctx: sys::SUNContext) -> Result<Self, SundialsError> {
        let ptr = unsafe {
            sys::SUNDenseMatrix(
                m as sys::sunindextype,
                n as sys::sunindextype,
                ctx,
            )
        };
        let ptr = NonNull::new(ptr).ok_or(SundialsError::MatrixCreate)?;
        Ok(Self { ptr })
    }
}

impl Drop for DenseMatrix {
    fn drop(&mut self) {
        unsafe { sys::SUNMatDestroy(self.ptr.as_ptr()) };
    }
}

unsafe impl Send for DenseMatrix {}