1 files changed, 84 insertions, 0 deletions
diff --git a/src/grt/grt-analog_solver.adb b/src/grt/grt-analog_solver.adb
new file mode 100644
index 000000000..c4d2c3b39
--- /dev/null
+++ b/src/grt/grt-analog_solver.adb
@@ -0,0 +1,84 @@
+--  GHDL Run Time (GRT) -  analog solver for sundials.
+--  Copyright (C) 2022 Tristan Gingold
+--
+--  This program is free software: you can redistribute it and/or modify
+--  it under the terms of the GNU General Public License as published by
+--  the Free Software Foundation, either version 2 of the License, or
+--  (at your option) any later version.
+--
+--  This program is distributed in the hope that it will be useful,
+--  but WITHOUT ANY WARRANTY; without even the implied warranty of
+--  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+--  GNU General Public License for more details.
+--
+--  You should have received a copy of the GNU General Public License
+--  along with this program.  If not, see <gnu.org/licenses>.
+--
+--  As a special exception, if other files instantiate generics from this
+--  unit, or you link this unit with other files to produce an executable,
+--  this unit does not by itself cause the resulting executable to be
+--  covered by the GNU General Public License. This exception does not
+--  however invalidate any other reasons why the executable file might be
+--  covered by the GNU Public License.
+
+with Grt.Errors; use Grt.Errors;
+
+package body Grt.Analog_Solver is
+   function Sundials_Init (Sz : Ghdl_I32) return Ghdl_I32;
+   pragma Import (C, Sundials_Init, "grt_sundials_init");
+
+   function Sundials_Start return Ghdl_I32;
+   pragma Import (C, Sundials_Start, "grt_sundials_start");
+
+   procedure Sundials_Solve (T0 : Ghdl_F64; Tn : in out Ghdl_F64;
+                                            Res : out Integer);
+   pragma Import (C, Sundials_Solve, "grt_sundials_solve");
+
+   procedure Init (Size : Ghdl_I32)
+   is
+      Res : Ghdl_I32;
+   begin
+      Res := Sundials_Init (Size);
+      if Res < 0 then
+         Internal_Error ("sundials initialization failure");
+      end if;
+   end Init;
+
+   function Sundials_Get_Yy_Vec return F64_C_Arr_Ptr;
+   pragma Import (C, Sundials_Get_Yy_Vec, "grt_sundials_get_yy_vec");
+
+   function Get_Init_Val_Ptr return F64_C_Arr_Ptr is
+   begin
+      return Sundials_Get_Yy_Vec;
+   end Get_Init_Val_Ptr;
+
+   function Sundials_Get_Yp_Vec return F64_C_Arr_Ptr;
+   pragma Import (C, Sundials_Get_Yp_Vec, "grt_sundials_get_yp_vec");
+
+   function Get_Init_Der_Ptr return F64_C_Arr_Ptr is
+   begin
+      return Sundials_Get_Yp_Vec;
+   end Get_Init_Der_Ptr;
+
+   procedure Start is
+   begin
+      if Sundials_Start /= 0 then
+         Internal_Error ("sundials start");
+      end if;
+   end Start;
+
+   procedure Set_Root_Size (Size : Ghdl_I32) is
+   begin
+      Internal_Error ("sundials set_root_size");
+   end Set_Root_Size;
+
+   procedure Solve (T : Ghdl_F64; Tn : in out Ghdl_F64; Res : out Integer) is
+   begin
+      Sundials_Solve (T, Tn, Res);
+   end Solve;
+
+   procedure Finish is
+   begin
+      Internal_Error ("sundials finish");
+   end Finish;
+end Grt.Analog_Solver;