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

author	Sascha L. Teichmann <sascha.teichmann@intevation.de>
date	Sun, 05 Jun 2011 18:24:46 +0000
parents
children	baea7981477a

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
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	3 import java.util.List;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	4 import java.util.Arrays;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	5 import java.util.Comparator;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	6 import java.util.Collections;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	7
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	8 import de.intevation.flys.utils.DoubleUtil;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	9
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	10 import de.intevation.flys.model.River;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	11 import de.intevation.flys.model.Gauge;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	12
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	13 import de.intevation.flys.artifacts.model.WstValueTable.QPosition;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	14
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	15 import de.intevation.flys.artifacts.math.Function;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	16 import de.intevation.flys.artifacts.math.Linear;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	17 import de.intevation.flys.artifacts.math.Identity;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	18 import de.intevation.flys.artifacts.math.BackJumpCorrector;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	19
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	20 import org.apache.log4j.Logger;
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 public class Calculation4
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 private static Logger logger = Logger.getLogger(Calculation4.class);
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 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	27
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	28 public static final Comparator<Segment> REF_CMP =
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	29 new Comparator<Segment>() {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	30 public int compare(Segment a, Segment b) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	31 double d = a.referencePoint - b.referencePoint;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	32 if (d < 0d) return -1;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	33 return d > 0d ? +1 : 0;
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 };
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	36
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	37 protected List<Segment> segments;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	38
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	39 protected boolean isQ;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	40
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	41 public Calculation4() {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	42 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	43
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	44 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	45
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	46 this.segments = segments;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	47 this.isQ = isQ;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	48
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	49 int numResults = 1;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	50
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	51 // assign reference points
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	52 for (Segment segment: segments) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	53 Gauge gauge = river.determineGauge(
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	54 segment.getFrom(), segment.getTo());
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 segment.setReferencePoint(gauge != null
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	57 ? gauge.getStation().doubleValue()
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	58 : 0.5*(segment.getFrom() + segment.getTo()));
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	59
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	60 double [] values = segment.values;
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 if (numResults == -1) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	63 numResults = values.length;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	64 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	65 else if (numResults != values.length) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	66 throw new IllegalArgumentException("wrong length of values");
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	67 }
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 // convert to Q if needed
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	70 if (!isQ && gauge != null) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	71 double [][] table = new DischargeTables(
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	72 river.getName(), gauge.getName()).getFirstTable();
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	73
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	74 // need the original values for naming
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	75 segment.backup();
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	76
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	77 for (int i = 0; i < values.length; ++i) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	78 values[i] = DischargeTables.getQForW(table, values[i]);
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 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	81 } // for all segments
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 Collections.sort(segments, REF_CMP);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	84 }
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 public WQKms [] calculate(
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	87 WstValueTable table,
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	88 double from, double to, double step
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 boolean debug = logger.isDebugEnabled();
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	91
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	92 if (debug) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	93 logger.debug(
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	94 "calculate from " + from + " to " + to + " step " + step);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	95 logger.debug("# segments: " + segments.size());
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
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	98 if (segments.isEmpty()) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	99 logger.debug("no segments found");
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	100 return new WQKms[0];
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	101 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	102
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	103 int numResults = segments.get(0).values.length;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	104
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	105 if (numResults < 1) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	106 logger.debug("no values given");
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	107 return new WQKms[0];
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	108 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	109
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	110
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	111 WQKms [] results = new WQKms[numResults];
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	112 if (debug) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	113 logger.debug("after last segment -> gleichwertig");
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	114 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	115 for (int i = 0; i < results.length; ++i) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	116 results[i] = new WQKms();
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	117 }
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 if (Math.abs(step) < MINIMAL_STEP_WIDTH) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	120 step = MINIMAL_STEP_WIDTH;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	121 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	122
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	123 if (from > to) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	124 step = -step;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	125 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	126
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	127 QPosition [] qPositions = new QPosition[numResults];
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	128
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	129 Function [] functions = new Function[numResults];
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 double [] out = new double[2];
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	132
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	133 Segment sentinel = new Segment(Double.MAX_VALUE);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	134 Segment s1 = sentinel, s2 = sentinel;
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 for (double pos = from;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	137 from < to ? pos <= to : pos >= to;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	138 pos = DoubleUtil.round(pos + step)
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	139 ) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	140 if (pos < s1.referencePoint \|\| pos > s2.referencePoint) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	141 if (debug) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	142 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	143 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	144 // find new interval
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	145 if (pos <= segments.get(0).referencePoint) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	146 // before first segment -> "gleichwertig"
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	147 if (debug) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	148 logger.debug("before first segment -> gleichwertig");
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	149 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	150 Segment first = segments.get(0);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	151 double [] values = first.values;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	152 double refPos = first.referencePoint;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	153 for (int i = 0; i < qPositions.length; ++i) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	154 qPositions[i] = table.getQPosition(
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	155 refPos, values[i]);
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 sentinel.setReferencePoint(-Double.MAX_VALUE);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	158 s1 = sentinel;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	159 s2 = segments.get(0);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	160 Arrays.fill(functions, Identity.IDENTITY);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	161 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	162 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	163 // after last segment -> "gleichwertig"
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	164 if (debug) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	165 logger.debug("after last segment -> gleichwertig");
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	166 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	167 Segment last = segments.get(segments.size()-1);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	168 double [] values = last.values;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	169 double refPos = last.referencePoint;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	170 for (int i = 0; i < qPositions.length; ++i) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	171 qPositions[i] = table.getQPosition(
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	172 refPos, values[i]);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	173 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	174 sentinel.setReferencePoint(Double.MAX_VALUE);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	175 s1 = last;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	176 s2 = sentinel;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	177 Arrays.fill(functions, Identity.IDENTITY);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	178 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	179 else { // "ungleichwertig"
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	180 // find matching interval
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	181 if (debug) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	182 logger.debug("in segments -> ungleichwertig");
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	183 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	184 s1 = s2 = null;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	185 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	186 Segment si1 = segments.get(i-1);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	187 Segment si = segments.get(i);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	188 if (debug) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	189 logger.debug("check " + pos + " in " +
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	190 si1.referencePoint + " - " + si.referencePoint);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	191 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	192 if (pos >= si1.referencePoint
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	193 && pos <= si. referencePoint) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	194 s1 = si1;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	195 s2 = si;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	196 break;
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 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	199
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	200 if (s1 == null) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	201 throw new IllegalStateException("no interval found");
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	202 }
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 Segment anchor, free;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	205
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	206 if (from > to) { anchor = s1; free = s2; }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	207 else { anchor = s2; free = s1; }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	208
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	209 // build transforms based on "gleichwertiger" phase
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	210 for (int i = 0; i < qPositions.length; ++i) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	211 QPosition qi = table.getQPosition(
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	212 anchor.referencePoint,
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	213 anchor.values[i]);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	214
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	215 if ((qPositions[i] = qi) == null) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	216 // TODO: error report
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	217 functions[i] = Identity.IDENTITY;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	218 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	219 else {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	220 double qA = table.getQ(qi, anchor.referencePoint);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	221 double qF = table.getQ(qi, free .referencePoint);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	222
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	223 functions[i] = Double.isNaN(qA) \|\| Double.isNaN(qF)
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	224 ? Identity.IDENTITY
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	225 : new Linear(
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	226 qA, qF,
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	227 anchor.values[i], free.values[i]);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	228 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	229 } // build transforms
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	230 } // "ungleichwertiges" interval
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	231 } // find matching interval
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 for (int i = 0; i < qPositions.length; ++i) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	234 QPosition qPosition = qPositions[i];
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	235
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	236 if (qPosition == null) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	237 continue;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	238 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	239
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	240 if (table.interpolate(pos, out, qPosition, functions[i])) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	241 results[i].add(out[0], out[1], pos);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	242 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	243 else {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	244 // TODO: error report
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	245 }
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 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	248
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	249 // Backjump correction
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	250 for (int i = 0; i < results.length; ++i) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	251 BackJumpCorrector bjc = new BackJumpCorrector();
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	252
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	253 double [] ws = results[i].getWs();
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	254 double [] kms = results[i].getKms();
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	255
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	256 if (bjc.doCorrection(kms, ws)) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	257 results[i] = new WQCKms(results[i], bjc.getCorrected());
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	258 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	259 // TODO: error report
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	260 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	261
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	262 // name the curves
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	263 for (int i = 0; i < results.length; ++i) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	264 results[i].setName(createName(i));
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	265 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	266
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	267 return results;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	268 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	269
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	270 protected String createName(int index) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	271 // TODO: i18n
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	272 StringBuilder sb = new StringBuilder(isQ ? "Q" : "W");
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	273 sb.append(" benutzerdefiniert (");
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	274 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	275 if (i > 0) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	276 sb.append("; ");
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	277 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	278 Segment segment = segments.get(i);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	279 sb.append((segment.backup != null
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	280 ? segment.backup
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	281 : segment.values)[index]);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	282 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	283 sb.append(')');
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	284 return sb.toString();
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	285 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	286 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	287 // vim:set ts=4 sw=4 si et sta sts=4 fenc=utf8 :

655

913b52064449 Refactored version of "Berechnung 4"