Mercurial > dive4elements > river
view flys-artifacts/src/main/java/de/intevation/flys/artifacts/model/WQCKms.java @ 4798:39885bdfc6fc
Added calculation of the "Umhuellende" to calculation of "W fuer ungleichwertige Abfluesse".
This is done by figuring out the WST columns that imfold the data and then do
simple "gleichwertige" calculations from the start of the interval. This
is too much because only the Qs are needed for the "Umhuellende".
author | Sascha L. Teichmann <teichmann@intevation.de> |
---|---|
date | Sun, 13 Jan 2013 16:18:28 +0100 |
parents | 5642a83420f2 |
children | bcf25d8c183e |
line wrap: on
line source
package de.intevation.flys.artifacts.model; import gnu.trove.TDoubleArrayList; /** * This class represents a pool of data triples that consists of 'W', 'Q' and * 'KM' data with corrected 'W' values computed by a BackJumpCorrector. * * @author <a href="mailto:ingo.weinzierl@intevation.de">Ingo Weinzierl</a> */ public class WQCKms extends WQKms { protected TDoubleArrayList cws; public WQCKms() { } public WQCKms(WQKms other, double [] cws) { this.ws = other.ws; this.qs = other.qs; this.kms = other.kms; this.cws = new TDoubleArrayList(cws); } public WQCKms(double[] kms, double[] qs, double[] ws, double[] cws) { super(kms, qs, ws); this.cws = new TDoubleArrayList(cws); } @Override public void removeNaNs() { removeNaNs(new TDoubleArrayList [] { ws, qs, cws, kms }); } /** * Adds a new row to this data pool with corrected W. * * @param w a W. * @param q a Q. * @param kms a Kms. * @param cw The corrected W. */ public void add(double w, double q, double kms, double cw) { super.add(w, q, kms); cws.add(cw); } @Override public double[] get(int idx) { return get(idx, new double[4]); } /** * This method returns a 4dim array of W, Q,Kms and corrected W. * * @param idx The position of the triple. * @param dst destination array * * @return a 4dim array of [W, Q, Kms, CW] in dst. */ @Override public double[] get(int idx, double[] dst) { dst = super.get(idx, dst); if (dst.length < 4) { return dst; } if (cws != null && cws.size() > idx) { dst[3] = cws.getQuick(idx); } return dst; } public double getC(int idx) { return cws.getQuick(idx); } /** * Returns the double array of corrected W values. * * @return the double array of corrected W values. */ public double[] getCWs() { return cws.toNativeArray(); } } // vim:set ts=4 sw=4 si et sta sts=4 fenc=utf8 :