dive4elements/river: flys-artifacts/src/main/java/de/intevation/flys/artifacts/model/fixings/FixCalculation.java annotate

annotate flys-artifacts/src/main/java/de/intevation/flys/artifacts/model/fixings/FixCalculation.java @ 2993:70469e3d34b9

FixA: Made Delta W/t calculation finally work flys-artifacts/trunk@4545 c6561f87-3c4e-4783-a992-168aeb5c3f6f

author	Sascha L. Teichmann <sascha.teichmann@intevation.de>
date	Mon, 28 May 2012 20:59:27 +0000
parents	0abdede5a0b8
children	a2b265b9b2f7

rev	line source
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	1 package de.intevation.flys.artifacts.model.fixings;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	2
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	3 import de.intevation.flys.artifacts.FixationArtifactAccess;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	4
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	5 import de.intevation.flys.artifacts.math.fitting.Function;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	6 import de.intevation.flys.artifacts.math.fitting.FunctionFactory;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	7
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	8 import de.intevation.flys.artifacts.model.Calculation;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	9 import de.intevation.flys.artifacts.model.CalculationResult;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	10 import de.intevation.flys.artifacts.model.FixingsColumn;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	11 import de.intevation.flys.artifacts.model.FixingsColumnFactory;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	12 import de.intevation.flys.artifacts.model.FixingsOverview;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	13 import de.intevation.flys.artifacts.model.FixingsOverviewFactory;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	14 import de.intevation.flys.artifacts.model.Parameters;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	15
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	16 import de.intevation.flys.artifacts.model.FixingsOverview.Range;
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	17 import de.intevation.flys.artifacts.model.FixingsOverview.Fixing;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	18 import de.intevation.flys.artifacts.model.FixingsOverview.IdFilter;
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	19 import de.intevation.flys.artifacts.model.FixingsOverview.AndFilter;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	20 import de.intevation.flys.artifacts.model.FixingsOverview.DateRangeFilter;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	21 import de.intevation.flys.artifacts.model.FixingsOverview.SectorRangeFilter;
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	22
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	23 import de.intevation.flys.utils.DoubleUtil;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	24
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	25 import java.util.ArrayList;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	26 import java.util.List;
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	27 import java.util.Date;
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	28
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	29 import org.apache.commons.math.MathException;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	30
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	31 import org.apache.commons.math.optimization.fitting.CurveFitter;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	32
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	33 import org.apache.commons.math.optimization.general.LevenbergMarquardtOptimizer;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	34
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	35 import org.apache.log4j.Logger;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	36
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	37 public class FixCalculation
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	38 extends Calculation
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	39 {
2786 f9e62e76ee1f FixA: Return correct data structure from calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2783 diff changeset	40 private static Logger log = Logger.getLogger(FixCalculation.class);
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	41
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	42 protected String river;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	43 protected double from;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	44 protected double to;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	45 protected double step;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	46 protected boolean preprocessing;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	47 protected String function;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	48 protected int [] events;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	49 protected long [][] analysisPeriods;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	50 protected int qSectorStart;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	51 protected int qSectorEnd;
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	52
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	53 public FixCalculation() {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	54 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	55
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	56 public FixCalculation(FixationArtifactAccess access) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	57
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	58 String river = access.getRiver();
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	59 Double from = access.getFrom();
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	60 Double to = access.getTo();
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	61 Double step = access.getStep();
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	62 String function = access.getFunction();
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	63 int [] events = access.getEvents();
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	64 long [][] analysisPeriods = access.getAnalysisPeriods();
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	65 Integer qSectorStart = access.getQSectorStart();
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	66 Integer qSectorEnd = access.getQSectorEnd();
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	67
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	68 if (river == null) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	69 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	70 addProblem("fix.missing.river");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	71 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	72
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	73 if (from == null) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	74 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	75 addProblem("fix.missing.from");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	76 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	77
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	78 if (to == null) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	79 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	80 addProblem("fix.missing.to");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	81 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	82
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	83 if (step == null) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	84 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	85 addProblem("fix.missing.step");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	86 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	87
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	88 if (function == null) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	89 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	90 addProblem("fix.missing.function");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	91 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	92
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	93 if (events == null \|\| events.length < 1) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	94 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	95 addProblem("fix.missing.events");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	96 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	97
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	98 if (analysisPeriods == null \|\| analysisPeriods.length < 1) {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	99 // TODO: i18n
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	100 addProblem("fix.missing.analysis.periods");
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	101 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	102
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	103 if (qSectorStart == null) {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	104 // TODO: i18n
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	105 addProblem("fix.missing.qstart.sector");
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	106 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	107
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	108 if (qSectorEnd == null) {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	109 // TODO: i18n
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	110 addProblem("fix.missing.qend.sector");
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	111 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	112
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	113 if (!hasProblems()) {
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	114 this.river = river;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	115 this.from = from;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	116 this.to = to;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	117 this.step = step;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	118 this.function = function;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	119 this.events = events;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	120 this.analysisPeriods = analysisPeriods;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	121 this.qSectorStart = qSectorStart;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	122 this.qSectorEnd = qSectorEnd;
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	123 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	124 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	125
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	126 public CalculationResult calculate() {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	127
2992 0abdede5a0b8 FixA: Added facet to return delta w/t as CSV Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2788 diff changeset	128 boolean debug = log.isDebugEnabled();
0abdede5a0b8 FixA: Added facet to return delta w/t as CSV Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2788 diff changeset	129
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	130 FixingsOverview overview =
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	131 FixingsOverviewFactory.getOverview(river);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	132
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	133 if (overview == null) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	134 addProblem("fix.no.overview.available");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	135 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	136
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	137 Function func = FunctionFactory.getInstance()
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	138 .getFunction(function);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	139
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	140 if (func == null) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	141 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	142 addProblem("fix.invalid.function.name");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	143 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	144
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	145 if (hasProblems()) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	146 return new CalculationResult(this);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	147 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	148
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	149 FixingsColumnFactory fcf = FixingsColumnFactory.getInstance();
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	150
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	151 List<FixingsColumn> dataColumns =
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	152 new ArrayList<FixingsColumn>(events.length);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	153
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	154 for (int eventId: events) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	155 IdFilter idFilter = new IdFilter(eventId);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	156
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	157 List<Fixing.Column> columns = overview.filter(null, idFilter);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	158 if (columns.isEmpty()) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	159 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	160 addProblem("fix.missing.column", eventId);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	161 continue;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	162 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	163 FixingsColumn dataColumn = fcf.getColumnData(columns.get(0));
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	164 if (dataColumn == null) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	165 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	166 addProblem("fix.cannot.load.data", eventId);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	167 continue;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	168 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	169 dataColumns.add(dataColumn);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	170 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	171
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	172 if (dataColumns.size() < 2) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	173 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	174 addProblem("fix.too.less.data.columns");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	175 return new CalculationResult(this);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	176 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	177
2992 0abdede5a0b8 FixA: Added facet to return delta w/t as CSV Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2788 diff changeset	178 double [] kms = DoubleUtil.explode(from, to, step / 1000.0);
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	179
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	180 double [] ws = new double[dataColumns.size()];
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	181 double [] qs = new double[ws.length];
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	182
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	183 String [] parameterNames = func.getParameterNames();
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	184
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	185 Parameters results =
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	186 new Parameters(createColumnNames(parameterNames));
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	187
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	188 boolean invalid = false;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	189
2992 0abdede5a0b8 FixA: Added facet to return delta w/t as CSV Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2788 diff changeset	190 if (debug) {
0abdede5a0b8 FixA: Added facet to return delta w/t as CSV Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2788 diff changeset	191 log.debug("number of kms: " + kms.length);
0abdede5a0b8 FixA: Added facet to return delta w/t as CSV Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2788 diff changeset	192 }
0abdede5a0b8 FixA: Added facet to return delta w/t as CSV Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2788 diff changeset	193
0abdede5a0b8 FixA: Added facet to return delta w/t as CSV Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2788 diff changeset	194 int numFailed = 0;
0abdede5a0b8 FixA: Added facet to return delta w/t as CSV Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2788 diff changeset	195
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	196 for (int i = 0; i < kms.length; ++i) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	197 double km = kms[i];
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	198
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	199 for (int j = 0; j < ws.length; ++j) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	200 FixingsColumn column = dataColumns.get(j);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	201 qs[j] = column.getQ(km);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	202 boolean interpolated = column.getW(km, ws, j);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	203 // TODO: mark as interpolated.
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	204 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	205
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	206 // TODO: Do preprocessing here!
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	207 double [] parameters = fit(func, km, ws, qs);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	208 if (parameters == null) { // Problems are reported already.
2992 0abdede5a0b8 FixA: Added facet to return delta w/t as CSV Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2788 diff changeset	209 ++numFailed;
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	210 continue;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	211 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	212
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	213 int row = results.newRow();
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	214
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	215 results.set(row, "km", km);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	216 for (int j = 0; j < parameters.length; ++j) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	217 if (Double.isNaN(parameters[j])) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	218 invalid = true;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	219 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	220 else {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	221 results.set(row, parameterNames[j], parameters[j]);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	222 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	223 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	224 // TODO: Calculate statistics, too!
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	225 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	226
2992 0abdede5a0b8 FixA: Added facet to return delta w/t as CSV Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2788 diff changeset	227 if (debug) {
0abdede5a0b8 FixA: Added facet to return delta w/t as CSV Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2788 diff changeset	228 log.debug("success: " + (kms.length - numFailed));
0abdede5a0b8 FixA: Added facet to return delta w/t as CSV Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2788 diff changeset	229 log.debug("failed: " + numFailed);
0abdede5a0b8 FixA: Added facet to return delta w/t as CSV Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2788 diff changeset	230 }
0abdede5a0b8 FixA: Added facet to return delta w/t as CSV Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2788 diff changeset	231
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	232 if (invalid) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	233 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	234 addProblem("fix.invalid.values");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	235 results.removeNaNs();
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	236 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	237
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	238 // Calculate Delta W/t
2783 cbb513a8f548 FixA: km indexed datastructure for delta w/t results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2744 diff changeset	239 DeltaWTsKM deltaWTsKM = calculateDeltaWTs(
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	240 func,
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	241 overview,
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	242 results);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	243
2786 f9e62e76ee1f FixA: Return correct data structure from calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2783 diff changeset	244 FixResult fr = new FixResult(results, deltaWTsKM);
f9e62e76ee1f FixA: Return correct data structure from calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2783 diff changeset	245
f9e62e76ee1f FixA: Return correct data structure from calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2783 diff changeset	246 return new CalculationResult(fr, this);
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	247 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	248
2783 cbb513a8f548 FixA: km indexed datastructure for delta w/t results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2744 diff changeset	249 public DeltaWTsKM calculateDeltaWTs(
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	250 Function function,
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	251 FixingsOverview overview,
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	252 Parameters results
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	253 ) {
2993 70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	254 boolean debug = log.isDebugEnabled();
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	255
2783 cbb513a8f548 FixA: km indexed datastructure for delta w/t results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2744 diff changeset	256 DeltaWTsKM deltaWTsKM = new DeltaWTsKM(results.size());
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	257
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	258 Column [][] analysisColumns = getAnalysisColumns(overview);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	259
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	260 int [] parameterIndices =
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	261 results.columnIndices(function.getParameterNames());
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	262
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	263 double [] parameterValues =
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	264 new double[parameterIndices.length];
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	265
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	266 double [] ow = new double[1];
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	267
2993 70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	268 int kmIdx = results.columnIndex("km");
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	269
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	270 for (int i = 0, N = results.size(); i < N; ++i) {
2993 70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	271 double km = results.get(i, kmIdx);
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	272 results.get(i, parameterIndices, parameterValues);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	273
2783 cbb513a8f548 FixA: km indexed datastructure for delta w/t results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2744 diff changeset	274 DeltaWTsKM.KM dwtkm = new DeltaWTsKM.KM(km);
cbb513a8f548 FixA: km indexed datastructure for delta w/t results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2744 diff changeset	275 deltaWTsKM.add(dwtkm);
cbb513a8f548 FixA: km indexed datastructure for delta w/t results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2744 diff changeset	276
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	277 // This is the paraterized function for a given km.
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	278 de.intevation.flys.artifacts.math.Function instance =
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	279 function.instantiate(parameterValues);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	280
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	281 // Evaluate all columns for all analysis periods.
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	282 for (int j = 0; j < analysisColumns.length; ++j) {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	283 Column [] periodColumns = analysisColumns[j];
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	284
2993 70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	285 int failedQ = 0;
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	286 int failedW = 0;
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	287 int failedC = 0;
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	288
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	289 for (int k = 0; k < periodColumns.length; ++k) {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	290 Column pc = periodColumns[k];
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	291
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	292 // Q from real data.
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	293 double q = pc.data.getQ(km);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	294 if (Double.isNaN(q)) {
2993 70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	295 ++failedQ;
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	296 continue;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	297 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	298
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	299 // Calculate W from function.
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	300 double nw = instance.value(q);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	301 if (Double.isNaN(nw)) {
2993 70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	302 ++failedC;
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	303 continue;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	304 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	305
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	306 // W from real data.
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	307 pc.data.getW(km, ow);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	308
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	309 if (Double.isNaN(ow[0])) {
2993 70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	310 ++failedW;
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	311 continue;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	312 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	313
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	314 double deltaW = ow[0] - nw;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	315
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	316 DeltaWT deltaWT = new DeltaWT(
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	317 deltaW,
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	318 pc.meta.getStartTime(),
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	319 pc.meta.getDescription());
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	320
2783 cbb513a8f548 FixA: km indexed datastructure for delta w/t results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2744 diff changeset	321 dwtkm.add(deltaWT);
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	322 }
2993 70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	323 if (debug) {
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	324 log.debug("failed W: " + failedW);
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	325 log.debug("failed Q: " + failedQ);
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	326 log.debug("failed C: " + failedC);
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	327 log.debug("input size: " + periodColumns.length);
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	328 log.debug("outpt size: " + dwtkm.size());
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	329 }
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	330 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	331 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	332
2783 cbb513a8f548 FixA: km indexed datastructure for delta w/t results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2744 diff changeset	333 return deltaWTsKM;
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	334 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	335
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	336 /** Helper class to bundle the meta information of a column
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	337 * and the real data.
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	338 */
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	339 protected static class Column {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	340
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	341 protected Fixing.Column meta;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	342 protected FixingsColumn data;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	343
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	344 public Column() {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	345 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	346
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	347 public Column(Fixing.Column meta, FixingsColumn data) {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	348 this.meta = meta;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	349 this.data = data;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	350 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	351 } // class Column
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	352
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	353 /** Fetch meta and data columns for analysis periods. */
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	354 protected Column [][] getAnalysisColumns(FixingsOverview overview) {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	355
2993 70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	356 boolean debug = log.isDebugEnabled();
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	357 if (debug) {
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	358 log.debug("number analysis periods: " + analysisPeriods.length);
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	359 }
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	360
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	361 Column columns [][] = new Column[analysisPeriods.length][];
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	362
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	363 Range range = new Range(from, to);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	364 SectorRangeFilter sectorRangeFilter =
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	365 new SectorRangeFilter(qSectorStart, qSectorEnd);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	366
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	367 FixingsColumnFactory fcf = FixingsColumnFactory.getInstance();
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	368
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	369 for (int i = 0; i < columns.length; ++i) {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	370
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	371 // Construct filter for period.
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	372 long [] period = analysisPeriods[i];
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	373
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	374 AndFilter filter = new AndFilter();
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	375
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	376 DateRangeFilter dateRangeFilter =
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	377 new DateRangeFilter(
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	378 new Date(period[0]),
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	379 new Date(period[1]));
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	380
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	381 filter.add(dateRangeFilter);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	382 filter.add(sectorRangeFilter);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	383
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	384 List<Fixing.Column> metaCols = overview.filter(range, filter);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	385
2993 70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	386 if (debug) {
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	387 log.debug("number of filtered columns: " + metaCols.size());
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	388 }
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	389
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	390 ArrayList<Column> cols = new ArrayList<Column>(metaCols.size());
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	391
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	392 // Only use columns which have data.
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	393 for (Fixing.Column meta: metaCols) {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	394 FixingsColumn data = fcf.getColumnData(meta);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	395 if (data != null) {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	396 cols.add(new Column(meta, data));
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	397 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	398 }
2993 70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	399
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	400 if (debug) {
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	401 log.debug("failed loading: " + (metaCols.size()-cols.size()));
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	402 }
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	403 columns[i] = cols.toArray(new Column[cols.size()]);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	404 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	405
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	406 return columns;
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	407 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	408
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	409 protected static String [] createColumnNames(String [] parameters) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	410 String [] result = new String[parameters.length + 1];
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	411 result[0] = "km";
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	412 // TODO: Add statistic columns, too.
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	413 System.arraycopy(parameters, 0, result, 1, parameters.length);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	414 return result;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	415 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	416
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	417 protected double [] fit(
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	418 Function function,
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	419 double km,
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	420 double [] ws,
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	421 double [] qs
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	422 ) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	423 LevenbergMarquardtOptimizer lmo = new LevenbergMarquardtOptimizer();
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	424 CurveFitter cf = new CurveFitter(lmo);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	425
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	426 boolean missingWs = false;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	427 boolean missingQs = false;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	428
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	429 for (int i = 0; i < ws.length; ++i) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	430 boolean ignore = false;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	431 if (Double.isNaN(ws[i])) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	432 ignore = true;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	433 if (!missingWs) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	434 missingWs = true;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	435 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	436 addProblem(km, "fix.missing.w");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	437 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	438 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	439 if (Double.isNaN(qs[i])) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	440 ignore = true;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	441 if (!missingQs) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	442 missingQs = true;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	443 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	444 addProblem(km, "fix.missing.q");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	445 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	446 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	447 if (!ignore) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	448 cf.addObservedPoint(ws[i], qs[i]);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	449 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	450 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	451
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	452 try {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	453 return cf.fit(function, function.getInitialGuess());
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	454 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	455 catch (MathException me) {
2992 0abdede5a0b8 FixA: Added facet to return delta w/t as CSV Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2788 diff changeset	456 log.warn(me, me);
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	457 addProblem(km, "fix.fitting.failed");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	458 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	459
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	460 return null;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	461 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	462 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	463 // vim:set ts=4 sw=4 si et sta sts=4 fenc=utf8 :

Mercurial > dive4elements > river

annotate flys-artifacts/src/main/java/de/intevation/flys/artifacts/model/fixings/FixCalculation.java @ 2993:70469e3d34b9