FilesExpand file tree

 master

/

NumInt.java

Copy path

BlameMore file actions

BlameMore file actions

Latest commit

 

History

History

History

54 lines (49 loc) · 1.33 KB

 master

/

NumInt.java

Top

File metadata and controls

• Code
• Blame

54 lines (49 loc) · 1.33 KB

Raw

Copy raw file

Download raw file

Open symbols panel

Edit and raw actions

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30

31

32

33

34

35

36

37

38

39

40

41

42

43

44

45

46

47

48

49

50

51

52

53

54

package symjava.numeric;

import symjava.bytecode.BytecodeFunc;

import symjava.symbolic.Func;

import symjava.symbolic.Integrate;

/**

* Numerical Integration

*

*/

public class NumInt {

Integrate integrate;

BytecodeFunc byteFunc;

public NumInt(Integrate integrate) {

this.integrate = integrate;

if(this.integrate.integrand instanceof Func) {

Func f = (Func) this.integrate.integrand;

byteFunc = f.toBytecodeFunc();

} else {

Func f = new Func(

this.getClass().getSimpleName()+java.util.UUID.randomUUID().toString().replaceAll("-", ""),

this.integrate.integrand);

byteFunc = f.toBytecodeFunc();

}

}

public double eval() {

double[][] pnts = integrate.domain.getIntWeightAndPoints(3);

int dim = integrate.domain.getDim();

double sum = 0.0;

for(int k=0; k<pnts.length; k++) {

sum += byteFunc.apply(pnts[k])*pnts[k][dim];

}

return sum;

}

public double eval(double ...params) {

double[][] pnts = integrate.domain.getIntWeightAndPoints(3);

int dim = integrate.domain.getDim();

double[] allArgs = new double[dim + params.length];

double sum = 0.0;

for(int k=0; k<pnts.length; k++) {

for(int i=0; i<dim; i++) {

allArgs[i] = pnts[k][i];

}

for(int i=dim; i<allArgs.length; i++) {

allArgs[i] = params[i-dim];

}

sum += byteFunc.apply(allArgs)*pnts[k][dim];

}

return sum;

}

}