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 @ 1860:ddd425858169

Improved logging to include exception in fail case of spline interpolation of wst. flys-artifacts/trunk@3209 c6561f87-3c4e-4783-a992-168aeb5c3f6f

author	Felix Wolfsteller <felix.wolfsteller@intevation.de>
date	Thu, 10 Nov 2011 10:13:57 +0000
parents	853cd2120d69
children	d2e217f438bc

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
1838 853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	23 /**
853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	24 * 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	25 */
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	26 public class WstValueTable
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	27 implements Serializable
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	28 {
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	29 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	30
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	31 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	32
1721 eb35570df0e8 Cosmetics, docs. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1681 diff changeset	33 /**
eb35570df0e8 Cosmetics, docs. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1681 diff changeset	34 * A Column in the table, typically representing one measurement session.
eb35570df0e8 Cosmetics, docs. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1681 diff changeset	35 */
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	36 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	37 implements Serializable
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	38 {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	39 protected String name;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	40
632 07640ab913fd First part of storing qs in ranges Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 621 diff changeset	41 protected QRangeTree qRangeTree;
07640ab913fd First part of storing qs in ranges Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 621 diff changeset	42
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	43 public Column() {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	44 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	45
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	46 public Column(String name) {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	47 this.name = name;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	48 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	49
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	50 public String getName() {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	51 return name;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	52 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	53
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	54 public void setName(String name) {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	55 this.name = name;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	56 }
632 07640ab913fd First part of storing qs in ranges Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 621 diff changeset	57
07640ab913fd First part of storing qs in ranges Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 621 diff changeset	58 public QRangeTree getQRangeTree() {
07640ab913fd First part of storing qs in ranges Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 621 diff changeset	59 return qRangeTree;
07640ab913fd First part of storing qs in ranges Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 621 diff changeset	60 }
07640ab913fd First part of storing qs in ranges Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 621 diff changeset	61
07640ab913fd First part of storing qs in ranges Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 621 diff changeset	62 public void setQRangeTree(QRangeTree qRangeTree) {
07640ab913fd First part of storing qs in ranges Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 621 diff changeset	63 this.qRangeTree = qRangeTree;
07640ab913fd First part of storing qs in ranges Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 621 diff changeset	64 }
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	65 } // class Column
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	66
1721 eb35570df0e8 Cosmetics, docs. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1681 diff changeset	67 /**
eb35570df0e8 Cosmetics, docs. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1681 diff changeset	68 * A (weighted) position used for interpolation.
eb35570df0e8 Cosmetics, docs. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1681 diff changeset	69 */
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	70 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	71
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	72 protected int index;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	73 protected double weight;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	74
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	75 public QPosition() {
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
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	78 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	79 this.index = index;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	80 this.weight = weight;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	81 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	82
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	83 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	84 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	85 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	86 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	87 }
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
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	89 } // class Position
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	90
1721 eb35570df0e8 Cosmetics, docs. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1681 diff changeset	91 /**
1838 853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	92 * A row, typically a position where measurements were taken.
1721 eb35570df0e8 Cosmetics, docs. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1681 diff changeset	93 */
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	94 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	95 implements Serializable, Comparable<Row>
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	96 {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	97 double km;
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	98 double [] ws;
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	99
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	100 public Row() {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	101 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	102
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	103 public Row(double km) {
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	104 this.km = km;
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	105 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	106
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	107 public Row(double km, double [] ws) {
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	108 this(km);
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	109 this.ws = ws;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	110 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	111
1838 853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	112 /**
853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	113 * Compare according to place of measurement (km).
853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	114 */
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	115 public int compareTo(Row other) {
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	116 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	117 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	118 if (d > 0.0001) return +1;
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	119 return 0;
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	120 }
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	121
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	122 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	123 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	124 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	125 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	126 double [] ows,
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	127 WstValueTable table,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	128 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	129 ) {
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	130 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	131
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	132 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	133
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	134 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	135 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	136 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	137 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	138 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	139 if (errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	140 // TODO: I18N
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	141 errors.addProblem(
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	142 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	143 }
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	144 ows[i] = Double.NaN;
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	145 }
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	146 else {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	147 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	148 }
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	149 }
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	150 else {
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	151 if (errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	152 // TODO: I18N
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	153 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	154 }
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	155 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	156 }
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	157 }
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	158 }
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	159
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	160 public double [][] interpolateWQ(
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	161 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	162 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	163 int steps,
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	164 WstValueTable table,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	165 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	166 ) {
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	167 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	168
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	169 if (W < 1) {
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	170 if (errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	171 // TODO: I18N
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	172 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	173 }
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	174 return new double[2][0];
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	175 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	176
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	177 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	178
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	179 double [] splineQ = new double[W];
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	180 double [] splineW = new double[W];
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	181
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	182 double minQ = Double.MAX_VALUE;
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	183 double maxQ = -Double.MAX_VALUE;
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	184
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	185 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	186 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	187 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	188 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	189 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	190 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	191
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	192 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	193 if (errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	194 // TODO: I18N
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	195 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	196 }
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	197 }
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	198 else {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	199 if (wqs < minQ) minQ = wqs;
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	200 if (wqs > maxQ) maxQ = wqs;
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	201 }
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	202
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	203 splineW[i] = wws;
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	204 splineQ[i] = wqs;
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	205 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	206
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	207 double stepWidth = (maxQ - minQ)/steps;
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	208
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	209 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	210 PolynomialSplineFunction spline;
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	211
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	212 try {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	213 spline = interpolator.interpolate(splineQ, splineW);
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	214 }
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	215 catch (MathIllegalArgumentException miae) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	216 if (errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	217 // TODO: I18N
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	218 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	219 }
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	220 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	221 return new double[2][0];
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	222 }
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	223
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	224 double [] outWs = new double[steps];
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	225 double [] outQs = new double[steps];
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	226
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	227 Arrays.fill(outWs, Double.NaN);
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	228 Arrays.fill(outQs, Double.NaN);
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	229
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	230 try {
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	231 double q = minQ;
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	232 for (int i = 0; i < outWs.length; ++i, q += stepWidth) {
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	233 outWs[i] = spline.value(outQs[i] = q);
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	234 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	235 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	236 catch (ArgumentOutsideDomainException aode) {
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	237 if (errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	238 // TODO: I18N
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	239 errors.addProblem(km, "spline interpolation failed");
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	240 }
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	241 log.error("spline interpolation failed", aode);
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	242 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	243
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	244 return new double [][] { outWs, outQs };
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	245 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	246
1721 eb35570df0e8 Cosmetics, docs. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1681 diff changeset	247
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	248 public double [][] interpolateWQ(
742 c09c9e05ecfa Removed trailing whitespace. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 704 diff changeset	249 int steps,
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	250 WstValueTable table,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	251 Calculation errors
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	252 ) {
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	253 int W = ws.length;
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	254
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	255 if (W < 1) {
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	256 if (errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	257 // TODO: I18N
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	258 errors.addProblem(km, "no ws found");
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	259 }
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	260 return new double[2][0];
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	261 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	262
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	263 double [] splineQ = new double[W];
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	264
742 c09c9e05ecfa Removed trailing whitespace. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 704 diff changeset	265 double minQ = Double.MAX_VALUE;
c09c9e05ecfa Removed trailing whitespace. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 704 diff changeset	266 double maxQ = -Double.MAX_VALUE;
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	267
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	268 for (int i = 0; i < W; ++i) {
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	269 double sq = table.getQIndex(i, km);
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	270 if (Double.isNaN(sq)) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	271 if (errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	272 // TODO: I18N
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	273 errors.addProblem(
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	274 km, "no q found in " + (i+1) + " column");
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	275 }
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	276 }
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	277 else {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	278 if (sq < minQ) minQ = sq;
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	279 if (sq > maxQ) maxQ = sq;
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	280 }
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	281 splineQ[i] = sq;
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	282 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	283
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	284 double stepWidth = (maxQ - minQ)/steps;
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	285
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	286 SplineInterpolator interpolator = new SplineInterpolator();
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	287
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	288 PolynomialSplineFunction spline;
742 c09c9e05ecfa Removed trailing whitespace. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 704 diff changeset	289
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	290 try {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	291 spline = interpolator.interpolate(splineQ, ws);
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	292 }
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	293 catch (MathIllegalArgumentException miae) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	294 if (errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	295 // TODO: I18N
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	296 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	297 }
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	298 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	299 return new double[2][0];
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	300 }
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	301
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	302 double [] outWs = new double[steps];
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	303 double [] outQs = new double[steps];
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	304
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	305 Arrays.fill(outWs, Double.NaN);
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	306 Arrays.fill(outQs, Double.NaN);
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	307
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	308 try {
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	309 double q = minQ;
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	310 for (int i = 0; i < outWs.length; ++i, q += stepWidth) {
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	311 outWs[i] = spline.value(outQs[i] = q);
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	312 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	313 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	314 catch (ArgumentOutsideDomainException aode) {
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	315 if (errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	316 // TODO: I18N
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	317 errors.addProblem(km, "spline interpolation failed");
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	318 }
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	319 log.error("spline interpolation failed.", aode);
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	320 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	321
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	322 return new double [][] { outWs, outQs };
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	323 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	324
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	325
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	326 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	327 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	328 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	329
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	330 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	331 ? ws[index]
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	332 : 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	333 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	334
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	335 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	336 Row other,
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	337 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	338 QPosition qPosition
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	339 ) {
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	340 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	341
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	342 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	343 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	344
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	345 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	346
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	347 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	348 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	349 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	350 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	351 else {
742 c09c9e05ecfa Removed trailing whitespace. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 704 diff changeset	352 tw = Linear.weight(weight, ws[index-1], ws[index]);
c09c9e05ecfa Removed trailing whitespace. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 704 diff changeset	353 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	354 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	355
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	356 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	357 }
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	358 } // class Row
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	359
1838 853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	360 /** Rows in table. */
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	361 protected List<Row> rows;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	362
1838 853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	363 /** Columns in table. */
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	364 protected Column [] columns;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	365
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	366 public WstValueTable() {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	367 rows = new ArrayList<Row>();
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	368 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	369
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	370 public WstValueTable(Column [] columns) {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	371 this();
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	372 this.columns = columns;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	373 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 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 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	376 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	377 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	378 }
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
1838 853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	380 /** Sort rows (by km). */
458 523a256451cd Use unsharp km lookup (epsilon = 0.0001). Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 443 diff changeset	381 public void sortRows() {
523a256451cd Use unsharp km lookup (epsilon = 0.0001). Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 443 diff changeset	382 Collections.sort(rows);
523a256451cd Use unsharp km lookup (epsilon = 0.0001). Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 443 diff changeset	383 }
523a256451cd Use unsharp km lookup (epsilon = 0.0001). Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 443 diff changeset	384
380 5f55047a17e8 Moved expensive allocation of result arrays out of km loop. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 374 diff changeset	385 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	386 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	387 }
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	388
1098 1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	389
1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	390 /**
1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	391 * @param ws (output parameter), gets returned.
1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	392 * @return output parameter ws.
1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	393 */
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	394 public double [] interpolateW(
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	395 double km,
742 c09c9e05ecfa Removed trailing whitespace. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 704 diff changeset	396 double [] qs,
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	397 double [] ws,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	398 Calculation errors
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	399 ) {
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	400 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	401
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	402 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	403
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	404 if (rowIndex >= 0) { // direct row match
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	405 Row row = rows.get(rowIndex);
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	406 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	407 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	408 if (errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	409 // TODO: I18N
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	410 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	411 }
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	412 ws[i] = Double.NaN;
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	413 }
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	414 else {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	415 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	416 && errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	417 // TODO: I18N
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	418 errors.addProblem(
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	419 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	420 }
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	421 }
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	422 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	423 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	424 else { // needs bilinear interpolation
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	425 rowIndex = -rowIndex -1;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	426
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	427 if (rowIndex < 1 \|\| rowIndex >= rows.size()) {
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	428 // do not extrapolate
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	429 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	430 if (errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	431 // TODO: I18N
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	432 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	433 }
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	434 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	435 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	436 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	437 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	438 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	439 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	440 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	441
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	442 return ws;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	443 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	444
1838 853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	445 /**
853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	446 * Interpolate W and Q values at a given km.
853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	447 */
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	448 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	449 return interpolateWQ(km, null);
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	450 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	451
1838 853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	452 /**
853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	453 * Interpolate W and Q values at a given km.
853cd2120d69 Cosmetics, doc. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1721 diff changeset	454 */
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	455 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	456 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	457 }
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	458
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	459 /**
ddd425858169 Improved logging to include exception in fail case of spline interpolation of wst. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 1838 diff changeset	460 * 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	461 */
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	462 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	463
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	464 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	465
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	466 if (rowIndex >= 0) { // direct row match
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	467 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	468 return row.interpolateWQ(steps, this, errors);
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	469 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	470
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	471 rowIndex = -rowIndex -1;
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	472
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	473 if (rowIndex < 1 \|\| rowIndex >= rows.size()) {
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	474 // do not extrapolate
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	475 if (errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	476 // TODO: I18N
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	477 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	478 }
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	479 return new double[2][0];
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	480 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	481
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	482 Row r1 = rows.get(rowIndex-1);
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	483 Row r2 = rows.get(rowIndex);
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	484
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	485 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	486 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	487
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	488 public boolean interpolate(
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	489 double km,
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	490 double [] out,
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	491 QPosition qPosition,
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	492 Function qFunction
649 44175d4720f8 Fix crashes in back jump correction. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 645 diff changeset	493 ) {
427 909196be11a0 First step to calculate "W fuer ungleichwertige Abfluesse" correctly. flys/issue55 Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 426 diff changeset	494 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	495
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	496 out[1] = qFunction.value(getQ(qPosition, km));
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	497
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	498 if (Double.isNaN(out[1])) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	499 return false;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	500 }
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	501
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	502 QPosition nPosition = new QPosition();
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	503 if (getQPosition(km, out[1], nPosition) == null) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	504 return false;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	505 }
427 909196be11a0 First step to calculate "W fuer ungleichwertige Abfluesse" correctly. flys/issue55 Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 426 diff changeset	506
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	507 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	508
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	509 if (rowIndex >= 0) {
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	510 // direct row match
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	511 out[0] = rows.get(rowIndex).getW(nPosition);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	512 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	513 }
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	514
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	515 rowIndex = -rowIndex -1;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	516
1681 e5f7f25a511c fixed flys/issue201 Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1098 diff changeset	517 if (rowIndex < 1 \|\| rowIndex > R1) {
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	518 // do not extrapolate
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	519 return false;
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	520 }
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	521
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	522 Row r1 = rows.get(rowIndex-1);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	523 Row r2 = rows.get(rowIndex);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	524 out[0] = r1.getW(r2, km, nPosition);
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	525
913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	526 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	527 }
909196be11a0 First step to calculate "W fuer ungleichwertige Abfluesse" correctly. flys/issue55 Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 426 diff changeset	528
1098 1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	529
1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	530 /**
1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	531 * Look up interpolation of a Q at given positions.
1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	532 *
1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	533 * @param q the non-interpolated Q value.
1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	534 * @param referenceKm the reference km (e.g. gauge position).
1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	535 * @param kms positions for which to interpolate.
1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	536 * @param ws (output) resulting interpolated ws.
1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	537 * @param qs (output) resulting interpolated qs.
1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	538 * @param errors calculation object to store errors.
1ea7eb72aaa6 Minor cosmetics. Felix Wolfsteller <felix.wolfsteller@intevation.de> parents: 742 diff changeset	539 */
427 909196be11a0 First step to calculate "W fuer ungleichwertige Abfluesse" correctly. flys/issue55 Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 426 diff changeset	540 public QPosition interpolate(
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	541 double q,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	542 double referenceKm,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	543 double [] kms,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	544 double [] ws,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	545 double [] qs,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	546 Calculation errors
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	547 ) {
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	548 return interpolate(
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	549 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	550 }
44175d4720f8 Fix crashes in back jump correction. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 645 diff changeset	551
44175d4720f8 Fix crashes in back jump correction. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 645 diff changeset	552 public QPosition interpolate(
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	553 double q,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	554 double referenceKm,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	555 double [] kms,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	556 double [] ws,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	557 double [] qs,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	558 int startIndex,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	559 int length,
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	560 Calculation errors
649 44175d4720f8 Fix crashes in back jump correction. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 645 diff changeset	561 ) {
645 433f67a076aa Use an explicit reference km for interpolation now. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 633 diff changeset	562 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	563
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	564 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	565 // 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	566 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	567 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	568 if (errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	569 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	570 }
427 909196be11a0 First step to calculate "W fuer ungleichwertige Abfluesse" correctly. flys/issue55 Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 426 diff changeset	571 return null;
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	572 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	573
645 433f67a076aa Use an explicit reference km for interpolation now. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 633 diff changeset	574 Row kmKey = new Row();
433f67a076aa Use an explicit reference km for interpolation now. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 633 diff changeset	575
433f67a076aa Use an explicit reference km for interpolation now. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 633 diff changeset	576 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	577
649 44175d4720f8 Fix crashes in back jump correction. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 645 diff changeset	578 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	579
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	580 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	581 if (errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	582 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	583 }
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	584 ws[i] = Double.NaN;
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	585 continue;
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	586 }
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	587
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	588 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	589 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	590
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	591 if (rowIndex >= 0) {
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	592 // direct row match
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	593 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	594 && errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	595 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	596 }
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	597 continue;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	598 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	599
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	600 rowIndex = -rowIndex -1;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	601
1681 e5f7f25a511c fixed flys/issue201 Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 1098 diff changeset	602 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	603 // do not extrapolate
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	604 if (errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	605 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	606 }
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	607 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	608 continue;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	609 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	610 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	611 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	612
686 3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	613 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	614 && errors != null) {
3dc61e00385e Merged with trunk and introduced hashing of computed values. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 655 diff changeset	615 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	616 }
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	617 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	618
427 909196be11a0 First step to calculate "W fuer ungleichwertige Abfluesse" correctly. flys/issue55 Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 426 diff changeset	619 return qPosition;
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	620 }
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	621
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	622 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	623 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	624 }
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	625
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	626 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	627
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	628 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	629 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	630 }
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	631
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	632 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	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	634 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	635 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	636 }
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	637
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	638 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	639 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	640 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	641 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	642 }
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	643
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	644 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	645 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	646 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	647
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	648 if (q > qMin && q < qMax) {
655 913b52064449 Refactored version of "Berechnung 4" Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 649 diff changeset	649 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	650 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	651 }
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	652 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	653 }
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 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	656 }
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	657
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	658 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	659 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	660 }
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	661
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	662 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	663 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	664 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	665
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	666 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	667 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	668 }
d08f77e7f7e8 WST value table: Qs are now stored in ranges for each column. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 632 diff changeset	669 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	670 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	671 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	672 }
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	673 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	674 // 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 @ 1860:ddd425858169