Changeset 16227

Timestamp:

10/10/18 12:00:19 (6 years ago)

Author:

gkronber

Message:

#2925: changed test case to linear oscillator

File:

• branches/2925_AutoDiffForDynamicalModels/AutoDiffForDynamicalModelsTest/TestCvodes.cs (modified) (3 diffs)

branches/2925_AutoDiffForDynamicalModels/AutoDiffForDynamicalModelsTest/TestCvodes.cs

@@ -30 +30 @@
 
 
-      // simple temperature annealing
-      int numberOfEquations = 1;
+      // linear oscillator
+      int numberOfEquations = 2;
       var y = Problem.N_VNew_Serial(numberOfEquations);
       // y must be initialized before calling CVodeInit
       // Problem.N_VConst_Serial(100.0, y);
-      Problem.NV_Set_Ith_S(y, 0, 100.0);
+      Problem.NV_Set_Ith_S(y, 0, 0.5); // x
+      Problem.NV_Set_Ith_S(y, 1, 1);  // v
 
       var cvode_mem = Problem.CVodeCreate(Problem.MultistepMethod.CV_ADAMS, Problem.NonlinearSolverIteration.CV_FUNCTIONAL);
@@ -65 +66 @@
         flag = Problem.CVode(cvode_mem, tout, y, ref t, Problem.CV_NORMAL);
         Assert.AreEqual(Problem.CV_SUCCESS, flag);
-        Console.WriteLine("t={0} y={1}", t, Problem.N_VL1Norm_Serial(y));
+        Console.WriteLine("{0} {1} {2}", t, Problem.NV_Get_Ith_S(y, 0), Problem.NV_Get_Ith_S(y, 1));
         // Problem.N_VPrint_Serial(y);
         tout += 0.1; 
@@ -87 +88 @@
       IntPtr user_data) {
 
-      Problem.NV_Set_Ith_S(ydot, 0, -0.1 * Problem.NV_Get_Ith_S(y, 0));
+      Problem.NV_Set_Ith_S(ydot, 0, Problem.NV_Get_Ith_S(y, 1));
+      Problem.NV_Set_Ith_S(ydot, 1, -0.3 * Problem.NV_Get_Ith_S(y, 0));
       return 0;
       ;