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

annotate flys-artifacts/src/main/java/de/intevation/flys/artifacts/model/WstValueTable.java @ 2072:4cdd9c4896f6

#393 Added a new Renderer and Options in Themes that allow displaying minimum and maximum of a chart series. flys-artifacts/trunk@3581 c6561f87-3c4e-4783-a992-168aeb5c3f6f

author	Ingo Weinzierl <ingo.weinzierl@intevation.de>
date	Wed, 04 Jan 2012 12:24:35 +0000
parents	2730d17df021
children	2898b1ff6013

rev	line source
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	1 package de.intevation.flys.artifacts.model;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	2
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	3 import java.io.Serializable;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	4
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	5 import de.intevation.flys.artifacts.math.Linear;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	6 import de.intevation.flys.artifacts.math.Function;
427 909196be11a0 First step to calculate "W fuer ungleichwertige Abfluesse" correctly. flys/issue55 Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 426 diff changeset	7
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	8 import java.util.Arrays;
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	9 import java.util.ArrayList;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	10 import java.util.List;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	11 import java.util.Collections;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	12
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	13 import org.apache.log4j.Logger;
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	14
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	15 import org.apache.commons.math.analysis.interpolation.SplineInterpolator;
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	16
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	17 import org.apache.commons.math.analysis.polynomials.PolynomialSplineFunction;
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	18
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	19 import org.apache.commons.math.ArgumentOutsideDomainException;
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	20
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	21 import org.apache.commons.math.exception.MathIllegalArgumentException;
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	22
1937 f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	23 import gnu.trove.TDoubleArrayList;
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	24
1838 853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	25 /**
853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	26 * W, Q and km data from database 'wsts' spiced with interpolation algorithms.
853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	27 */
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	28 public class WstValueTable
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	29 implements Serializable
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	30 {
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	31 private static Logger log = Logger.getLogger(WstValueTable.class);
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	32
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	33 public static final int DEFAULT_Q_STEPS = 500;
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	34
1721 eb35570df0e8 Cosmetics, docs. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1681 diff changeset	35 /**
eb35570df0e8 Cosmetics, docs. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1681 diff changeset	36 * A Column in the table, typically representing one measurement session.
eb35570df0e8 Cosmetics, docs. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1681 diff changeset	37 */
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	38 public static final class Column
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	39 implements Serializable
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	40 {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	41 protected String name;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	42
632 07640ab913fd First part of storing qs in ranges Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 621 diff changeset	43 protected QRangeTree qRangeTree;
07640ab913fd First part of storing qs in ranges Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 621 diff changeset	44
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	45 public Column() {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	46 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	47
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	48 public Column(String name) {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	49 this.name = name;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	50 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	51
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	52 public String getName() {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	53 return name;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	54 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	55
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	56 public void setName(String name) {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	57 this.name = name;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	58 }
632 07640ab913fd First part of storing qs in ranges Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 621 diff changeset	59
1907 6c7fdfd262ac Do not re-evaluate constant size() in for-loops. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1899 diff changeset	60 public QRangeTree getQRangeTree() {
632 07640ab913fd First part of storing qs in ranges Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 621 diff changeset	61 return qRangeTree;
07640ab913fd First part of storing qs in ranges Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 621 diff changeset	62 }
07640ab913fd First part of storing qs in ranges Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 621 diff changeset	63
07640ab913fd First part of storing qs in ranges Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 621 diff changeset	64 public void setQRangeTree(QRangeTree qRangeTree) {
07640ab913fd First part of storing qs in ranges Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 621 diff changeset	65 this.qRangeTree = qRangeTree;
07640ab913fd First part of storing qs in ranges Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 621 diff changeset	66 }
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	67 } // class Column
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	68
1721 eb35570df0e8 Cosmetics, docs. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1681 diff changeset	69 /**
eb35570df0e8 Cosmetics, docs. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1681 diff changeset	70 * A (weighted) position used for interpolation.
eb35570df0e8 Cosmetics, docs. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1681 diff changeset	71 */
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	72 public static final class QPosition {
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	73
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	74 protected int index;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	75 protected double weight;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	76
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	77 public QPosition() {
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	78 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	79
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	80 public QPosition(int index, double weight) {
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	81 this.index = index;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	82 this.weight = weight;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	83 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	84
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	85 public QPosition set(int index, double weight) {
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	86 this.index = index;
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	87 this.weight = weight;
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	88 return this;
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	89 }
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	90
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	91 } // class Position
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	92
1939 2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	93 public static final class SplineFunction {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	94
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	95 public PolynomialSplineFunction spline;
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	96 public double [] splineQs;
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	97 public double [] splineWs;
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	98
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	99 public SplineFunction(
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	100 PolynomialSplineFunction spline,
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	101 double [] splineQs,
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	102 double [] splineWs
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	103 ) {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	104 this.spline = spline;
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	105 this.splineQs = splineQs;
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	106 this.splineWs = splineWs;
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	107 }
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	108
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	109 public double [][] sample(
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	110 int numSamples,
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	111 double km,
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	112 Calculation errors
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	113 ) {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	114 double minQ = getQMin();
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	115 double maxQ = getQMax();
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	116
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	117 double [] outWs = new double[numSamples];
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	118 double [] outQs = new double[numSamples];
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	119
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	120 Arrays.fill(outWs, Double.NaN);
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	121 Arrays.fill(outQs, Double.NaN);
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	122
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	123 double stepWidth = (maxQ - minQ)/numSamples;
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	124
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	125 try {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	126 double q = minQ;
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	127 for (int i = 0; i < outWs.length; ++i, q += stepWidth) {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	128 outWs[i] = spline.value(outQs[i] = q);
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	129 }
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	130 }
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	131 catch (ArgumentOutsideDomainException aode) {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	132 if (errors != null) {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	133 // TODO: I18N
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	134 errors.addProblem(km, "spline interpolation failed");
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	135 }
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	136 log.error("spline interpolation failed.", aode);
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	137 }
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	138
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	139 return new double [][] { outWs, outQs };
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	140 }
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	141
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	142 public double getQMin() {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	143 return Math.min(splineQs[0], splineQs[splineQs.length-1]);
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	144 }
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	145
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	146 public double getQMax() {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	147 return Math.max(splineQs[0], splineQs[splineQs.length-1]);
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	148 }
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	149
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	150 /** Constructs a continues index between the columns to Qs. */
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	151 public PolynomialSplineFunction createIndexQRelation() {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	152
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	153 double [] indices = new double[splineQs.length];
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	154 for (int i = 0; i < indices.length; ++i) {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	155 indices[i] = i;
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	156 }
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	157
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	158 try {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	159 SplineInterpolator interpolator = new SplineInterpolator();
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	160 return interpolator.interpolate(indices, splineQs);
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	161 }
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	162 catch (MathIllegalArgumentException miae) {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	163 // Ignore me!
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	164 }
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	165 return null;
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	166 }
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	167 } // class SplineFunction
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	168
1721 eb35570df0e8 Cosmetics, docs. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1681 diff changeset	169 /**
1838 853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	170 * A row, typically a position where measurements were taken.
1721 eb35570df0e8 Cosmetics, docs. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1681 diff changeset	171 */
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	172 public static final class Row
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	173 implements Serializable, Comparable<Row>
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	174 {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	175 double km;
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	176 double [] ws;
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	177
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	178 public Row() {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	179 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	180
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	181 public Row(double km) {
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	182 this.km = km;
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	183 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	184
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	185 public Row(double km, double [] ws) {
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	186 this(km);
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	187 this.ws = ws;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	188 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	189
1838 853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	190 /**
853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	191 * Compare according to place of measurement (km).
853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	192 */
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	193 public int compareTo(Row other) {
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	194 double d = km - other.km;
458 523a256451cd Use unsharp km lookup (epsilon = 0.0001). Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 443 diff changeset	195 if (d < -0.0001) return -1;
523a256451cd Use unsharp km lookup (epsilon = 0.0001). Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 443 diff changeset	196 if (d > 0.0001) return +1;
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	197 return 0;
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	198 }
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	199
1887 d2e217f438bc Minor doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1860 diff changeset	200 /**
d2e217f438bc Minor doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1860 diff changeset	201 * Interpolate Ws, given Qs and a km.
d2e217f438bc Minor doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1860 diff changeset	202 *
d2e217f438bc Minor doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1860 diff changeset	203 * @param iqs Given ("input") Qs.
d2e217f438bc Minor doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1860 diff changeset	204 * @param ows Resulting ("output") Ws.
d2e217f438bc Minor doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1860 diff changeset	205 * @param table Table of which to use data for interpolation.
d2e217f438bc Minor doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1860 diff changeset	206 */
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	207 public void interpolateW(
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	208 Row other,
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	209 double km,
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	210 double [] iqs,
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	211 double [] ows,
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	212 WstValueTable table,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	213 Calculation errors
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	214 ) {
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	215 double kmWeight = Linear.factor(km, this.km, other.km);
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	216
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	217 QPosition qPosition = new QPosition();
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	218
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	219 for (int i = 0; i < iqs.length; ++i) {
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	220 if (table.getQPosition(km, iqs[i], qPosition) != null) {
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	221 double wt = getW(qPosition);
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	222 double wo = other.getW(qPosition);
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	223 if (Double.isNaN(wt) \|\| Double.isNaN(wo)) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	224 if (errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	225 // TODO: I18N
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	226 errors.addProblem(
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	227 km, "cannot find w for q = " + iqs[i]);
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	228 }
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	229 ows[i] = Double.NaN;
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	230 }
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	231 else {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	232 ows[i] = Linear.weight(kmWeight, wt, wo);
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	233 }
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	234 }
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	235 else {
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	236 if (errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	237 // TODO: I18N
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	238 errors.addProblem(km, "cannot find q = " + iqs[i]);
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	239 }
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	240 ows[i] = Double.NaN;
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	241 }
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	242 }
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	243 }
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	244
1938 1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	245
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	246 public SplineFunction createSpline(
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	247 WstValueTable table,
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	248 Calculation errors
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	249 ) {
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	250 int W = ws.length;
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	251
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	252 if (W < 1) {
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	253 if (errors != null) {
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	254 // TODO: I18N
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	255 errors.addProblem(km, "no ws found");
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	256 }
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	257 return null;
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	258 }
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	259
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	260 double [] splineQs = new double[W];
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	261
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	262 for (int i = 0; i < W; ++i) {
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	263 double sq = table.getQIndex(i, km);
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	264 if (Double.isNaN(sq) && errors != null) {
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	265 // TODO: I18N
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	266 errors.addProblem(
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	267 km, "no q found in " + (i+1) + " column");
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	268 }
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	269 splineQs[i] = sq;
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	270 }
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	271
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	272 try {
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	273 SplineInterpolator interpolator = new SplineInterpolator();
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	274 PolynomialSplineFunction spline =
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	275 interpolator.interpolate(splineQs, ws);
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	276
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	277 return new SplineFunction(spline, splineQs, ws);
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	278 }
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	279 catch (MathIllegalArgumentException miae) {
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	280 if (errors != null) {
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	281 // TODO: I18N
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	282 errors.addProblem(km, "spline creation failed");
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	283 }
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	284 log.error("spline creation failed", miae);
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	285 }
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	286 return null;
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	287 }
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	288
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	289 public SplineFunction createSpline(
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	290 Row other,
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	291 double km,
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	292 WstValueTable table,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	293 Calculation errors
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	294 ) {
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	295 int W = Math.min(ws.length, other.ws.length);
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	296
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	297 if (W < 1) {
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	298 if (errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	299 // TODO: I18N
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	300 errors.addProblem("no ws found");
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	301 }
1938 1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	302 return null;
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	303 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	304
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	305 double factor = Linear.factor(km, this.km, other.km);
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	306
1938 1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	307 double [] splineQs = new double[W];
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	308 double [] splineWs = new double[W];
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	309
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	310 for (int i = 0; i < W; ++i) {
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	311 double wws = Linear.weight(factor, ws[i], other.ws[i]);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	312 double wqs = Linear.weight(
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	313 factor,
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	314 table.getQIndex(i, km),
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	315 table.getQIndex(i, other.km));
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	316
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	317 if (Double.isNaN(wws) \|\| Double.isNaN(wqs)) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	318 if (errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	319 // TODO: I18N
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	320 errors.addProblem(km, "cannot find w or q");
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	321 }
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	322 }
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	323
1938 1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	324 splineWs[i] = wws;
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	325 splineQs[i] = wqs;
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	326 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	327
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	328 SplineInterpolator interpolator = new SplineInterpolator();
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	329
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	330 try {
1938 1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	331 PolynomialSplineFunction spline =
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	332 interpolator.interpolate(splineQs, splineWs);
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	333
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	334 return new SplineFunction(spline, splineQs, splineWs);
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	335 }
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	336 catch (MathIllegalArgumentException miae) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	337 if (errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	338 // TODO: I18N
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	339 errors.addProblem(km, "spline creation failed");
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	340 }
1860 ddd425858169 Improved logging to include exception in fail case of spline interpolation of wst. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1838 diff changeset	341 log.error("spline creation failed", miae);
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	342 }
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	343
1938 1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	344 return null;
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	345 }
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	346
1938 1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	347 public double [][] interpolateWQ(
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	348 Row other,
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	349 double km,
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	350 int steps,
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	351 WstValueTable table,
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	352 Calculation errors
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	353 ) {
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	354 SplineFunction sf = createSpline(other, km, table, errors);
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	355
1938 1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	356 return sf != null
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	357 ? sf.sample(steps, km, errors)
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	358 : new double[2][0];
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	359 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	360
1721 eb35570df0e8 Cosmetics, docs. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1681 diff changeset	361
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	362 public double [][] interpolateWQ(
742 c09c9e05ecfa Removed trailing whitespace. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 704 diff changeset	363 int steps,
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	364 WstValueTable table,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	365 Calculation errors
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	366 ) {
1938 1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	367 SplineFunction sf = createSpline(table, errors);
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	368
1938 1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	369 return sf != null
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	370 ? sf.sample(steps, km, errors)
1d991c91285b Refactored the code of the 'berechnete Abflusskurve' to be reusable in the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1937 diff changeset	371 : new double[2][0];
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	372 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	373
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	374
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	375 public double getW(QPosition qPosition) {
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	376 int index = qPosition.index;
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	377 double weight = qPosition.weight;
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	378
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	379 return weight == 1.0
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	380 ? ws[index]
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	381 : Linear.weight(weight, ws[index-1], ws[index]);
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	382 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	383
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	384 public double getW(
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	385 Row other,
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	386 double km,
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	387 QPosition qPosition
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	388 ) {
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	389 double kmWeight = Linear.factor(km, this.km, other.km);
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	390
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	391 int index = qPosition.index;
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	392 double weight = qPosition.weight;
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	393
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	394 double tw, ow;
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	395
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	396 if (weight == 1.0) {
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	397 tw = ws[index];
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	398 ow = other.ws[index];
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	399 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	400 else {
742 c09c9e05ecfa Removed trailing whitespace. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 704 diff changeset	401 tw = Linear.weight(weight, ws[index-1], ws[index]);
c09c9e05ecfa Removed trailing whitespace. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 704 diff changeset	402 ow = Linear.weight(weight, other.ws[index-1], other.ws[index]);
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	403 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	404
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	405 return Linear.weight(kmWeight, tw, ow);
427 909196be11a0 First step to calculate "W fuer ungleichwertige Abfluesse" correctly. flys/issue55 Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 426 diff changeset	406 }
1937 f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	407
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	408 public double [] findQsForW(double w, WstValueTable table) {
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	409
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	410 TDoubleArrayList qs = new TDoubleArrayList();
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	411
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	412 if (ws.length > 0 && Math.abs(ws[0]-w) < 0.000001) {
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	413 double q = table.getQIndex(0, km);
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	414 if (!Double.isNaN(q)) {
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	415 qs.add(q);
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	416 }
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	417 }
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	418
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	419 for (int i = 1; i < ws.length; ++i) {
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	420 double w2 = ws[i];
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	421 if (Double.isNaN(w2)) {
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	422 continue;
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	423 }
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	424 if (Math.abs(w2-w) < 0.000001) {
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	425 double q = table.getQIndex(i, km);
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	426 if (!Double.isNaN(q)) {
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	427 qs.add(q);
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	428 }
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	429 continue;
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	430 }
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	431 double w1 = ws[i-1];
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	432 if (Double.isNaN(w1)) {
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	433 continue;
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	434 }
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	435
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	436 if (w < Math.min(w1, w2) \|\| w > Math.max(w1, w2)) {
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	437 continue;
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	438 }
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	439
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	440 double q1 = table.getQIndex(i-1, km);
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	441 double q2 = table.getQIndex(i, km);
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	442 if (Double.isNaN(q1) \|\| Double.isNaN(q2)) {
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	443 continue;
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	444 }
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	445
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	446 double q = Linear.linear(w, w1, w2, q1, q2);
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	447 qs.add(q);
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	448 }
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	449
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	450 return qs.toNativeArray();
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	451 }
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	452
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	453 public double [] findQsForW(
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	454 Row other,
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	455 double w,
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	456 double km,
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	457 WstValueTable table
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	458 ) {
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	459 TDoubleArrayList qs = new TDoubleArrayList();
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	460
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	461 double factor = Linear.factor(km, this.km, other.km);
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	462
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	463 if (ws.length > 0) {
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	464 double wt = Linear.weight(factor, ws[0], other.ws[0]);
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	465 if (!Double.isNaN(wt)) {
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	466 double q = table.getQIndex(0, km);
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	467 if (!Double.isNaN(q)) {
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	468 qs.add(q);
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	469 }
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	470 }
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	471 }
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	472
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	473 for (int i = 1; i < ws.length; ++i) {
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	474 double w2 = Linear.weight(factor, ws[i], other.ws[i]);
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	475 if (Double.isNaN(w2)) {
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	476 continue;
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	477 }
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	478 if (Math.abs(w2-w) < 0.000001) {
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	479 double q = table.getQIndex(i, km);
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	480 if (!Double.isNaN(q)) {
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	481 qs.add(q);
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	482 }
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	483 continue;
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	484 }
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	485 double w1 = Linear.weight(factor, ws[i-1], other.ws[i-1]);
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	486 if (Double.isNaN(w1)) {
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	487 continue;
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	488 }
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	489
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	490 if (w < Math.min(w1, w2) \|\| w > Math.max(w1, w2)) {
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	491 continue;
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	492 }
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	493
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	494 double q1 = table.getQIndex(i-1, km);
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	495 double q2 = table.getQIndex(i, km);
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	496 if (Double.isNaN(q1) \|\| Double.isNaN(q2)) {
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	497 continue;
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	498 }
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	499
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	500 double q = Linear.linear(w, w1, w2, q1, q2);
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	501 qs.add(q);
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	502 }
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	503
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	504 return qs.toNativeArray();
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	505 }
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	506 } // class Row
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	507
1838 853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	508 /** Rows in table. */
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	509 protected List<Row> rows;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	510
1838 853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	511 /** Columns in table. */
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	512 protected Column [] columns;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	513
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	514 public WstValueTable() {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	515 rows = new ArrayList<Row>();
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	516 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	517
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	518 public WstValueTable(Column [] columns) {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	519 this();
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	520 this.columns = columns;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	521 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	522
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	523 public WstValueTable(Column [] columns, List<Row> rows) {
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	524 this.columns = columns;
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	525 this.rows = rows;
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	526 }
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	527
1838 853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	528 /** Sort rows (by km). */
458 523a256451cd Use unsharp km lookup (epsilon = 0.0001). Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 443 diff changeset	529 public void sortRows() {
523a256451cd Use unsharp km lookup (epsilon = 0.0001). Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 443 diff changeset	530 Collections.sort(rows);
523a256451cd Use unsharp km lookup (epsilon = 0.0001). Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 443 diff changeset	531 }
523a256451cd Use unsharp km lookup (epsilon = 0.0001). Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 443 diff changeset	532
1887 d2e217f438bc Minor doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1860 diff changeset	533 /**
d2e217f438bc Minor doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1860 diff changeset	534 * @param km Given kilometer.
d2e217f438bc Minor doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1860 diff changeset	535 * @param qs Given Q values.
d2e217f438bc Minor doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1860 diff changeset	536 * @param ws output parameter.
d2e217f438bc Minor doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1860 diff changeset	537 */
380 5f55047a17e8 Moved expensive allocation of result arrays out of km loop. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 374 diff changeset	538 public double [] interpolateW(double km, double [] qs, double [] ws) {
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	539 return interpolateW(km, qs, ws, null);
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	540 }
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	541
1098 1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	542
1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	543 /**
1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	544 * @param ws (output parameter), gets returned.
1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	545 * @return output parameter ws.
1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	546 */
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	547 public double [] interpolateW(
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	548 double km,
742 c09c9e05ecfa Removed trailing whitespace. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 704 diff changeset	549 double [] qs,
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	550 double [] ws,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	551 Calculation errors
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	552 ) {
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	553 int rowIndex = Collections.binarySearch(rows, new Row(km));
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	554
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	555 QPosition qPosition = new QPosition();
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	556
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	557 if (rowIndex >= 0) { // direct row match
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	558 Row row = rows.get(rowIndex);
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	559 for (int i = 0; i < qs.length; ++i) {
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	560 if (getQPosition(km, qs[i], qPosition) == null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	561 if (errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	562 // TODO: I18N
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	563 errors.addProblem(km, "cannot find q = " + qs[i]);
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	564 }
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	565 ws[i] = Double.NaN;
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	566 }
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	567 else {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	568 if (Double.isNaN(ws[i] = row.getW(qPosition))
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	569 && errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	570 // TODO: I18N
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	571 errors.addProblem(
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	572 km, "cannot find w for q = " + qs[i]);
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	573 }
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	574 }
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	575 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	576 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	577 else { // needs bilinear interpolation
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	578 rowIndex = -rowIndex -1;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	579
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	580 if (rowIndex < 1 \|\| rowIndex >= rows.size()) {
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	581 // do not extrapolate
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	582 Arrays.fill(ws, Double.NaN);
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	583 if (errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	584 // TODO: I18N
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	585 errors.addProblem(km, "km not found");
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	586 }
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	587 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	588 else {
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	589 Row r1 = rows.get(rowIndex-1);
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	590 Row r2 = rows.get(rowIndex);
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	591 r1.interpolateW(r2, km, qs, ws, this, errors);
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	592 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	593 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	594
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	595 return ws;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	596 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	597
1838 853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	598 /**
853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	599 * Interpolate W and Q values at a given km.
853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	600 */
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	601 public double [][] interpolateWQ(double km) {
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	602 return interpolateWQ(km, null);
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	603 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	604
1838 853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	605 /**
853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	606 * Interpolate W and Q values at a given km.
853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	607 */
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	608 public double [][] interpolateWQ(double km, Calculation errors) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	609 return interpolateWQ(km, DEFAULT_Q_STEPS, errors);
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	610 }
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	611
1860 ddd425858169 Improved logging to include exception in fail case of spline interpolation of wst. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1838 diff changeset	612 /**
ddd425858169 Improved logging to include exception in fail case of spline interpolation of wst. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1838 diff changeset	613 * Interpolate W and Q values at a given km.
ddd425858169 Improved logging to include exception in fail case of spline interpolation of wst. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1838 diff changeset	614 */
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	615 public double [][] interpolateWQ(double km, int steps, Calculation errors) {
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	616
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	617 int rowIndex = Collections.binarySearch(rows, new Row(km));
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	618
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	619 if (rowIndex >= 0) { // direct row match
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	620 Row row = rows.get(rowIndex);
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	621 return row.interpolateWQ(steps, this, errors);
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	622 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	623
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	624 rowIndex = -rowIndex -1;
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	625
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	626 if (rowIndex < 1 \|\| rowIndex >= rows.size()) {
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	627 // do not extrapolate
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	628 if (errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	629 // TODO: I18N
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	630 errors.addProblem(km, "km not found");
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	631 }
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	632 return new double[2][0];
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	633 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	634
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	635 Row r1 = rows.get(rowIndex-1);
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	636 Row r2 = rows.get(rowIndex);
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	637
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	638 return r1.interpolateWQ(r2, km, steps, this, errors);
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	639 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	640
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	641 public boolean interpolate(
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	642 double km,
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	643 double [] out,
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	644 QPosition qPosition,
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	645 Function qFunction
649 44175d4720f8 Fix crashes in back jump correction. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 645 diff changeset	646 ) {
427 909196be11a0 First step to calculate "W fuer ungleichwertige Abfluesse" correctly. flys/issue55 Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 426 diff changeset	647 int R1 = rows.size()-1;
909196be11a0 First step to calculate "W fuer ungleichwertige Abfluesse" correctly. flys/issue55 Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 426 diff changeset	648
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	649 out[1] = qFunction.value(getQ(qPosition, km));
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	650
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	651 if (Double.isNaN(out[1])) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	652 return false;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	653 }
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	654
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	655 QPosition nPosition = new QPosition();
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	656 if (getQPosition(km, out[1], nPosition) == null) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	657 return false;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	658 }
427 909196be11a0 First step to calculate "W fuer ungleichwertige Abfluesse" correctly. flys/issue55 Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 426 diff changeset	659
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	660 int rowIndex = Collections.binarySearch(rows, new Row(km));
427 909196be11a0 First step to calculate "W fuer ungleichwertige Abfluesse" correctly. flys/issue55 Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 426 diff changeset	661
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	662 if (rowIndex >= 0) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	663 // direct row match
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	664 out[0] = rows.get(rowIndex).getW(nPosition);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	665 return !Double.isNaN(out[0]);
427 909196be11a0 First step to calculate "W fuer ungleichwertige Abfluesse" correctly. flys/issue55 Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 426 diff changeset	666 }
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	667
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	668 rowIndex = -rowIndex -1;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	669
1681 e5f7f25a511c fixed flys/issue201 Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1098 diff changeset	670 if (rowIndex < 1 \|\| rowIndex > R1) {
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	671 // do not extrapolate
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	672 return false;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	673 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	674
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	675 Row r1 = rows.get(rowIndex-1);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	676 Row r2 = rows.get(rowIndex);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	677 out[0] = r1.getW(r2, km, nPosition);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	678
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	679 return !Double.isNaN(out[0]);
427 909196be11a0 First step to calculate "W fuer ungleichwertige Abfluesse" correctly. flys/issue55 Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 426 diff changeset	680 }
909196be11a0 First step to calculate "W fuer ungleichwertige Abfluesse" correctly. flys/issue55 Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 426 diff changeset	681
1098 1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	682
1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	683 /**
1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	684 * Look up interpolation of a Q at given positions.
1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	685 *
1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	686 * @param q the non-interpolated Q value.
1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	687 * @param referenceKm the reference km (e.g. gauge position).
1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	688 * @param kms positions for which to interpolate.
1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	689 * @param ws (output) resulting interpolated ws.
1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	690 * @param qs (output) resulting interpolated qs.
1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	691 * @param errors calculation object to store errors.
1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	692 */
427 909196be11a0 First step to calculate "W fuer ungleichwertige Abfluesse" correctly. flys/issue55 Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 426 diff changeset	693 public QPosition interpolate(
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	694 double q,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	695 double referenceKm,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	696 double [] kms,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	697 double [] ws,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	698 double [] qs,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	699 Calculation errors
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	700 ) {
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	701 return interpolate(
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	702 q, referenceKm, kms, ws, qs, 0, kms.length, errors);
649 44175d4720f8 Fix crashes in back jump correction. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 645 diff changeset	703 }
44175d4720f8 Fix crashes in back jump correction. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 645 diff changeset	704
44175d4720f8 Fix crashes in back jump correction. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 645 diff changeset	705 public QPosition interpolate(
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	706 double q,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	707 double referenceKm,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	708 double [] kms,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	709 double [] ws,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	710 double [] qs,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	711 int startIndex,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	712 int length,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	713 Calculation errors
649 44175d4720f8 Fix crashes in back jump correction. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 645 diff changeset	714 ) {
645 433f67a076aa Use an explicit reference km for interpolation now. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 633 diff changeset	715 QPosition qPosition = getQPosition(referenceKm, q);
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	716
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	717 if (qPosition == null) {
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	718 // we cannot locate q at km
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	719 Arrays.fill(ws, Double.NaN);
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	720 Arrays.fill(qs, Double.NaN);
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	721 if (errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	722 errors.addProblem(referenceKm, "cannot find q");
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	723 }
427 909196be11a0 First step to calculate "W fuer ungleichwertige Abfluesse" correctly. flys/issue55 Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 426 diff changeset	724 return null;
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	725 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	726
645 433f67a076aa Use an explicit reference km for interpolation now. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 633 diff changeset	727 Row kmKey = new Row();
433f67a076aa Use an explicit reference km for interpolation now. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 633 diff changeset	728
433f67a076aa Use an explicit reference km for interpolation now. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 633 diff changeset	729 int R1 = rows.size()-1;
433f67a076aa Use an explicit reference km for interpolation now. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 633 diff changeset	730
649 44175d4720f8 Fix crashes in back jump correction. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 645 diff changeset	731 for (int i = startIndex, end = startIndex+length; i < end; ++i) {
458 523a256451cd Use unsharp km lookup (epsilon = 0.0001). Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 443 diff changeset	732
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	733 if (Double.isNaN(qs[i] = getQ(qPosition, kms[i]))) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	734 if (errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	735 errors.addProblem(kms[i], "cannot find q");
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	736 }
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	737 ws[i] = Double.NaN;
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	738 continue;
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	739 }
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	740
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	741 kmKey.km = kms[i];
645 433f67a076aa Use an explicit reference km for interpolation now. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 633 diff changeset	742 int rowIndex = Collections.binarySearch(rows, kmKey);
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	743
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	744 if (rowIndex >= 0) {
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	745 // direct row match
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	746 if (Double.isNaN(ws[i] = rows.get(rowIndex).getW(qPosition))
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	747 && errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	748 errors.addProblem(kms[i], "cannot find w");
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	749 }
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	750 continue;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	751 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	752
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	753 rowIndex = -rowIndex -1;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	754
1681 e5f7f25a511c fixed flys/issue201 Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1098 diff changeset	755 if (rowIndex < 1 \|\| rowIndex > R1) {
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	756 // do not extrapolate
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	757 if (errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	758 errors.addProblem(kms[i], "cannot find km");
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	759 }
1888 abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	760 ws[i] = Double.NaN;
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	761 continue;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	762 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	763 Row r1 = rows.get(rowIndex-1);
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	764 Row r2 = rows.get(rowIndex);
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	765
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	766 if (Double.isNaN(ws[i] = r1.getW(r2, kms[i], qPosition))
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	767 && errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	768 errors.addProblem(kms[i], "cannot find w");
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	769 }
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	770 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	771
427 909196be11a0 First step to calculate "W fuer ungleichwertige Abfluesse" correctly. flys/issue55 Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 426 diff changeset	772 return qPosition;
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	773 }
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	774
1888 abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	775 /**
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	776 * Linearly interpolate w at a km at a column of two rows.
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	777 *
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	778 * @param km position for which to interpolate.
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	779 * @param row1 first row.
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	780 * @param row2 second row.
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	781 * @param col column-index at which to look.
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	782 *
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	783 * @return Linearly interpolated w, NaN if one of the given rows was null.
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	784 */
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	785 public static double linearW(double km, Row row1, Row row2, int col) {
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	786 if (row1 == null \|\| row2 == null) {
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	787 return Double.NaN;
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	788 }
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	789
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	790 return Linear.linear(km,
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	791 row1.km, row2.km,
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	792 row1.ws[col], row2.ws[col]);
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	793 }
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	794
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	795 /**
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	796 * Do interpolation/lookup of W and Q within columns (i.e. ignoring values
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	797 * of other columns).
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	798 * @param km position (km) at which to interpolate/lookup.
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	799 * @return [[q0, q1, .. qx] , [w0, w1, .. wx]] (can contain NaNs)
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	800 */
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	801 public double [][] interpolateWQColumnwise(double km) {
1919 9bec7d2f8c35 Prevent ArrayIndexOutOfBounds, log method entry. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1907 diff changeset	802 log.debug("WstValueTable.interpolateWQColumnwise");
1888 abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	803 double [] qs = new double[columns.length];
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	804 double [] ws = new double[columns.length];
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	805
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	806 // Find out row from where we will start searching.
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	807 int rowIndex = Collections.binarySearch(rows, new Row(km));
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	808
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	809 if (rowIndex < 0) {
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	810 rowIndex = -rowIndex -1;
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	811 }
1937 f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	812 if (rowIndex >= rows.size()) {
1919 9bec7d2f8c35 Prevent ArrayIndexOutOfBounds, log method entry. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1907 diff changeset	813 rowIndex = rows.size() -1;
1937 f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	814 }
1888 abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	815
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	816 Row startRow = rows.get(rowIndex);
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	817
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	818 for (int col = 0; col < columns.length; col++) {
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	819 qs[col] = columns[col].getQRangeTree().findQ(km);
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	820 if (startRow.km == km && startRow.ws[col] != Double.NaN) {
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	821 // Great. W is defined at km.
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	822 ws[col] = startRow.ws[col];
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	823 continue;
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	824 }
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	825
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	826 // Search neighbouring rows that define w at this col.
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	827 Row rowBefore = null;
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	828 Row rowAfter = null;
1899 d002198c64a2 Removed logging noise, find better suited rows for W interpolation. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1888 diff changeset	829 for (int before = rowIndex -1; before >= 0; before--) {
1888 abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	830 if (!Double.isNaN(rows.get(before).ws[col])) {
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	831 rowBefore = rows.get(before);
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	832 break;
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	833 }
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	834 }
1907 6c7fdfd262ac Do not re-evaluate constant size() in for-loops. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1899 diff changeset	835 for (int after = rowIndex, R = rows.size(); after < R; after++) {
1888 abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	836 if (!Double.isNaN(rows.get(after).ws[col])) {
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	837 rowAfter = rows.get(after);
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	838 break;
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	839 }
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	840 }
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	841
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	842 ws[col] = linearW(km, rowBefore, rowAfter, col);
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	843 }
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	844
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	845 return new double [][] {qs, ws};
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	846 }
abf5da3b285c Added new (at fixed columns) interpolation mechanism. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1887 diff changeset	847
1937 f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	848 public double [] findQsForW(double km, double w) {
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	849
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	850 int rowIndex = Collections.binarySearch(rows, new Row(km));
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	851
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	852 if (rowIndex >= 0) {
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	853 return rows.get(rowIndex).findQsForW(w, this);
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	854 }
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	855
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	856 rowIndex = -rowIndex - 1;
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	857
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	858 if (rowIndex < 1 \|\| rowIndex >= rows.size()) {
1939 2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	859 // Do not extrapolate.
1937 f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	860 return new double[0];
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	861 }
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	862
1939 2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	863 // Needs bilinear interpolation.
1937 f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	864 Row r1 = rows.get(rowIndex-1);
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	865 Row r2 = rows.get(rowIndex);
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	866
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	867 return r1.findQsForW(r2, w, km, this);
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	868 }
f07d64d5cbe1 'W auf freier Strecke' calculation. Fetch corresponding Qs for given Ws from the WST model Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1919 diff changeset	869
1939 2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	870 protected SplineFunction createSpline(double km, Calculation errors) {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	871
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	872 int rowIndex = Collections.binarySearch(rows, new Row(km));
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	873
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	874 if (rowIndex >= 0) {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	875 SplineFunction sf = rows.get(rowIndex).createSpline(this, errors);
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	876 if (sf == null && errors != null) {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	877 // TODO: I18N
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	878 errors.addProblem(km, "cannot create w/q relation");
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	879 }
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	880 return sf;
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	881 }
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	882
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	883 rowIndex = -rowIndex - 1;
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	884
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	885 if (rowIndex < 1 \|\| rowIndex >= rows.size()) {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	886 // Do not extrapolate.
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	887 if (errors != null) {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	888 // TODO: I18N
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	889 errors.addProblem(km, "km not found");
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	890 }
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	891 return null;
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	892 }
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	893
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	894 // Needs bilinear interpolation.
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	895 Row r1 = rows.get(rowIndex-1);
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	896 Row r2 = rows.get(rowIndex);
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	897
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	898 SplineFunction sf = r1.createSpline(r2, km, this, errors);
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	899 if (sf == null && errors != null) {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	900 // TODO: I18N
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	901 errors.addProblem(km, "cannot create w/q relation");
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	902 }
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	903
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	904 return sf;
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	905 }
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	906
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	907 /** 'Bezugslinienverfahren' */
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	908 public double [][] relateWs(
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	909 double km1,
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	910 double km2,
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	911 int numSamples,
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	912 Calculation errors
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	913 ) {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	914 SplineFunction sf1 = createSpline(km1, errors);
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	915 if (sf1 == null) {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	916 return new double[2][0];
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	917 }
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	918
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	919 SplineFunction sf2 = createSpline(km1, errors);
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	920 if (sf2 == null) {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	921 return new double[2][0];
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	922 }
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	923
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	924 PolynomialSplineFunction iQ1 = sf1.createIndexQRelation();
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	925 if (iQ1 == null) {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	926 if (errors != null) {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	927 // TODO: I18N
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	928 errors.addProblem(km1, "cannot create index/q relation");
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	929 }
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	930 return new double[2][0];
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	931 }
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	932
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	933 PolynomialSplineFunction iQ2 = sf1.createIndexQRelation();
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	934 if (iQ2 == null) {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	935 if (errors != null) {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	936 // TODO: I18N
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	937 errors.addProblem(km2, "cannot create index/q relation");
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	938 }
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	939 return new double[2][0];
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	940 }
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	941
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	942 int N = sf1.splineQs.length;
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	943 double stepWidth = N/(double)numSamples;
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	944
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	945 PolynomialSplineFunction qW1 = sf1.spline;
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	946 PolynomialSplineFunction qW2 = sf2.spline;
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	947
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	948 double [] ws1 = new double[numSamples];
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	949 double [] ws2 = new double[numSamples];
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	950
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	951 Arrays.fill(ws1, Double.NaN);
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	952 Arrays.fill(ws2, Double.NaN);
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	953
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	954 boolean hadErrors = false;
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	955
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	956 double p = 0d;
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	957 for (int i = 0; i < numSamples; ++i, p += stepWidth) {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	958 try {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	959 double q1 = iQ1.value(p);
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	960 double w1 = qW1.value(q1);
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	961 double q2 = iQ2.value(p);
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	962 double w2 = qW2.value(q2);
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	963 ws1[i] = w1;
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	964 ws2[i] = w2;
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	965 }
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	966 catch (ArgumentOutsideDomainException aode) {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	967 if (!hadErrors) {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	968 // XXX: I'm not sure if this really can happen
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	969 // and if we should report this more than once.
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	970 hadErrors = true;
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	971 if (errors != null) {
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	972 // TODO: I18N
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	973 log.debug("W~W failed", aode);
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	974 errors.addProblem("W~W failed");
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	975 }
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	976 }
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	977 }
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	978 }
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	979
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	980 return new double [][] { ws1, ws2 };
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	981 }
2730d17df021 Added the implementation of the 'Bezugslinienverfahren'. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1938 diff changeset	982
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	983 public QPosition getQPosition(double km, double q) {
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	984 return getQPosition(km, q, new QPosition());
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	985 }
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	986
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	987 public QPosition getQPosition(double km, double q, QPosition qPosition) {
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	988
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	989 if (columns.length == 0) {
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	990 return null;
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	991 }
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	992
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	993 double qLast = columns[0].getQRangeTree().findQ(km);
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	994
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	995 if (Math.abs(qLast - q) < 0.00001) {
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	996 return qPosition.set(0, 1d);
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	997 }
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	998
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	999 for (int i = 1; i < columns.length; ++i) {
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1000 double qCurrent = columns[i].getQRangeTree().findQ(km);
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1001 if (Math.abs(qCurrent - q) < 0.00001) {
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1002 return qPosition.set(i, 1d);
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1003 }
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1004
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1005 double qMin, qMax;
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1006 if (qLast < qCurrent) { qMin = qLast; qMax = qCurrent; }
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1007 else { qMin = qCurrent; qMax = qLast; }
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1008
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1009 if (q > qMin && q < qMax) {
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	1010 double weight = Linear.factor(q, qLast, qCurrent);
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1011 return qPosition.set(i, weight);
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1012 }
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1013 qLast = qCurrent;
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1014 }
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1015
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1016 return null;
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1017 }
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1018
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1019 public double getQIndex(int index, double km) {
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1020 return columns[index].getQRangeTree().findQ(km);
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1021 }
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1022
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1023 public double getQ(QPosition qPosition, double km) {
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1024 int index = qPosition.index;
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1025 double weight = qPosition.weight;
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1026
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1027 if (weight == 1d) {
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1028 return columns[index].getQRangeTree().findQ(km);
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1029 }
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1030 double q1 = columns[index-1].getQRangeTree().findQ(km);
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1031 double q2 = columns[index ].getQRangeTree().findQ(km);
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	1032 return Linear.weight(weight, q1, q2);
633 d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	1033 }
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	1034 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	1035 // 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/WstValueTable.java @ 2072:4cdd9c4896f6