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 @ 426:c6b7ca0febcb

Fix step to fix issue69: Do the right calculation flys-artifacts/trunk@1924 c6561f87-3c4e-4783-a992-168aeb5c3f6f

author	Sascha L. Teichmann <sascha.teichmann@intevation.de>
date	Sun, 15 May 2011 10:59:08 +0000
parents	aae8f327425e
children	909196be11a0

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
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	5 import de.intevation.flys.model.River;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	6 import de.intevation.flys.model.Wst;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	7 import de.intevation.flys.model.WstColumn;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	8
374 91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	9 import de.intevation.flys.artifacts.cache.CacheFactory;
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	10
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	11 import de.intevation.flys.backend.SessionHolder;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	12
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	13 import java.util.Arrays;
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	14 import java.util.ArrayList;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	15 import java.util.Comparator;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	16 import java.util.List;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	17 import java.util.Collections;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	18
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	19 import org.apache.log4j.Logger;
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	20
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	21 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	22
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	23 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	24
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	25 import org.apache.commons.math.ArgumentOutsideDomainException;
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	26
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	27 import org.hibernate.Session;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	28 import org.hibernate.Query;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	29 import org.hibernate.SQLQuery;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	30
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	31 import org.hibernate.type.StandardBasicTypes;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	32
374 91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	33 import net.sf.ehcache.Cache;
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	34 import net.sf.ehcache.Element;
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	35
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	36 public class WstValueTable
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	37 implements Serializable
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	38 {
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	39 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	40
374 91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	41 public static final String CACHE_NAME = "wst-value-table";
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	42
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	43 // TODO: put this into a property file
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	44 public static final String SQL_POS_WQ =
327 e09634fbf6bc Moved the query complexity of creating a w/q value table into view Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 326 diff changeset	45 "SELECT position, w, q, column_pos" +
e09634fbf6bc Moved the query complexity of creating a w/q value table into view Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 326 diff changeset	46 " FROM wst_value_table" +
e09634fbf6bc Moved the query complexity of creating a w/q value table into view Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 326 diff changeset	47 " WHERE wst_id = :wst_id";
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	48
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	49 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	50
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	51 public static class Column
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	52 implements Serializable
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 protected String name;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	55
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	56 public Column() {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	57 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	58
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	59 public Column(String name) {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	60 this.name = name;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	61 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	62
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	63 public String getName() {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	64 return name;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	65 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	66
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	67 public void setName(String name) {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	68 this.name = name;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	69 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	70 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	71 // class Column
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	72
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	73 public static class QPosition {
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 protected double q;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	76 protected int index;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	77 protected double weight;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	78
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	79 public QPosition() {
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	80 }
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 public QPosition(double q, int index, double weight) {
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	83 this.index = index;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	84 this.weight = weight;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	85 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	86
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	87 } // class Position
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	88
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	89 public static class Row
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	90 implements Serializable, Comparable<Row>
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	91 {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	92 double km;
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	93 double [] ws;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	94 double [] qs;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	95 boolean qSorted;
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 public Row() {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	98 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	99
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	100 public Row(double km) {
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	101 this.km = km;
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	102 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	103
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	104 public Row(double km, double [] ws, double [] qs) {
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	105 this(km);
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	106 this.ws = ws;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	107 this.qs = qs;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	108 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	109
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	110 public static final double linear(
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	111 double x,
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	112 double x1, double x2,
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	113 double y1, double y2
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	114 ) {
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	115 // y1 = m*x1 + b
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	116 // y2 = m*x2 + b
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	117 // y2 - y1 = m*(x2 - x1)
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	118 // m = (y2 - y1)/(x2 - x1) # x2 != x1
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	119 // b = y1 - m*x1
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	120
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	121 if (x1 == x2) {
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	122 return 0.5*(y1 + y2);
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	123 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	124 double m = (y2 - y1)/(x2 - x1);
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	125 double b = y1 - m*x1;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	126 return x*m + b;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	127 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	128
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	129 public static final double factor(double x, double p1, double p2) {
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	130 // 0 = mp1 + b <=> b = -mp1
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	131 // 1 = m*p2 + b
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	132 // 1 = m*(p2 - p1)
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	133 // m = 1/(p2 - p1) # p1 != p2
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	134 // f(x) = x/(p2-p1) - p1/(p2-p1) <=> (x-p1)/(p2-p1)
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	135
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	136 return p1 == p2 ? 0.0 : (x-p1)/(p2-p1);
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	137 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	138
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	139 public static final double weight(double factor, double a, double b) {
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	140 return (1.0-factor)a + factorb;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	141 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	142
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	143 public double getWForKM(Row other, int index, double km) {
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	144 double w1 = ws[index];
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	145 double w2 = other.ws[index];
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	146 double km1 = this.km;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	147 double km2 = other.km;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	148 return linear(km, km1, km2, w1, w2);
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	149 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	150
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	151 public void interpolateW(
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	152 Row other,
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	153 double km,
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	154 double [] input,
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	155 double [] output
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	156 ) {
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	157 double factor = factor(km, this.km, other.km);
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	158
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	159 for (int i = 0; i < input.length; ++i) {
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	160 double q = input[i];
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	161 int idx1 = getQIndex(q);
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	162 int idx2 = other.getQIndex(q);
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	163
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	164 double w1 = idx1 >= 0
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	165 ? ws[idx1]
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	166 : interpolateW(-idx1-1, q);
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	167
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	168 double w2 = idx2 >= 0
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	169 ? other.ws[idx2]
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	170 : other.interpolateW(-idx2-1, q);
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	171
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	172 output[i] = weight(factor, w1, w2);
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	173 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	174 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	175
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	176 public double interpolateW(int idx, double q) {
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	177 return idx < 1 \|\| idx >= qs.length
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	178 ? Double.NaN // do not extrapolate
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	179 : linear(q, qs[idx-1], qs[idx], ws[idx-1], ws[idx]);
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	180 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	181
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	182 public double interpolateW(double q) {
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	183 if (Double.isNaN(q)) return Double.NaN;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	184 int index = getQIndex(q);
360 3571357c85a7 Bugfix: Removed endless loop and a bug while iterating over Hibernate results. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 344 diff changeset	185 return index >= 0 ? ws[index] : interpolateW(-index -1, q);
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	186 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	187
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	188 public int getQIndex(double q) {
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	189 return qSorted ? binaryQIndex(q) : linearQIndex(q);
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	190 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	191
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	192 public int binaryQIndex(double q) {
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	193 return Arrays.binarySearch(qs, q);
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	194 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	195
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	196 public int linearQIndex(double q) {
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	197 switch (qs.length) {
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	198 case 0: break;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	199 case 1:
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	200 if (qs[0] == q) return 0;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	201 if (qs[0] < q) return -(1+1);
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	202 return -(0+1);
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	203 default:
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	204 for (int i = 1; i < qs.length; ++i) {
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	205 double qa = qs[i-1];
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	206 double qb = qs[i];
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	207 if (qa == q) return i-1;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	208 if (qb == q) return i;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	209 if (qa > qb) { double t = qa; qa = qb; qb = t; }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	210 if (q > qa && q < qb) return -(i+1);
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	211 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	212 return -qs.length;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	213 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	214
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	215 return -1;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	216 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	217
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	218 public int compareTo(Row other) {
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	219 double d = km - other.km;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	220 if (d < 0.0) return -1;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	221 if (d > 0.0) return +1;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	222 return 0;
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	223 }
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	224
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	225 public double [][] cloneWQs() {
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	226 return new double [][] {
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	227 (double [])ws.clone(),
344 79401797f4e1 Corrected a silly c&p mistake. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 339 diff changeset	228 (double [])qs.clone() };
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	229 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	230
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	231 public double [][] interpolateWQ(Row other, double km, int steps) {
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	232
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	233 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	234
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	235 if (W < 1) {
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	236 return new double[2][0];
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	237 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	238
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	239 double factor = factor(km, this.km, other.km);
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	240
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	241 double [] splineQ = new double[W];
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	242 double [] splineW = new double[W];
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	243
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	244 double minQ = Double.MAX_VALUE;
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	245 double maxQ = -Double.MAX_VALUE;
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	246
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	247 for (int i = 0; i < W; ++i) {
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	248 double wws = weight(factor, ws[i], other.ws[i]);
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	249 double wqs = weight(factor, qs[i], other.qs[i]);
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	250 splineW[i] = wws;
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	251 splineQ[i] = wqs;
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	252 if (wqs < minQ) minQ = wqs;
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	253 if (wqs > maxQ) maxQ = wqs;
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
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	256 double stepWidth = (maxQ - minQ)/steps;
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	257
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	258 SplineInterpolator interpolator = new SplineInterpolator();
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	259 PolynomialSplineFunction spline =
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	260 interpolator.interpolate(splineQ, splineW);
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	261
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	262 double [] outWs = new double[steps];
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	263 double [] outQs = new double[steps];
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	264
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	265 try {
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	266 double q = minQ;
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	267 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	268 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	269 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	270 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	271 catch (ArgumentOutsideDomainException aode) {
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	272 log.error("Spline interpolation failed.", aode);
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	273 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	274
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	275 return new double [][] { outWs, outQs };
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	276 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	277
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	278 public double [][] interpolateWQ(int steps) {
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	279
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	280 int W = ws.length;
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	281
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	282 if (W < 1) {
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	283 return new double[2][0];
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	284 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	285
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	286 double [] splineW = new double[W];
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	287 double [] splineQ = new double[W];
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	288
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	289 double minQ = Double.MAX_VALUE;
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	290 double maxQ = -Double.MAX_VALUE;
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	291
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	292 for (int i = 0; i < W; ++i) {
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	293 splineW[i] = ws[i];
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	294 splineQ[i] = qs[i];
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	295 if (qs[i] < minQ) minQ = qs[i];
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	296 if (qs[i] > maxQ) maxQ = qs[i];
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	297 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	298
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	299 double stepWidth = (maxQ - minQ)/steps;
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	300
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	301 SplineInterpolator interpolator = new SplineInterpolator();
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	302
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	303 PolynomialSplineFunction spline =
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	304 interpolator.interpolate(splineQ, splineW);
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	305
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	306 double [] outWs = new double[steps];
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	307 double [] outQs = new double[steps];
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	308
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	309 try {
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	310 double q = minQ;
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	311 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	312 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	313 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	314 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	315 catch (ArgumentOutsideDomainException aode) {
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	316 log.error("Spline interpolation failed.", aode);
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	317 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	318
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	319 return new double [][] { outWs, outQs };
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 public int ascendingWs() {
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	323 if (ws.length < 2) {
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	324 return ws.length;
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	325 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	326
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	327 int idx = 1;
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	328
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	329 for (; idx < ws.length; ++idx) {
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	330 if (Double.isNaN(ws[idx]) \|\| ws[idx] < ws[idx-1]) {
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	331 return idx;
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	332 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	333 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	334
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	335 return idx;
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	336 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	337
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	338 public double [][] weightWQs(Row other, double km) {
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	339 int W = Math.min(ws.length, other.ws.length);
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	340 double factor = factor(km, this.km, other.km);
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	341
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	342 double [] outWs = new double[W];
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	343 double [] outQs = new double[W];
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	344
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	345 for (int i = 0; i < W; ++i) {
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	346 outWs[i] = weight(factor, ws[i], other.ws[i]);
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	347 outQs[i] = weight(factor, qs[i], other.qs[i]);
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	348 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	349
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	350 return new double [][] { outWs, outQs };
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	351 }
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	352
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	353 public QPosition getQPosition(double q) {
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	354 int qIndex = getQIndex(q);
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	355
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	356 if (qIndex >= 0) {
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	357 // direct column match
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	358 return new QPosition(q, qIndex, 1.0);
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	359 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	360
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	361 qIndex = -qIndex -1;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	362
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	363 if (qIndex < 0 \|\| qIndex >= qs.length-1) {
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	364 // do not extrapolate
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	365 return null;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	366 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	367
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	368 return new QPosition(
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	369 q, qIndex, factor(q, qs[qIndex-1], qs[qIndex]));
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	370 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	371
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	372 public QPosition getQPosition(Row other, double km, double q) {
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	373
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	374 QPosition qpt = getQPosition(q);
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	375 QPosition qpo = other.getQPosition(q);
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	376
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	377 if (qpt == null \|\| qpo == null) {
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	378 return null;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	379 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	380
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	381 double kmWeight = factor(km, this.km, other.km);
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	382
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	383 // XXX: Index interpolation is a bit sticky here!
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	384
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	385 int index = (int)Math.round(
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	386 weight(kmWeight, qpt.index, qpo.index));
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	387
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	388 double weight = weight(kmWeight, qpt.weight, qpo.weight);
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	389
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	390 return new QPosition(q, index, weight);
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	391 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	392
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	393 public void storeWQ(
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	394 QPosition qPosition,
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	395 double [] ows,
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	396 double [] oqs,
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	397 int index
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	398 ) {
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	399 int qIdx = qPosition.index;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	400 double qWeight = qPosition.weight;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	401
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	402 if (qWeight == 1.0) {
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	403 oqs[index] = qs[qIdx];
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	404 ows[index] = ws[qIdx];
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	405 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	406 else {
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	407 oqs[index] = weight(qWeight, qs[qIdx-1], qs[qIdx]);
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	408 ows[index] = weight(qWeight, ws[qIdx-1], ws[qIdx]);
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	409 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	410 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	411
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	412 public void storeWQ(
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	413 Row other,
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	414 double km,
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	415 QPosition qPosition,
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	416 double [] ows,
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	417 double [] oqs,
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	418 int index
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	419 ) {
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	420 double kmWeight = factor(km, this.km, other.km);
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	421
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	422 double tq, tw;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	423 double oq, ow;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	424
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	425 int qIdx = qPosition.index;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	426 double qWeight = qPosition.weight;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	427
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	428 if (qWeight == 1.0) {
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	429 tq = qs[qIdx];
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	430 tw = ws[qIdx];
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	431 oq = other.qs[qIdx];
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	432 ow = other.ws[qIdx];
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	433 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	434 else {
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	435 tq = weight(qWeight, qs[qIdx-1], qs[qIdx]);
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	436 tw = weight(qWeight, ws[qIdx-1], ws[qIdx]);
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	437 oq = weight(qWeight, other.qs[qIdx-1], other.qs[qIdx]);
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	438 ow = weight(qWeight, other.ws[qIdx-1], other.ws[qIdx]);
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	439 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	440
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	441 oqs[index] = weight(kmWeight, oq, tq);
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	442 ows[index] = weight(kmWeight, ow, tw);
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	443 }
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	444 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	445 // class Row
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	446
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	447 protected List<Row> rows;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	448
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	449 protected Column [] columns;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	450
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	451 public WstValueTable() {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	452 rows = new ArrayList<Row>();
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	453 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	454
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	455 public WstValueTable(Column [] columns) {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	456 this();
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	457 this.columns = columns;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	458 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	459
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	460 public double [] interpolateW(double km, double [] qs) {
380 5f55047a17e8 Moved expensive allocation of result arrays out of km loop. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 374 diff changeset	461 return interpolateW(km, qs, new double[qs.length]);
5f55047a17e8 Moved expensive allocation of result arrays out of km loop. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 374 diff changeset	462 }
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	463
380 5f55047a17e8 Moved expensive allocation of result arrays out of km loop. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 374 diff changeset	464 public double [] interpolateW(double km, double [] qs, double [] ws) {
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	465
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	466 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	467
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	468 if (rowIndex >= 0) { // direct row match
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	469 Row row = rows.get(rowIndex);
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	470 for (int i = 0; i < qs.length; ++i) {
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	471 ws[i] = row.interpolateW(qs[i]);
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	472 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	473 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	474 else { // needs bilinear interpolation
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	475 rowIndex = -rowIndex -1;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	476
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	477 if (rowIndex < 1 \|\| rowIndex >= rows.size()) {
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	478 // do not extrapolate
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	479 Arrays.fill(ws, Double.NaN);
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	480 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	481 else {
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	482 rows.get(rowIndex-1).interpolateW(
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	483 rows.get(rowIndex),
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	484 km, qs, ws);
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	485 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	486 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	487
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	488 return ws;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	489 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	490
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	491
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	492 public double [][] interpolateWQ(double km) {
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	493 return interpolateWQ(km, DEFAULT_Q_STEPS);
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	494 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	495
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	496 public double [][] interpolateWQ(double km, int steps) {
aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	497
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	498 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	499
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	500 if (rowIndex >= 0) { // direct row match
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	501 Row row = rows.get(rowIndex);
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	502 return row.interpolateWQ(steps);
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	503 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	504
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	505 rowIndex = -rowIndex -1;
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	506
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	507 if (rowIndex < 1 \|\| rowIndex >= rows.size()) {
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	508 // do not extrapolate
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	509 return new double[2][0];
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	510 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	511
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	512 Row r1 = rows.get(rowIndex-1);
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	513 Row r2 = rows.get(rowIndex);
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	514
395 aae8f327425e 'berechnete Abflusskurve': generate ws from qs Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 380 diff changeset	515 return r1.interpolateWQ(r2, km, steps);
339 4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	516 }
4509ba8fae68 Implementation of "Abflusskurve/Abflusstafel" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 335 diff changeset	517
426 c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	518 public boolean interpolate(
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	519 double q,
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	520 int referenceIndex,
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	521 double [] kms,
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	522 double [] ws,
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	523 double [] qs
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	524 ) {
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	525 Row kmKey = new Row(kms[referenceIndex]);
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	526
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	527 int rowIndex = Collections.binarySearch(rows, kmKey);
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	528
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	529 int R1 = rows.size()-1;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	530
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	531 QPosition qPosition = null;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	532
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	533 if (rowIndex >= 0) { // direct row match
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	534 qPosition = rows.get(rowIndex).getQPosition(q);
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	535 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	536 else {
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	537 rowIndex = -rowIndex -1;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	538
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	539 if (rowIndex < 1 \|\| rowIndex >= R1) {
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	540 // do not extrapolate
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	541 return false;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	542 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	543
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	544 // interpolation case
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	545 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	546 Row r2 = rows.get(rowIndex);
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	547
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	548 qPosition = r1.getQPosition(r2, kms[referenceIndex], q);
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	549 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	550
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	551 if (qPosition == null) {
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	552 // we cannot locate q in reference row
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	553 return false;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	554 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	555
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	556 for (int i = 0; i < kms.length; ++i) {
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	557 kmKey.km = kms[i];
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	558 rowIndex = Collections.binarySearch(rows, kmKey);
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	559
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	560 if (rowIndex >= 0) {
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	561 // direct row match
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	562 rows.get(rowIndex).storeWQ(qPosition, ws, qs, i);
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	563 continue;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	564 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	565
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	566 rowIndex = -rowIndex -1;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	567
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	568 if (rowIndex < 1 \|\| rowIndex >= R1) {
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	569 // do not extrapolate
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	570 ws[i] = qs[i] = Double.NaN;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	571 continue;
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 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	574 Row r2 = rows.get(rowIndex);
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	575 r1.storeWQ(r2, kms[i], qPosition, ws, qs, i);
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	576 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	577
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	578 return true;
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	579 }
c6b7ca0febcb Fix step to fix issue69: Do the right calculation Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 395 diff changeset	580
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	581 public static WstValueTable getTable(River river) {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	582 return getTable(river, 0);
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	583 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	584
374 91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	585 public static final class CacheKey
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	586 implements Serializable
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	587 {
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	588 private int riverId;
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	589 private int kind;
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	590
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	591 public CacheKey(int riverId, int kind) {
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	592 this.riverId = riverId;
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	593 this.kind = kind;
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	594 }
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	595
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	596 public int hashCode() {
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	597 return (riverId << 8) \| kind;
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	598 }
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	599
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	600 public boolean equals(Object other) {
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	601 if (!(other instanceof CacheKey)) {
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	602 return false;
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	603 }
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	604 CacheKey o = (CacheKey)other;
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	605 return riverId == o.riverId && kind == o.kind;
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	606 }
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	607 } // class CacheKey
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	608
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	609 public static WstValueTable getTable(River river, int kind) {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	610
374 91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	611 Cache cache = CacheFactory.getCache(CACHE_NAME);
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	612
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	613 CacheKey cacheKey;
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	614
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	615 if (cache != null) {
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	616 cacheKey = new CacheKey(river.getId(), kind);
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	617 Element element = cache.get(cacheKey);
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	618 if (element != null) {
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	619 log.debug("got wst value table from cache");
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	620 return (WstValueTable)element.getValue();
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	621 }
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	622 }
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	623 else {
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	624 cacheKey = null;
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	625 }
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	626
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	627 WstValueTable valueTable = getTableUncached(river, kind);
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	628
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	629 if (cacheKey != null) {
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	630 log.debug("store wst value in cache");
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	631 Element element = new Element(cacheKey, valueTable);
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	632 cache.put(element);
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	633 }
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	634
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	635 return valueTable;
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	636 }
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	637
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	638 public static WstValueTable getTableUncached(River river, int kind) {
91fbaa2744bf Added caching support for the wst value tables. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 372 diff changeset	639
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	640 Session session = SessionHolder.HOLDER.get();
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	641
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	642 Query query = session.createQuery(
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	643 "from Wst where river=:river and kind=:kind");
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	644 query.setParameter("river", river);
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	645 query.setInteger("kind", kind);
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	646
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	647 List<Wst> wsts = query.list();
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	648
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	649 if (wsts.isEmpty()) {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	650 return null;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	651 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	652
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	653 Wst wst = wsts.get(0);
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	654
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	655 // TODO: Do this sorting at database level
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	656 List<WstColumn> wstColumns = new ArrayList(wst.getColumns());
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	657 Collections.sort(wstColumns, new Comparator<WstColumn>() {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	658 public int compare(WstColumn a, WstColumn b) {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	659 int pa = a.getPosition();
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	660 int pb = b.getPosition();
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	661 if (pa < pb) return -1;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	662 if (pa > pb) return +1;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	663 return 0;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	664 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	665 });
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	666
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	667 Column [] columns = new Column[wstColumns.size()];
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	668 for (int i = 0; i < columns.length; ++i) {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	669 columns[i] = new Column(wstColumns.get(i).getName());
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	670 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	671
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	672 // using native SQL here to avoid myriad of small objects.
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	673 SQLQuery sqlQuery = session.createSQLQuery(SQL_POS_WQ)
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	674 .addScalar("position", StandardBasicTypes.DOUBLE)
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	675 .addScalar("w", StandardBasicTypes.DOUBLE)
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	676 .addScalar("q", StandardBasicTypes.DOUBLE)
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	677 .addScalar("column_pos", StandardBasicTypes.INTEGER);
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	678
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	679 sqlQuery.setInteger("wst_id", wst.getId());
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	680
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	681 WstValueTable valueTable = new WstValueTable(columns);
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	682
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	683 int lastColumnNo = -1;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	684 Row row = null;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	685
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	686 Double lastQ = -Double.MAX_VALUE;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	687 boolean qSorted = true;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	688
360 3571357c85a7 Bugfix: Removed endless loop and a bug while iterating over Hibernate results. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 344 diff changeset	689 for (Object r: sqlQuery.list()) {
3571357c85a7 Bugfix: Removed endless loop and a bug while iterating over Hibernate results. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 344 diff changeset	690 Object[] result = (Object[]) r;
3571357c85a7 Bugfix: Removed endless loop and a bug while iterating over Hibernate results. Ingo Weinzierl <ingo.weinzierl@intevation.de> parents: 344 diff changeset	691
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	692 double km = (Double) result[0];
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	693 Double w = (Double) result[1];
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	694 Double q = (Double) result[2];
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	695 int columnNo = (Integer)result[3];
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	696
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	697 if (columnNo > lastColumnNo) { // new row
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	698 if (row != null) {
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	699 row.qSorted = qSorted;
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	700 valueTable.rows.add(row);
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	701 }
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	702 row = new Row(
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	703 km,
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	704 new double[columnNo+1],
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	705 new double[columnNo+1]);
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	706 lastQ = -Double.MAX_VALUE;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	707 qSorted = true;
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	708 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	709
335 64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	710 row.ws[columnNo] = w != null ? w : Double.NaN;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	711 row.qs[columnNo] = q != null ? q : Double.NaN;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	712
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	713 if (qSorted && (q == null \|\| lastQ > q)) {
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	714 qSorted = false;
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	715 }
64cfbd631f29 Implemented the "Wasserstand/Wasserspiegellagen" calculation. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: 327 diff changeset	716 lastQ = q;
326 07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	717
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	718 lastColumnNo = columnNo;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	719 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	720
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	721 if (row != null) {
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	722 valueTable.rows.add(row);
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	723 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	724
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	725 return valueTable;
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	726 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	727 }
07eba6d3b4a5 Fetches w/q value tables from the backend. Sascha L. Teichmann <sascha.teichmann@intevation.de> parents: diff changeset	728 // 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 @ 426:c6b7ca0febcb