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;

}

}