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 @ 3003:e5e6363e6cba

FixA: Explicit model for date ranges. flys-artifacts/trunk@4559 c6561f87-3c4e-4783-a992-168aeb5c3f6f

author	Sascha L. Teichmann <sascha.teichmann@intevation.de>
date	Wed, 30 May 2012 14:45:36 +0000
parents	a2b265b9b2f7
children	9e0500d64524

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;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	27
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	28 import org.apache.commons.math.MathException;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	29
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	30 import org.apache.commons.math.optimization.fitting.CurveFitter;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	31
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	32 import org.apache.commons.math.optimization.general.LevenbergMarquardtOptimizer;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	33
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	34 import org.apache.log4j.Logger;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	35
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	36 public class FixCalculation
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	37 extends Calculation
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	38 {
2786 f9e62e76ee1f FixA: Return correct data structure from calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2783 diff changeset	39 private static Logger log = Logger.getLogger(FixCalculation.class);
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	40
3003 e5e6363e6cba FixA: Explicit model for date ranges. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2994 diff changeset	41 protected String river;
e5e6363e6cba FixA: Explicit model for date ranges. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2994 diff changeset	42 protected double from;
e5e6363e6cba FixA: Explicit model for date ranges. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2994 diff changeset	43 protected double to;
e5e6363e6cba FixA: Explicit model for date ranges. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2994 diff changeset	44 protected double step;
e5e6363e6cba FixA: Explicit model for date ranges. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2994 diff changeset	45 protected boolean preprocessing;
e5e6363e6cba FixA: Explicit model for date ranges. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2994 diff changeset	46 protected String function;
e5e6363e6cba FixA: Explicit model for date ranges. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2994 diff changeset	47 protected int [] events;
e5e6363e6cba FixA: Explicit model for date ranges. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2994 diff changeset	48 protected DateRange [] analysisPeriods;
e5e6363e6cba FixA: Explicit model for date ranges. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2994 diff changeset	49 protected int qSectorStart;
e5e6363e6cba FixA: Explicit model for date ranges. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2994 diff changeset	50 protected int qSectorEnd;
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	51
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	52 public FixCalculation() {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	53 }
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 public FixCalculation(FixationArtifactAccess access) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	56
3003 e5e6363e6cba FixA: Explicit model for date ranges. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2994 diff changeset	57 String river = access.getRiver();
e5e6363e6cba FixA: Explicit model for date ranges. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2994 diff changeset	58 Double from = access.getFrom();
e5e6363e6cba FixA: Explicit model for date ranges. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2994 diff changeset	59 Double to = access.getTo();
e5e6363e6cba FixA: Explicit model for date ranges. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2994 diff changeset	60 Double step = access.getStep();
e5e6363e6cba FixA: Explicit model for date ranges. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2994 diff changeset	61 String function = access.getFunction();
e5e6363e6cba FixA: Explicit model for date ranges. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2994 diff changeset	62 int [] events = access.getEvents();
e5e6363e6cba FixA: Explicit model for date ranges. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2994 diff changeset	63 DateRange [] analysisPeriods = access.getAnalysisPeriods();
e5e6363e6cba FixA: Explicit model for date ranges. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2994 diff changeset	64 Integer qSectorStart = access.getQSectorStart();
e5e6363e6cba FixA: Explicit model for date ranges. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2994 diff changeset	65 Integer qSectorEnd = access.getQSectorEnd();
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	66
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	67 if (river == null) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	68 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	69 addProblem("fix.missing.river");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	70 }
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 if (from == null) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	73 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	74 addProblem("fix.missing.from");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	75 }
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 if (to == null) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	78 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	79 addProblem("fix.missing.to");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	80 }
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 if (step == null) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	83 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	84 addProblem("fix.missing.step");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	85 }
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 if (function == null) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	88 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	89 addProblem("fix.missing.function");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	90 }
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 if (events == null \|\| events.length < 1) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	93 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	94 addProblem("fix.missing.events");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	95 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	96
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	97 if (analysisPeriods == null \|\| analysisPeriods.length < 1) {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	98 // TODO: i18n
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	99 addProblem("fix.missing.analysis.periods");
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	100 }
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 if (qSectorStart == null) {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	103 // TODO: i18n
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	104 addProblem("fix.missing.qstart.sector");
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	105 }
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 if (qSectorEnd == null) {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	108 // TODO: i18n
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	109 addProblem("fix.missing.qend.sector");
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	110 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	111
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	112 if (!hasProblems()) {
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	113 this.river = river;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	114 this.from = from;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	115 this.to = to;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	116 this.step = step;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	117 this.function = function;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	118 this.events = events;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	119 this.analysisPeriods = analysisPeriods;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	120 this.qSectorStart = qSectorStart;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	121 this.qSectorEnd = qSectorEnd;
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	122 }
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 public CalculationResult calculate() {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	126
2992 0abdede5a0b8 FixA: Added facet to return delta w/t as CSV Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2788 diff changeset	127 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	128
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	129 FixingsOverview overview =
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	130 FixingsOverviewFactory.getOverview(river);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	131
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	132 if (overview == null) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	133 addProblem("fix.no.overview.available");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	134 }
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 Function func = FunctionFactory.getInstance()
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	137 .getFunction(function);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	138
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	139 if (func == null) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	140 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	141 addProblem("fix.invalid.function.name");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	142 }
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 if (hasProblems()) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	145 return new CalculationResult(this);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	146 }
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 FixingsColumnFactory fcf = FixingsColumnFactory.getInstance();
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	149
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	150 List<FixingsColumn> dataColumns =
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	151 new ArrayList<FixingsColumn>(events.length);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	152
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	153 for (int eventId: events) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	154 IdFilter idFilter = new IdFilter(eventId);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	155
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	156 List<Fixing.Column> columns = overview.filter(null, idFilter);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	157 if (columns.isEmpty()) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	158 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	159 addProblem("fix.missing.column", eventId);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	160 continue;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	161 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	162 FixingsColumn dataColumn = fcf.getColumnData(columns.get(0));
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	163 if (dataColumn == null) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	164 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	165 addProblem("fix.cannot.load.data", eventId);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	166 continue;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	167 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	168 dataColumns.add(dataColumn);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	169 }
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 if (dataColumns.size() < 2) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	172 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	173 addProblem("fix.too.less.data.columns");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	174 return new CalculationResult(this);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	175 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	176
2992 0abdede5a0b8 FixA: Added facet to return delta w/t as CSV Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2788 diff changeset	177 double [] kms = DoubleUtil.explode(from, to, step / 1000.0);
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	178
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	179 double [] ws = new double[dataColumns.size()];
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	180 double [] qs = new double[ws.length];
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	181
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	182 String [] parameterNames = func.getParameterNames();
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	183
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	184 Parameters results =
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	185 new Parameters(createColumnNames(parameterNames));
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	186
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	187 boolean invalid = false;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	188
2992 0abdede5a0b8 FixA: Added facet to return delta w/t as CSV Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2788 diff changeset	189 if (debug) {
0abdede5a0b8 FixA: Added facet to return delta w/t as CSV Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2788 diff changeset	190 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	191 }
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 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	194
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	195 for (int i = 0; i < kms.length; ++i) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	196 double km = kms[i];
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	197
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	198 for (int j = 0; j < ws.length; ++j) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	199 FixingsColumn column = dataColumns.get(j);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	200 qs[j] = column.getQ(km);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	201 boolean interpolated = column.getW(km, ws, j);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	202 // TODO: mark as interpolated.
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	203 }
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 // TODO: Do preprocessing here!
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	206 double [] parameters = fit(func, km, ws, qs);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	207 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	208 ++numFailed;
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	209 continue;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	210 }
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 int row = results.newRow();
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	213
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	214 results.set(row, "km", km);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	215 for (int j = 0; j < parameters.length; ++j) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	216 if (Double.isNaN(parameters[j])) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	217 invalid = true;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	218 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	219 else {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	220 results.set(row, parameterNames[j], parameters[j]);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	221 }
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 // TODO: Calculate statistics, too!
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	224 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	225
2992 0abdede5a0b8 FixA: Added facet to return delta w/t as CSV Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2788 diff changeset	226 if (debug) {
0abdede5a0b8 FixA: Added facet to return delta w/t as CSV Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2788 diff changeset	227 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	228 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	229 }
0abdede5a0b8 FixA: Added facet to return delta w/t as CSV Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2788 diff changeset	230
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	231 if (invalid) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	232 // TODO: i18n
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	233 addProblem("fix.invalid.values");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	234 results.removeNaNs();
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	235 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	236
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	237 // 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	238 DeltaWTsKM deltaWTsKM = calculateDeltaWTs(
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	239 func,
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	240 overview,
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	241 results);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	242
2786 f9e62e76ee1f FixA: Return correct data structure from calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2783 diff changeset	243 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	244
f9e62e76ee1f FixA: Return correct data structure from calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2783 diff changeset	245 return new CalculationResult(fr, this);
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	246 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	247
2783 cbb513a8f548 FixA: km indexed datastructure for delta w/t results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2744 diff changeset	248 public DeltaWTsKM calculateDeltaWTs(
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	249 Function function,
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	250 FixingsOverview overview,
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	251 Parameters results
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	252 ) {
2993 70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	253 boolean debug = log.isDebugEnabled();
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	254
2783 cbb513a8f548 FixA: km indexed datastructure for delta w/t results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2744 diff changeset	255 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	256
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	257 Column [][] analysisColumns = getAnalysisColumns(overview);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	258
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	259 int [] parameterIndices =
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	260 results.columnIndices(function.getParameterNames());
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	261
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	262 double [] parameterValues =
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	263 new double[parameterIndices.length];
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	264
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	265 double [] ow = new double[1];
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	266
2993 70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	267 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	268
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	269 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	270 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	271 results.get(i, parameterIndices, parameterValues);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	272
2783 cbb513a8f548 FixA: km indexed datastructure for delta w/t results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2744 diff changeset	273 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	274 deltaWTsKM.add(dwtkm);
cbb513a8f548 FixA: km indexed datastructure for delta w/t results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2744 diff changeset	275
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	276 // 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	277 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	278 function.instantiate(parameterValues);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	279
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	280 // 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	281 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	282 Column [] periodColumns = analysisColumns[j];
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	283
2993 70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	284 int failedQ = 0;
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	285 int failedW = 0;
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	286 int failedC = 0;
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	287
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	288 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	289 Column pc = periodColumns[k];
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	290
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	291 // Q from real data.
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	292 double q = pc.data.getQ(km);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	293 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	294 ++failedQ;
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	295 continue;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	296 }
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 // Calculate W from function.
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	299 double nw = instance.value(q);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	300 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	301 ++failedC;
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	302 continue;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	303 }
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 // W from real data.
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	306 pc.data.getW(km, ow);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	307
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	308 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	309 ++failedW;
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	310 continue;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	311 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	312
2994 a2b265b9b2f7 FixA: corrected fitting (Q->W instead W->Q). Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2993 diff changeset	313 double deltaW = (ow[0] - nw)*100.0; // in cm
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	314
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	315 DeltaWT deltaWT = new DeltaWT(
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	316 deltaW,
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	317 pc.meta.getStartTime(),
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	318 pc.meta.getDescription());
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	319
2783 cbb513a8f548 FixA: km indexed datastructure for delta w/t results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2744 diff changeset	320 dwtkm.add(deltaWT);
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	321 }
2993 70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	322 if (debug) {
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	323 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	324 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	325 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	326 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	327 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	328 }
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	329 }
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
2783 cbb513a8f548 FixA: km indexed datastructure for delta w/t results. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2744 diff changeset	332 return deltaWTsKM;
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	333 }
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 /** 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	336 * and the real data.
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	337 */
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	338 protected static class Column {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	339
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	340 protected Fixing.Column meta;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	341 protected FixingsColumn data;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	342
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	343 public Column() {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	344 }
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 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	347 this.meta = meta;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	348 this.data = data;
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	349 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	350 } // class Column
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	351
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	352 /** 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	353 protected Column [][] getAnalysisColumns(FixingsOverview overview) {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	354
2993 70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	355 boolean debug = log.isDebugEnabled();
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	356 if (debug) {
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	357 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	358 }
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	359
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	360 Column columns [][] = new Column[analysisPeriods.length][];
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	361
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	362 Range range = new Range(from, to);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	363 SectorRangeFilter sectorRangeFilter =
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	364 new SectorRangeFilter(qSectorStart, qSectorEnd);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	365
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	366 FixingsColumnFactory fcf = FixingsColumnFactory.getInstance();
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	367
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	368 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	369
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	370 // Construct filter for period.
3003 e5e6363e6cba FixA: Explicit model for date ranges. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2994 diff changeset	371 DateRange period = analysisPeriods[i];
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	372
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	373 AndFilter filter = new AndFilter();
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	374
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	375 DateRangeFilter dateRangeFilter =
3003 e5e6363e6cba FixA: Explicit model for date ranges. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2994 diff changeset	376 new DateRangeFilter(period.getFrom(), period.getTo());
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	377
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	378 filter.add(dateRangeFilter);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	379 filter.add(sectorRangeFilter);
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 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	382
2993 70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	383 if (debug) {
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	384 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	385 }
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	386
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	387 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	388
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	389 // Only use columns which have data.
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	390 for (Fixing.Column meta: metaCols) {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	391 FixingsColumn data = fcf.getColumnData(meta);
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	392 if (data != null) {
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	393 cols.add(new Column(meta, data));
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	394 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	395 }
2993 70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	396
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	397 if (debug) {
70469e3d34b9 FixA: Made Delta W/t calculation finally work Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2992 diff changeset	398 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	399 }
2744 c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	400 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	401 }
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	402
c1f2e792704a FixA: Calculate Delta W/t, too. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2729 diff changeset	403 return columns;
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	404 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	405
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	406 protected static String [] createColumnNames(String [] parameters) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	407 String [] result = new String[parameters.length + 1];
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	408 result[0] = "km";
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	409 // TODO: Add statistic columns, too.
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	410 System.arraycopy(parameters, 0, result, 1, parameters.length);
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	411 return result;
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	412 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	413
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	414 protected double [] fit(
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	415 Function function,
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	416 double km,
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	417 double [] ws,
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	418 double [] qs
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	419 ) {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	420 LevenbergMarquardtOptimizer lmo = new LevenbergMarquardtOptimizer();
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	421 CurveFitter cf = new CurveFitter(lmo);
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 for (int i = 0; i < ws.length; ++i) {
2994 a2b265b9b2f7 FixA: corrected fitting (Q->W instead W->Q). Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2993 diff changeset	424 if (!Double.isNaN(ws[i]) && !Double.isNaN(qs[i])) {
a2b265b9b2f7 FixA: corrected fitting (Q->W instead W->Q). Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 2993 diff changeset	425 cf.addObservedPoint(qs[i], ws[i]);
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	426 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	427 }
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 try {
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	430 return cf.fit(function, function.getInitialGuess());
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	431 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	432 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	433 log.warn(me, me);
2729 a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	434 addProblem(km, "fix.fitting.failed");
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	435 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	436
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	437 return null;
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 }
a441be7f1589 Added Fix calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	440 // 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 @ 3003:e5e6363e6cba