Skip to content

Navigation Menu

Sign in
Appearance settings

Search code, repositories, users, issues, pull requests...

Provide feedback

We read every piece of feedback, and take your input very seriously.

Saved searches

Use saved searches to filter your results more quickly

Appearance settings

Commit f36eb2b

Browse filesBrowse files
committed
add scipy optimize examples
1 parent 3a40398 commit f36eb2b
Copy full SHA for f36eb2b

File tree

Expand file treeCollapse file tree

1 file changed

+36
-0
lines changed
Open diff view settings
Filter options
Expand file treeCollapse file tree

1 file changed

+36
-0
lines changed
Open diff view settings
Collapse file
+36Lines changed: 36 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,36 @@
1+
package symjava.examples;
2+
3+
import symjava.relational.Eq;
4+
import symjava.symbolic.Expr;
5+
import symjava.symbolic.Sum;
6+
import symjava.symbolic.Symbols;
7+
import symjava.symbolic.utils.AddList;
8+
import static symjava.math.SymMath.*;
9+
import static symjava.symbolic.Symbol.C0;
10+
11+
public class OptimizationExamples {
12+
13+
public static void main(String[] args) {
14+
test1();
15+
}
16+
17+
public static void test1() {
18+
//Rosenbrock function of N variables
19+
Symbols x = new Symbols("x");
20+
AddList sum = new AddList();
21+
int N = 5;
22+
for(int i=2; i<=N; i++) {
23+
sum.add( 100*pow(x[i]-x[i-1]*x[i-1],2) + pow(1-x[i-1],2) );
24+
}
25+
double[] x0 = {1.3, 0.7, 0.8, 1.9, 1.2};
26+
27+
Expr[] freeVars = x.getRange(1, N);
28+
Eq eq = new Eq(sum.toExpr(), C0, freeVars);
29+
System.out.println(eq);
30+
31+
NewtonOptimization.solve(eq, x0, 1000, 1e-6, false);
32+
//[ 1. 1. 1. 1. 1.]
33+
34+
}
35+
36+
}

0 commit comments

Comments
0 (0)
Morty Proxy This is a proxified and sanitized view of the page, visit original site.