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 @ 4176:0e5544c54be8

Fixed flys/issue929. Made zero baseline in historical discharge curve charts visible.

author	Ingo Weinzierl <ingo.weinzierl@intevation.de>
date	Thu, 18 Oct 2012 15:14:53 +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 @ 4176:0e5544c54be8