private void fitglobal() {
int nparams = 11;
int nsel = 0;
for (int i = 0; i < ncurves; i++) {
if (include[i]) {
nsel++;
}
}
double[][] params = new double[nsel][nparams];
String[][] tempformulas = new String[nsel][nparams];
double[][][] constraints = new double[2][nsel][nparams];
int[][] vflmatrix = new int[nsel][nparams];
int counter = 0;
for (int i = 0; i < ncurves; i++) {
if (include[i]) {
for (int j = 0; j < nparams; j++) {
params[counter][j] = globalparams[i][j];
tempformulas[counter][j] = globalformulas[i][j];
constraints[0][counter][j] = globalconstraints[0][i][j];
constraints[1][counter][j] = globalconstraints[1][i][j];
vflmatrix[counter][j] = globalvflmatrix[i][j];
}
counter++;
}
for (int j = 0; j < nparams; j++) {
undoparams[i][j] = globalparams[i][j];
undoformulas[i][j] = globalformulas[i][j];
undovflmatrix[i][j] = globalvflmatrix[i][j];
}
for (int j = 0; j < xpts; j++) {
for (int k = 0; k < ypts; k++) {
undofit[i][j][k] = fit[i][j][k];
}
}
undoc2[i] = c2[i];
}
undoglobalc2 = globalc2;
if (showglobalfitdialog(params, tempformulas, vflmatrix)) {
counter = 0;
for (int i = 0; i < ncurves; i++) {
if (include[i]) {
for (int j = 0; j < nparams; j++) {
globalparams[i][j] = params[counter][j];
globalformulas[i][j] = tempformulas[counter][j];
globalvflmatrix[i][j] = vflmatrix[counter][j];
}
counter++;
}
}
double[] stats = new double[2];
float[][] tempdata = new float[nsel][xpts * ypts];
float[][] tempweights = new float[nsel][xpts * ypts];
counter = 0;
for (int i = 0; i < ncurves; i++) {
if (include[i]) {
for (int j = 0; j < xpts; j++) {
for (int k = 0; k < ypts; k++) {
tempdata[counter][j + k * xpts] = (float) ((double) pch[i][j][k] / (double) nmeas[i]);
tempweights[counter][j + k * xpts] = weights[i][j][k];
}
}
counter++;
}
}
int tempmaxiter = globalfitclass.maxiter;
if (checkc2) {
globalfitclass.changemaxiter(0);
}
double[] tempc2vals = new double[nsel];
IJ.showStatus("Fitting Globally");
float[][] tempfit =
globalfitclass.fitdata(
params,
vflmatrix,
tempformulas,
paramsnames,
constraints,
tempdata,
tempweights,
stats,
tempc2vals,
false);
IJ.showStatus("Fit Complete");
globalfitclass.changemaxiter(tempmaxiter);
globalc2 = stats[1];
globalc2label.setText("Global chi^2 = " + (float) globalc2);
counter = 0;
for (int i = 0; i < ncurves; i++) {
if (include[i]) {
for (int j = 0; j < xpts; j++) {
for (int k = 0; k < ypts; k++) {
fit[i][j][k] = tempfit[counter][j + xpts * k] * (float) nmeas[i];
}
}
if (i == dispcurve) {
pwfit.updateSeries(fit[dispcurve], 1, true);
}
for (int j = 0; j < nparams; j++) {
globalparams[i][j] = params[counter][j];
}
c2[i] = tempc2vals[counter];
c2array[i].setText("" + (float) c2[i]);
counter++;
}
}
float[] temp = pwfit.getLimits();
temp[4] = 1.0f;
pwfit.setLimits(temp);
}
}
public void showresults(String results) {
IJ.showStatus(results);
IJ.log(results);
}
