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

annotate flys-artifacts/src/main/java/de/intevation/flys/artifacts/model/Calculation4.java @ 4187:21f4e4b79121

Refactor GaugeDischargeCurveFacet to be able to set a facet name For adding another output of the GaugeDischargeCurveArtifact it is necessary to provide to facet instances with different names. Therefore the GaugeDischargeCurveFacet is extended to set the facet name in the constructor.

author	Björn Ricks <bjoern.ricks@intevation.de>
date	Fri, 19 Oct 2012 13:25:49 +0200
parents	22790758b132
children	39885bdfc6fc

rev	line source
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	1 package de.intevation.flys.artifacts.model;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	2
3441 bfbd478bd607 FixA: Moved Q -> W conversion code into segments to be reusable. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3439 diff changeset	3 import de.intevation.flys.artifacts.math.BackJumpCorrector;
bfbd478bd607 FixA: Moved Q -> W conversion code into segments to be reusable. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3439 diff changeset	4 import de.intevation.flys.artifacts.math.Function;
bfbd478bd607 FixA: Moved Q -> W conversion code into segments to be reusable. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3439 diff changeset	5 import de.intevation.flys.artifacts.math.Identity;
bfbd478bd607 FixA: Moved Q -> W conversion code into segments to be reusable. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3439 diff changeset	6 import de.intevation.flys.artifacts.math.Linear;
bfbd478bd607 FixA: Moved Q -> W conversion code into segments to be reusable. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3439 diff changeset	7
bfbd478bd607 FixA: Moved Q -> W conversion code into segments to be reusable. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3439 diff changeset	8 import de.intevation.flys.artifacts.model.WstValueTable.QPosition;
bfbd478bd607 FixA: Moved Q -> W conversion code into segments to be reusable. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3439 diff changeset	9
bfbd478bd607 FixA: Moved Q -> W conversion code into segments to be reusable. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3439 diff changeset	10 import de.intevation.flys.model.River;
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	11
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	12 import de.intevation.flys.utils.DoubleUtil;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	13
3441 bfbd478bd607 FixA: Moved Q -> W conversion code into segments to be reusable. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3439 diff changeset	14 import java.util.Arrays;
bfbd478bd607 FixA: Moved Q -> W conversion code into segments to be reusable. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3439 diff changeset	15 import java.util.List;
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	16
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	17 import org.apache.log4j.Logger;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	18
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	19 public class Calculation4
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 656 diff changeset	20 extends Calculation
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	21 {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	22 private static Logger logger = Logger.getLogger(Calculation4.class);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	23
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	24 public static final double MINIMAL_STEP_WIDTH = 1e-5;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	25
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	26 protected List<Segment> segments;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	27
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	28 protected boolean isQ;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	29
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	30 public Calculation4() {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	31 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	32
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	33 public Calculation4(List<Segment> segments, River river, boolean isQ) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	34
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	35 this.segments = segments;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	36 this.isQ = isQ;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	37
3441 bfbd478bd607 FixA: Moved Q -> W conversion code into segments to be reusable. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3439 diff changeset	38 Segment.setReferencePointConvertQ(segments, river, isQ, this);
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	39 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	40
709 3b7e9ddf6bb1 New model to transport data and error reports of calculations. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 704 diff changeset	41 public CalculationResult calculate(
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	42 WstValueTable table,
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	43 double from, double to, double step
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	44 ) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	45 boolean debug = logger.isDebugEnabled();
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	46
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	47 if (debug) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	48 logger.debug(
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	49 "calculate from " + from + " to " + to + " step " + step);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	50 logger.debug("# segments: " + segments.size());
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 656 diff changeset	51 for (Segment segment: segments) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 656 diff changeset	52 logger.debug(" " + segment);
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 656 diff changeset	53 }
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	54 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	55
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	56 if (segments.isEmpty()) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	57 logger.debug("no segments found");
2166 2898b1ff6013 I18N for WINFO calculation result messages. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 742 diff changeset	58 addProblem("no.segments.found");
709 3b7e9ddf6bb1 New model to transport data and error reports of calculations. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 704 diff changeset	59 return new CalculationResult(new WQKms[0], this);
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	60 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	61
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	62 int numResults = segments.get(0).values.length;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	63
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	64 if (numResults < 1) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	65 logger.debug("no values given");
2166 2898b1ff6013 I18N for WINFO calculation result messages. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 742 diff changeset	66 addProblem("no.values.given");
709 3b7e9ddf6bb1 New model to transport data and error reports of calculations. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 704 diff changeset	67 return new CalculationResult(new WQKms[0], this);
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	68 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	69
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	70
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	71 WQKms [] results = new WQKms[numResults];
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	72 for (int i = 0; i < results.length; ++i) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	73 results[i] = new WQKms();
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	74 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	75
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	76 if (Math.abs(step) < MINIMAL_STEP_WIDTH) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	77 step = MINIMAL_STEP_WIDTH;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	78 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	79
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	80 if (from > to) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	81 step = -step;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	82 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	83
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	84 QPosition [] qPositions = new QPosition[numResults];
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	85
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	86 Function [] functions = new Function[numResults];
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	87
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	88 double [] out = new double[2];
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	89
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	90 Segment sentinel = new Segment(Double.MAX_VALUE);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	91 Segment s1 = sentinel, s2 = sentinel;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	92
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	93 for (double pos = from;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	94 from < to ? pos <= to : pos >= to;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	95 pos = DoubleUtil.round(pos + step)
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	96 ) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	97 if (pos < s1.referencePoint \|\| pos > s2.referencePoint) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	98 if (debug) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	99 logger.debug("need to find new interval for " + pos);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	100 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	101 // find new interval
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	102 if (pos <= segments.get(0).referencePoint) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	103 // before first segment -> "gleichwertig"
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	104 if (debug) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	105 logger.debug("before first segment -> gleichwertig");
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	106 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	107 Segment first = segments.get(0);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	108 double [] values = first.values;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	109 double refPos = first.referencePoint;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	110 for (int i = 0; i < qPositions.length; ++i) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	111 qPositions[i] = table.getQPosition(
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	112 refPos, values[i]);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	113 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	114 sentinel.setReferencePoint(-Double.MAX_VALUE);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	115 s1 = sentinel;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	116 s2 = segments.get(0);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	117 Arrays.fill(functions, Identity.IDENTITY);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	118 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	119 else if (pos >= segments.get(segments.size()-1).referencePoint) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	120 // after last segment -> "gleichwertig"
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	121 if (debug) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	122 logger.debug("after last segment -> gleichwertig");
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	123 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	124 Segment last = segments.get(segments.size()-1);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	125 double [] values = last.values;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	126 double refPos = last.referencePoint;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	127 for (int i = 0; i < qPositions.length; ++i) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	128 qPositions[i] = table.getQPosition(
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	129 refPos, values[i]);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	130 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	131 sentinel.setReferencePoint(Double.MAX_VALUE);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	132 s1 = last;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	133 s2 = sentinel;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	134 Arrays.fill(functions, Identity.IDENTITY);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	135 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	136 else { // "ungleichwertig"
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	137 // find matching interval
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	138 if (debug) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	139 logger.debug("in segments -> ungleichwertig");
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	140 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	141 s1 = s2 = null;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	142 for (int i = 1, N = segments.size(); i < N; ++i) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	143 Segment si1 = segments.get(i-1);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	144 Segment si = segments.get(i);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	145 if (debug) {
742 c09c9e05ecfa Removed trailing whitespace. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 739 diff changeset	146 logger.debug("check " + pos + " in " +
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	147 si1.referencePoint + " - " + si.referencePoint);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	148 }
742 c09c9e05ecfa Removed trailing whitespace. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 739 diff changeset	149 if (pos >= si1.referencePoint
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	150 && pos <= si. referencePoint) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	151 s1 = si1;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	152 s2 = si;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	153 break;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	154 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	155 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	156
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	157 if (s1 == null) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	158 throw new IllegalStateException("no interval found");
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	159 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	160
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	161 Segment anchor, free;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	162
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	163 if (from > to) { anchor = s1; free = s2; }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	164 else { anchor = s2; free = s1; }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	165
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	166 // build transforms based on "gleichwertiger" phase
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	167 for (int i = 0; i < qPositions.length; ++i) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	168 QPosition qi = table.getQPosition(
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	169 anchor.referencePoint,
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	170 anchor.values[i]);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	171
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	172 if ((qPositions[i] = qi) == null) {
2166 2898b1ff6013 I18N for WINFO calculation result messages. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 742 diff changeset	173 addProblem(pos, "cannot.find.q", anchor.values[i]);
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	174 functions[i] = Identity.IDENTITY;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	175 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	176 else {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	177 double qA = table.getQ(qi, anchor.referencePoint);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	178 double qF = table.getQ(qi, free .referencePoint);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	179
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	180 functions[i] = Double.isNaN(qA) \|\| Double.isNaN(qF)
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	181 ? Identity.IDENTITY
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	182 : new Linear(
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	183 qA, qF,
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	184 anchor.values[i], free.values[i]);
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 656 diff changeset	185
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 656 diff changeset	186 if (debug) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 656 diff changeset	187 logger.debug(
742 c09c9e05ecfa Removed trailing whitespace. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 739 diff changeset	188 anchor.referencePoint + ": " +
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 656 diff changeset	189 qA + " -> " + functions[i].value(qA) +
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 656 diff changeset	190 " / " + free.referencePoint + ": " +
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 656 diff changeset	191 qF + " -> " + functions[i].value(qF));
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 656 diff changeset	192 }
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	193 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	194 } // build transforms
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	195 } // "ungleichwertiges" interval
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	196 } // find matching interval
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	197
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	198 for (int i = 0; i < qPositions.length; ++i) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	199 QPosition qPosition = qPositions[i];
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	200
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	201 if (qPosition == null) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	202 continue;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	203 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	204
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	205 if (table.interpolate(pos, out, qPosition, functions[i])) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	206 results[i].add(out[0], out[1], pos);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	207 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	208 else {
2166 2898b1ff6013 I18N for WINFO calculation result messages. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 742 diff changeset	209 addProblem(pos, "cannot.interpolate.w.q");
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	210 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	211 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	212 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	213
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	214 // Backjump correction
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	215 for (int i = 0; i < results.length; ++i) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	216 BackJumpCorrector bjc = new BackJumpCorrector();
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	217
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	218 double [] ws = results[i].getWs();
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	219 double [] kms = results[i].getKms();
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	220
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 656 diff changeset	221 if (bjc.doCorrection(kms, ws, this)) {
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	222 results[i] = new WQCKms(results[i], bjc.getCorrected());
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	223 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	224 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	225
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	226 // name the curves
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	227 for (int i = 0; i < results.length; ++i) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	228 results[i].setName(createName(i));
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	229 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	230
709 3b7e9ddf6bb1 New model to transport data and error reports of calculations. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 704 diff changeset	231 return new CalculationResult(results, this);
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	232 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	233
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	234 protected String createName(int index) {
3450 22790758b132 FixA/Vollmer: return the result now. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 3441 diff changeset	235 // TODO: i18n
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	236 StringBuilder sb = new StringBuilder(isQ ? "Q" : "W");
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	237 sb.append(" benutzerdefiniert (");
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	238 for (int i = 0, N = segments.size(); i < N; ++i) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	239 if (i > 0) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	240 sb.append("; ");
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	241 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	242 Segment segment = segments.get(i);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	243 sb.append((segment.backup != null
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	244 ? segment.backup
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	245 : segment.values)[index]);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	246 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	247 sb.append(')');
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	248 return sb.toString();
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	249 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	250 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	251 // 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/Calculation4.java @ 4187:21f4e4b79121