pydae 0.57.1__py3-none-any.whl → 0.57.3__py3-none-any.whl

pydae/bmapu/lines/temp_trap_cffi.c

@@ -573,52 +573,60 @@ static void (*_cffi_call_python_org)(struct _cffi_externpy_s *, char *);
 void de_jac_trap_up_eval(double *data,double *x,double *y,double *u,double *p,double Dt){
 
 data[0] = 0.5*Dt + 1; 
-data[13] = 0.5*Dt*p[8] + 1; 
-data[18] = 0.5*Dt; 
-data[133] = p[7]; 
-data[138] = -p[6]; 
+data[15] = 0.5*Dt*p[8] + 1; 
+data[26] = 0.5*Dt; 
+data[183] = p[7]; 
+data[194] = -p[6]; 
 
 }
 
 void de_jac_trap_xy_eval(double *data,double *x,double *y,double *u,double *p,double Dt){
 
-data[50] = y[2]*(p[2]*sin(y[1] - y[3]) - p[1]*cos(y[1] - y[3])); 
-data[51] = y[0]*y[2]*(p[2]*cos(y[1] - y[3]) + p[1]*sin(y[1] - y[3])); 
-data[52] = y[0]*(p[2]*sin(y[1] - y[3]) - p[1]*cos(y[1] - y[3])) + 2*y[2]*p[1]; 
-data[53] = y[0]*y[2]*(-p[2]*cos(y[1] - y[3]) - p[1]*sin(y[1] - y[3])); 
-data[62] = y[2]*(p[2]*cos(y[1] - y[3]) + p[1]*sin(y[1] - y[3])); 
-data[63] = y[0]*y[2]*(-p[2]*sin(y[1] - y[3]) + p[1]*cos(y[1] - y[3])); 
-data[64] = y[0]*(p[2]*cos(y[1] - y[3]) + p[1]*sin(y[1] - y[3])) + 2*y[2]*(-p[2] - 1.0/2.0*p[3]); 
-data[65] = y[0]*y[2]*(p[2]*sin(y[1] - y[3]) - p[1]*cos(y[1] - y[3])); 
-data[86] = 2*y[0]*p[1] - y[2]*(p[2]*sin(y[1] - y[3]) + p[1]*cos(y[1] - y[3])); 
-data[87] = -y[0]*y[2]*(p[2]*cos(y[1] - y[3]) - p[1]*sin(y[1] - y[3])); 
-data[88] = -y[0]*(p[2]*sin(y[1] - y[3]) + p[1]*cos(y[1] - y[3])); 
-data[89] = -y[0]*y[2]*(-p[2]*cos(y[1] - y[3]) + p[1]*sin(y[1] - y[3])); 
-data[98] = -2*y[0]*p[2] - y[0]*p[3] - y[2]*(-p[2]*cos(y[1] - y[3]) + p[1]*sin(y[1] - y[3])); 
-data[99] = -y[0]*y[2]*(p[2]*sin(y[1] - y[3]) + p[1]*cos(y[1] - y[3])); 
-data[100] = -y[0]*(-p[2]*cos(y[1] - y[3]) + p[1]*sin(y[1] - y[3])); 
+data[58] = y[2]*(p[2]*sin(y[1] - y[3]) - p[1]*cos(y[1] - y[3])); 
+data[59] = y[0]*y[2]*(p[2]*cos(y[1] - y[3]) + p[1]*sin(y[1] - y[3])); 
+data[60] = y[0]*(p[2]*sin(y[1] - y[3]) - p[1]*cos(y[1] - y[3])) + 2*y[2]*p[1]; 
+data[61] = y[0]*y[2]*(-p[2]*cos(y[1] - y[3]) - p[1]*sin(y[1] - y[3])); 
+data[72] = y[2]*(p[2]*cos(y[1] - y[3]) + p[1]*sin(y[1] - y[3])); 
+data[73] = y[0]*y[2]*(-p[2]*sin(y[1] - y[3]) + p[1]*cos(y[1] - y[3])); 
+data[74] = y[0]*(p[2]*cos(y[1] - y[3]) + p[1]*sin(y[1] - y[3])) + 2*y[2]*(-p[2] - 1.0/2.0*p[3]); 
+data[75] = y[0]*y[2]*(p[2]*sin(y[1] - y[3]) - p[1]*cos(y[1] - y[3])); 
+data[86] = -2*y[0]*p[1] - y[2]*(-p[2]*sin(y[1] - y[3]) - p[1]*cos(y[1] - y[3])); 
+data[87] = -y[0]*y[2]*(-p[2]*cos(y[1] - y[3]) + p[1]*sin(y[1] - y[3])); 
+data[88] = -y[0]*(-p[2]*sin(y[1] - y[3]) - p[1]*cos(y[1] - y[3])); 
+data[89] = -y[0]*y[2]*(p[2]*cos(y[1] - y[3]) - p[1]*sin(y[1] - y[3])); 
+data[100] = 2*y[0]*p[2] - y[2]*(p[2]*cos(y[1] - y[3]) - p[1]*sin(y[1] - y[3])); 
 data[101] = -y[0]*y[2]*(-p[2]*sin(y[1] - y[3]) - p[1]*cos(y[1] - y[3])); 
-data[110] = -y[2]*(-p[2]*sin(y[1] - y[3]) + p[1]*cos(y[1] - y[3])); 
-data[111] = -y[0]*y[2]*(-p[2]*cos(y[1] - y[3]) - p[1]*sin(y[1] - y[3])); 
-data[112] = -y[0]*(-p[2]*sin(y[1] - y[3]) + p[1]*cos(y[1] - y[3])) + 2*y[2]*p[1]; 
-data[113] = -y[0]*y[2]*(p[2]*cos(y[1] - y[3]) + p[1]*sin(y[1] - y[3])); 
-data[122] = -y[2]*(-p[2]*cos(y[1] - y[3]) - p[1]*sin(y[1] - y[3])); 
-data[123] = -y[0]*y[2]*(p[2]*sin(y[1] - y[3]) - p[1]*cos(y[1] - y[3])); 
-data[124] = -y[0]*(-p[2]*cos(y[1] - y[3]) - p[1]*sin(y[1] - y[3])) - 2*y[2]*p[2] - y[2]*p[3]; 
-data[125] = -y[0]*y[2]*(-p[2]*sin(y[1] - y[3]) + p[1]*cos(y[1] - y[3])); 
+data[102] = -y[0]*(p[2]*cos(y[1] - y[3]) - p[1]*sin(y[1] - y[3])); 
+data[103] = -y[0]*y[2]*(p[2]*sin(y[1] - y[3]) + p[1]*cos(y[1] - y[3])); 
+data[114] = -y[2]*(p[2]*sin(y[1] - y[3]) - p[1]*cos(y[1] - y[3])); 
+data[115] = -y[0]*y[2]*(p[2]*cos(y[1] - y[3]) + p[1]*sin(y[1] - y[3])); 
+data[116] = -y[0]*(p[2]*sin(y[1] - y[3]) - p[1]*cos(y[1] - y[3])) - 2*y[2]*p[1]; 
+data[117] = -y[0]*y[2]*(-p[2]*cos(y[1] - y[3]) - p[1]*sin(y[1] - y[3])); 
+data[128] = -y[2]*(p[2]*cos(y[1] - y[3]) + p[1]*sin(y[1] - y[3])); 
+data[129] = -y[0]*y[2]*(-p[2]*sin(y[1] - y[3]) + p[1]*cos(y[1] - y[3])); 
+data[130] = -y[0]*(p[2]*cos(y[1] - y[3]) + p[1]*sin(y[1] - y[3])) + 2*y[2]*p[2]; 
+data[131] = -y[0]*y[2]*(p[2]*sin(y[1] - y[3]) - p[1]*cos(y[1] - y[3])); 
+data[142] = -1.4433756729740647e-6*p[0]*sqrt(pow(y[4], 2) + pow(y[5], 2))/pow(y[0], 2); 
+data[146] = 1.4433756729740647e-6*p[0]*y[4]*pow(pow(y[4], 2) + pow(y[5], 2), -0.5)/y[0]; 
+data[147] = 1.4433756729740647e-6*p[0]*y[5]*pow(pow(y[4], 2) + pow(y[5], 2), -0.5)/y[0]; 
+data[158] = -1.4433756729740647e-6*p[0]*sqrt(pow(y[6], 2) + pow(y[7], 2))/pow(y[2], 2); 
+data[162] = 1.4433756729740647e-6*p[0]*y[6]*pow(pow(y[6], 2) + pow(y[7], 2), -0.5)/y[2]; 
+data[163] = 1.4433756729740647e-6*p[0]*y[7]*pow(pow(y[6], 2) + pow(y[7], 2), -0.5)/y[2]; 
 
 }
 
 void de_jac_trap_num_eval(double *data,double *x,double *y,double *u,double *p,double Dt){
 
-data[26] = 1; 
-data[39] = 1; 
-data[78] = -1; 
-data[91] = 1; 
-data[104] = 1; 
-data[117] = 1; 
-data[130] = 1; 
-data[143] = -1; 
+data[30] = 1; 
+data[45] = 1; 
+data[90] = 1; 
+data[105] = 1; 
+data[120] = 1; 
+data[135] = 1; 
+data[150] = -1; 
+data[165] = -1; 
+data[180] = -1; 
+data[195] = -1; 
 
 }
 

pydae/bmapu/lines/temp_xy_0.json

@@ -3,5 +3,9 @@
     "theta_1": 0.0,
     "V_2": 1.0,
     "theta_2": 0.0,
+    "p_line_pu_1_2": 0.1,
+    "p_line_pu_2_1": -0.1,
+    "I_1_2": 0.1,
+    "I_2_1": -0.1,
     "omega_coi": 1.0
 }

pydae/bmapu/lines/xy_0.json

@@ -3,5 +3,9 @@
     "theta_1": 0.0,
     "V_2": 1.0,
     "theta_2": 0.0,
+    "p_line_pu_1_2": 0.1,
+    "p_line_pu_2_1": -0.1,
+    "I_1_2": 0.1,
+    "I_2_1": -0.1,
     "omega_coi": 1.0
 }

pydae/templates/class_dae_template_v2.py

@@ -263,6 +263,10 @@ class model:
     #     jac_run_ss_eval_xy(self.jac_run,self.x,self.y_run,self.u_run,self.p)
     #     jac_run_ss_eval_up(self.jac_run,self.x,self.y_run,self.u_run,self.p)
         
+        
+    def jac_ini_eval(self):
+        de_jac_ini_eval(self.jac_ini,self.x,self.y_ini,self.u_ini,self.p,self.Dt)
+      
     def jac_run_eval(self):
         de_jac_run_eval(self.jac_run,self.x,self.y_run,self.u_run,self.p,self.Dt)
       

pydae/utils/algebra/algebra.py

@@ -0,0 +1,61 @@
+import numpy as np
+import sympy as sym
+from pydae.bmapu import bmapu_builder
+from pydae.build_v2 import build_mkl,build_numba
+
+
+
+def build():
+    L,G,M,K_d = sym.symbols('L,G,M,K_d', real=True)
+    p_x,p_y,v_x,v_y = sym.symbols('p_x,p_y,v_x,v_y', real=True) 
+    lam,f_x,theta = sym.symbols('lam,f_x,theta', real=True) 
+
+
+    dp_x = v_x
+    dp_y = v_y
+    dv_x = (-2*p_x*lam + f_x - K_d*v_x)/M
+    dv_y = (-M*G - 2*p_y*lam - K_d*v_y)/M   
+
+    g_1 = p_x**2 + p_y**2 - L**2 -lam*1e-6
+    g_2 = -theta + sym.atan2(p_x,-p_y)
+
+    params_dict = {'L':5.21,'G':9.81,'M':10.0,'K_d':1e-3}  # parameters with default values
+
+    u_ini_dict = {'theta':np.deg2rad(5.0)}  # input for the initialization problem
+    u_run_dict = {'f_x':0}                  # input for the running problem, its value is updated 
+
+    sys_dict = {'name':'pendulum',
+                'params_dict':params_dict,
+                'f_list':[dp_x,dp_y,dv_x,dv_y],
+                'g_list':[g_1,g_2],
+                'x_list':[ p_x, p_y, v_x, v_y],
+                'y_ini_list':[lam,f_x],
+                'y_run_list':[lam,theta],
+                'u_ini_dict':u_ini_dict,
+                'u_run_dict':u_run_dict,
+                'h_dict':{'E_p':M*G*(p_y+L),'E_k':0.5*M*(v_x**2+v_y**2),'f_x':f_x}}
+
+    build_numba(sys_dict)
+
+
+
+if __name__ == "__main__":
+
+    build()
+    import pendulum
+    model = pendulum.model()
+
+
+    M = 30.0  # mass of the bob (kg)
+    L = 5.21  # length of the pendulum (m)
+
+    model.ini({'M':M,'L':L,           # parameters setting
+            'theta':np.deg2rad(0)  # initial desired angle = 0º
+          },-1)                  # here -1 means that -1 is considered as initial gess for
+                                 # dynamic and algebraic states
+
+    model.report_x()  # obtained dynamic states
+    model.report_y()  # obtained algebraic states
+    model.report_z()  # obtained outputs
+    model.report_u()  # obtained algebraic states (theta is both state and output; f_x is both input and output)
+    model.report_params()  # considered parameters

pydae/utils/algebra/dev_algebra.ipynb

@@ -0,0 +1,191 @@
+{
+ "cells": [
+  {
+   "cell_type": "code",
+   "execution_count": 1,
+   "id": "3cebf691",
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "text/plain": [
+       "4"
+      ]
+     },
+     "execution_count": 1,
+     "metadata": {},
+     "output_type": "execute_result"
+    }
+   ],
+   "source": [
+    "2+2"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 2,
+   "id": "d0478c7d",
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "p_x   =  0.00\n",
+      "p_y   = -5.21\n",
+      "v_x   =  0.00\n",
+      "v_y   =  0.00\n",
+      "lam   = 28.24\n",
+      "theta =  0.00\n",
+      "E_p   = -0.00\n",
+      "E_k   =  0.00\n",
+      "f_x   = -0.00\n",
+      "f_x   =-0.00\n",
+      "L     = 5.21\n",
+      "G     = 9.81\n",
+      "M     =30.00\n",
+      "K_d   = 0.00\n"
+     ]
+    }
+   ],
+   "source": [
+    "import pendulum\n",
+    "import numpy as np\n",
+    "model = pendulum.model()\n",
+    "\n",
+    "\n",
+    "M = 30.0  # mass of the bob (kg)\n",
+    "L = 5.21  # length of the pendulum (m)\n",
+    "\n",
+    "model.ini({'M':M,'L':L,           # parameters setting\n",
+    "        'theta':np.deg2rad(0)  # initial desired angle = 0º\n",
+    "        },-1)                  # here -1 means that -1 is considered as initial gess for\n",
+    "                                # dynamic and algebraic states\n",
+    "\n",
+    "model.report_x()  # obtained dynamic states\n",
+    "model.report_y()  # obtained algebraic states\n",
+    "model.report_z()  # obtained outputs\n",
+    "model.report_u()  # obtained algebraic states (theta is both state and output; f_x is both input and output)\n",
+    "model.report_params()  # considered parameters\n"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 3,
+   "id": "fc98ab21",
+   "metadata": {},
+   "outputs": [
+    {
+     "ename": "LinAlgError",
+     "evalue": "Matrix is singular to machine precision.",
+     "output_type": "error",
+     "traceback": [
+      "\u001b[31m---------------------------------------------------------------------------\u001b[39m",
+      "\u001b[31mLinAlgError\u001b[39m                               Traceback (most recent call last)",
+      "\u001b[36mCell\u001b[39m\u001b[36m \u001b[39m\u001b[32mIn[3]\u001b[39m\u001b[32m, line 10\u001b[39m\n\u001b[32m      1\u001b[39m xy_0_dict = {\n\u001b[32m      2\u001b[39m     \u001b[33m\"\u001b[39m\u001b[33mp_x\u001b[39m\u001b[33m\"\u001b[39m: \u001b[32m1.943753270464573e-27\u001b[39m,\n\u001b[32m      3\u001b[39m     \u001b[33m\"\u001b[39m\u001b[33mp_y\u001b[39m\u001b[33m\"\u001b[39m: \u001b[32m0\u001b[39m,\n\u001b[32m   (...)\u001b[39m\u001b[32m      7\u001b[39m     \u001b[33m\"\u001b[39m\u001b[33mf_x\u001b[39m\u001b[33m\"\u001b[39m: -\u001b[32m1.0210510686201435e-24\u001b[39m\n\u001b[32m      8\u001b[39m }\n\u001b[32m---> \u001b[39m\u001b[32m10\u001b[39m \u001b[43mmodel\u001b[49m\u001b[43m.\u001b[49m\u001b[43mini\u001b[49m\u001b[43m(\u001b[49m\u001b[43m{\u001b[49m\u001b[33;43m'\u001b[39;49m\u001b[33;43mM\u001b[39;49m\u001b[33;43m'\u001b[39;49m\u001b[43m:\u001b[49m\u001b[43mM\u001b[49m\u001b[43m,\u001b[49m\u001b[33;43m'\u001b[39;49m\u001b[33;43mL\u001b[39;49m\u001b[33;43m'\u001b[39;49m\u001b[43m:\u001b[49m\u001b[43mL\u001b[49m\u001b[43m,\u001b[49m\u001b[43m           \u001b[49m\u001b[38;5;66;43;03m# parameters setting\u001b[39;49;00m\n\u001b[32m     11\u001b[39m \u001b[43m        \u001b[49m\u001b[33;43m'\u001b[39;49m\u001b[33;43mtheta\u001b[39;49m\u001b[33;43m'\u001b[39;49m\u001b[43m:\u001b[49m\u001b[43mnp\u001b[49m\u001b[43m.\u001b[49m\u001b[43mdeg2rad\u001b[49m\u001b[43m(\u001b[49m\u001b[32;43m0\u001b[39;49m\u001b[43m)\u001b[49m\u001b[43m  \u001b[49m\u001b[38;5;66;43;03m# initial desired angle = 0º\u001b[39;49;00m\n\u001b[32m     12\u001b[39m \u001b[43m        \u001b[49m\u001b[43m}\u001b[49m\u001b[43m,\u001b[49m\u001b[43mxy_0_dict\u001b[49m\u001b[43m)\u001b[49m   \n",
+      "\u001b[36mFile \u001b[39m\u001b[32mc:\\Users\\jmmau\\workspace\\pydae\\src\\pydae\\utils\\algebra\\pendulum.py:501\u001b[39m, in \u001b[36mmodel.ini\u001b[39m\u001b[34m(self, up_dict, xy_0)\u001b[39m\n\u001b[32m    498\u001b[39m \u001b[38;5;28;01mif\u001b[39;00m \u001b[38;5;28mtype\u001b[39m(xy_0) == \u001b[38;5;28mfloat\u001b[39m \u001b[38;5;129;01mor\u001b[39;00m \u001b[38;5;28mtype\u001b[39m(xy_0) == \u001b[38;5;28mint\u001b[39m:\n\u001b[32m    499\u001b[39m     \u001b[38;5;28mself\u001b[39m.xy_0 = np.ones(\u001b[38;5;28mself\u001b[39m.N_x+\u001b[38;5;28mself\u001b[39m.N_y,dtype=np.float64)*xy_0\n\u001b[32m--> \u001b[39m\u001b[32m501\u001b[39m xy_ini,it = \u001b[43msstate\u001b[49m\u001b[43m(\u001b[49m\u001b[38;5;28;43mself\u001b[39;49m\u001b[43m.\u001b[49m\u001b[43mxy_0\u001b[49m\u001b[43m,\u001b[49m\u001b[38;5;28;43mself\u001b[39;49m\u001b[43m.\u001b[49m\u001b[43mu_ini\u001b[49m\u001b[43m,\u001b[49m\u001b[38;5;28;43mself\u001b[39;49m\u001b[43m.\u001b[49m\u001b[43mp\u001b[49m\u001b[43m,\u001b[49m\n\u001b[32m    502\u001b[39m \u001b[43m                   \u001b[49m\u001b[38;5;28;43mself\u001b[39;49m\u001b[43m.\u001b[49m\u001b[43mjac_ini\u001b[49m\u001b[43m,\u001b[49m\n\u001b[32m    503\u001b[39m \u001b[43m                   \u001b[49m\u001b[38;5;28;43mself\u001b[39;49m\u001b[43m.\u001b[49m\u001b[43mN_x\u001b[49m\u001b[43m,\u001b[49m\u001b[38;5;28;43mself\u001b[39;49m\u001b[43m.\u001b[49m\u001b[43mN_y\u001b[49m\u001b[43m,\u001b[49m\n\u001b[32m    504\u001b[39m \u001b[43m                   \u001b[49m\u001b[43mmax_it\u001b[49m\u001b[43m=\u001b[49m\u001b[38;5;28;43mself\u001b[39;49m\u001b[43m.\u001b[49m\u001b[43mmax_it\u001b[49m\u001b[43m,\u001b[49m\u001b[43mtol\u001b[49m\u001b[43m=\u001b[49m\u001b[38;5;28;43mself\u001b[39;49m\u001b[43m.\u001b[49m\u001b[43mitol\u001b[49m\u001b[43m)\u001b[49m\n\u001b[32m    506\u001b[39m \u001b[38;5;28;01mif\u001b[39;00m it < \u001b[38;5;28mself\u001b[39m.max_it-\u001b[32m1\u001b[39m:\n\u001b[32m    508\u001b[39m     \u001b[38;5;28mself\u001b[39m.xy_ini = xy_ini\n",
+      "\u001b[36mFile \u001b[39m\u001b[32mc:\\Users\\jmmau\\miniconda3\\Lib\\site-packages\\numba\\np\\linalg.py:913\u001b[39m, in \u001b[36m_inv_err_handler\u001b[39m\u001b[34m()\u001b[39m\n\u001b[32m    911\u001b[39m     \u001b[38;5;28;01massert\u001b[39;00m \u001b[32m0\u001b[39m   \u001b[38;5;66;03m# unreachable\u001b[39;00m\n\u001b[32m    912\u001b[39m \u001b[38;5;28;01mif\u001b[39;00m r > \u001b[32m0\u001b[39m:\n\u001b[32m--> \u001b[39m\u001b[32m913\u001b[39m     \u001b[38;5;28;01mraise\u001b[39;00m np.linalg.LinAlgError(\n\u001b[32m    914\u001b[39m         \u001b[33m\"\u001b[39m\u001b[33mMatrix is singular to machine precision.\u001b[39m\u001b[33m\"\u001b[39m)\n",
+      "\u001b[31mLinAlgError\u001b[39m: Matrix is singular to machine precision."
+     ]
+    }
+   ],
+   "source": [
+    "xy_0_dict = {\n",
+    "    \"p_x\": 1.943753270464573e-27,\n",
+    "    \"p_y\": 0,\n",
+    "    \"v_x\": 0.0,\n",
+    "    \"v_y\": 0.0,\n",
+    "    \"lam\": 0,\n",
+    "    \"f_x\": -1.0210510686201435e-24\n",
+    "}\n",
+    "\n",
+    "model.ini({'M':M,'L':L,           # parameters setting\n",
+    "        'theta':np.deg2rad(0)  # initial desired angle = 0º\n",
+    "        },xy_0_dict)   \n",
+    "\n"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 6,
+   "id": "20241fd3",
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "text/plain": [
+       "array([[ 0.00000000e+00,  0.00000000e+00,  1.00000000e+00,\n",
+       "         0.00000000e+00,  0.00000000e+00,  0.00000000e+00],\n",
+       "       [ 0.00000000e+00,  0.00000000e+00,  0.00000000e+00,\n",
+       "         1.00000000e+00,  0.00000000e+00,  0.00000000e+00],\n",
+       "       [-0.00000000e+00,  0.00000000e+00, -3.33333333e-05,\n",
+       "         0.00000000e+00, -0.00000000e+00,  3.33333333e-02],\n",
+       "       [ 0.00000000e+00, -0.00000000e+00,  0.00000000e+00,\n",
+       "        -3.33333333e-05, -0.00000000e+00,  0.00000000e+00],\n",
+       "       [ 0.00000000e+00,  0.00000000e+00,  0.00000000e+00,\n",
+       "         0.00000000e+00, -1.00000000e-06,  0.00000000e+00],\n",
+       "       [            nan,             nan,  0.00000000e+00,\n",
+       "         0.00000000e+00,  0.00000000e+00,  0.00000000e+00]])"
+      ]
+     },
+     "execution_count": 6,
+     "metadata": {},
+     "output_type": "execute_result"
+    }
+   ],
+   "source": [
+    "xy_0_dict = {\n",
+    "    \"p_x\": 0.0,\n",
+    "    \"p_y\": 0.0,\n",
+    "    \"v_x\": 0.0,\n",
+    "    \"v_y\": 0.0,\n",
+    "    \"lam\": 0.0,\n",
+    "    \"f_x\": 0.0\n",
+    "}\n",
+    "\n",
+    "model.dict2xy0(xy_0_dict)\n",
+    "model.jac_ini_eval()\n",
+    "model.jac_ini"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "758508ea",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "model.save_xy_0('xy_0.json')"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "81147185",
+   "metadata": {},
+   "outputs": [],
+   "source": []
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "base",
+   "language": "python",
+   "name": "python3"
+  },
+  "language_info": {
+   "codemirror_mode": {
+    "name": "ipython",
+    "version": 3
+   },
+   "file_extension": ".py",
+   "mimetype": "text/x-python",
+   "name": "python",
+   "nbconvert_exporter": "python",
+   "pygments_lexer": "ipython3",
+   "version": "3.13.9"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 5
+}