svn commit: r780510 - in /commons/proper/math/trunk/src/test/org/apache/commons/math/ode/nonstiff: TestProblem6.java TestProblemFactory.java

[lu...@apache.org]

Author: luc
Date: Sun May 31 21:56:52 2009
New Revision: 780510

URL: http://svn.apache.org/viewvc?rev=780510&view=rev
Log:
added a new test problem with polynomial behavior

Added:
    commons/proper/math/trunk/src/test/org/apache/commons/math/ode/nonstiff/TestProblem6.java   (with props)
Modified:
    commons/proper/math/trunk/src/test/org/apache/commons/math/ode/nonstiff/TestProblemFactory.java

Added: commons/proper/math/trunk/src/test/org/apache/commons/math/ode/nonstiff/TestProblem6.java
URL: http://svn.apache.org/viewvc/commons/proper/math/trunk/src/test/org/apache/commons/math/ode/nonstiff/TestProblem6.java?rev=780510&view=auto
==============================================================================
--- commons/proper/math/trunk/src/test/org/apache/commons/math/ode/nonstiff/TestProblem6.java (added)
+++ commons/proper/math/trunk/src/test/org/apache/commons/math/ode/nonstiff/TestProblem6.java Sun May 31 21:56:52 2009
@@ -0,0 +1,96 @@
+/*
+ * Licensed to the Apache Software Foundation (ASF) under one or more
+ * contributor license agreements.  See the NOTICE file distributed with
+ * this work for additional information regarding copyright ownership.
+ * The ASF licenses this file to You under the Apache License, Version 2.0
+ * (the "License"); you may not use this file except in compliance with
+ * the License.  You may obtain a copy of the License at
+ *
+ *      http://www.apache.org/licenses/LICENSE-2.0
+ *
+ * Unless required by applicable law or agreed to in writing, software
+ * distributed under the License is distributed on an "AS IS" BASIS,
+ * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+ * See the License for the specific language governing permissions and
+ * limitations under the License.
+ */
+
+package org.apache.commons.math.ode.nonstiff;
+
+/**
+ * This class is used in the junit tests for the ODE integrators.
+
+ * <p>This specific problem is the following differential equation :
+ * <pre>
+ *    y' = 3x^5 - y
+ * </pre>
+ * when the initial condition is y(0) = -360, the solution of this
+ * equation degenerates to a simple quintic polynomial function :
+ * <pre>
+ *   y (t) = 3x^5 - 15x^4 + 60x^3 - 180x^2 + 360x - 360
+ * </pre>
+ * </p>
+
+ */
+public class TestProblem6
+  extends TestProblemAbstract {
+
+    /** Serializable version identifier. */
+    private static final long serialVersionUID = 1353409119804352378L;
+
+    /** theoretical state */
+    private double[] y;
+
+    /**
+     * Simple constructor.
+     */
+    public TestProblem6() {
+        super();
+        double[] y0 = { -360.0 };
+        setInitialConditions(0.0, y0);
+        setFinalConditions(1.0);
+        double[] errorScale = { 1.0 };
+        setErrorScale(errorScale);
+        y = new double[y0.length];
+    }
+
+    /**
+     * Copy constructor.
+     * @param problem problem to copy
+     */
+    public TestProblem6(TestProblem6 problem) {
+        super(problem);
+        y = problem.y.clone();
+    }
+
+    /**
+     * Clone operation.
+     * @return a copy of the instance
+     */
+    @Override
+    public Object clone() {
+        return new TestProblem6(this);
+    }
+
+    @Override
+    public void doComputeDerivatives(double t, double[] y, double[] yDot) {
+
+        // compute the derivatives
+        double t2 = t  * t;
+        double t4 = t2 * t2;
+        double t5 = t4 * t;
+        for (int i = 0; i < n; ++i) {
+            yDot[i] = 3 * t5 - y[i];
+        }
+
+    }
+
+    @Override
+    public double[] computeTheoreticalState(double t) {
+        for (int i = 0; i < n; ++i) {
+            y[i] = ((((3 * t - 15) * t + 60) * t - 180) * t + 360) * t - 360;
+        }
+        return y;
+    }
+
+}

Propchange: commons/proper/math/trunk/src/test/org/apache/commons/math/ode/nonstiff/TestProblem6.java
------------------------------------------------------------------------------
    svn:eol-style = native

Propchange: commons/proper/math/trunk/src/test/org/apache/commons/math/ode/nonstiff/TestProblem6.java
------------------------------------------------------------------------------
    svn:keywords = Author Date Id Revision

Modified: commons/proper/math/trunk/src/test/org/apache/commons/math/ode/nonstiff/TestProblemFactory.java
URL: http://svn.apache.org/viewvc/commons/proper/math/trunk/src/test/org/apache/commons/math/ode/nonstiff/TestProblemFactory.java?rev=780510&r1=780509&r2=780510&view=diff
==============================================================================
--- commons/proper/math/trunk/src/test/org/apache/commons/math/ode/nonstiff/TestProblemFactory.java (original)
+++ commons/proper/math/trunk/src/test/org/apache/commons/math/ode/nonstiff/TestProblemFactory.java Sun May 31 21:56:52 2009
@@ -28,7 +28,8 @@
     new TestProblem2(),
     new TestProblem3(),
     new TestProblem4(),
-    new TestProblem5()
+    new TestProblem5(),
+    new TestProblem6()
   };
 
   /**