View Javadoc
1   /*
2    * Licensed to the Apache Software Foundation (ASF) under one or more
3    * contributor license agreements.  See the NOTICE file distributed with
4    * this work for additional information regarding copyright ownership.
5    * The ASF licenses this file to You under the Apache License, Version 2.0
6    * (the "License"); you may not use this file except in compliance with
7    * the License.  You may obtain a copy of the License at
8    *
9    *      http://www.apache.org/licenses/LICENSE-2.0
10   *
11   * Unless required by applicable law or agreed to in writing, software
12   * distributed under the License is distributed on an "AS IS" BASIS,
13   * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
14   * See the License for the specific language governing permissions and
15   * limitations under the License.
16   */
17  
18  package org.apache.commons.math4.legacy.ode.nonstiff;
19  
20  import org.apache.commons.math4.legacy.ode.FirstOrderDifferentialEquations;
21  import org.apache.commons.math4.legacy.ode.events.EventHandler;
22  
23  
24  public class StepProblem
25    implements FirstOrderDifferentialEquations, EventHandler {
26  
27    public StepProblem(double rateBefore, double rateAfter,
28                       double switchTime) {
29      this.rateAfter  = rateAfter;
30      this.switchTime = switchTime;
31      setRate(rateBefore);
32    }
33  
34    @Override
35  public void computeDerivatives(double t, double[] y, double[] yDot) {
36      yDot[0] = rate;
37    }
38  
39    @Override
40  public int getDimension() {
41      return 1;
42    }
43  
44    public void setRate(double rate) {
45      this.rate = rate;
46    }
47  
48    @Override
49  public void init(double t0, double[] y0, double t) {
50    }
51  
52    @Override
53  public Action eventOccurred(double t, double[] y, boolean increasing) {
54      setRate(rateAfter);
55      return Action.RESET_DERIVATIVES;
56    }
57  
58    @Override
59  public double g(double t, double[] y) {
60      return t - switchTime;
61    }
62  
63    @Override
64  public void resetState(double t, double[] y) {
65    }
66  
67    private double rate;
68    private double rateAfter;
69    private double switchTime;
70  }