Commit 1d55470a authored by Christoph Junghans's avatar Christoph Junghans

dpot_crop.pl: obsolete

parent 59233ebe
......@@ -101,7 +101,6 @@ convergence_check default convergence_check_default.sh
# table preparation
resample target resample_target.sh
dpot crop dpot_crop.pl
# ibi specific stuff
update ibi update_ibi.sh
......
#! /usr/bin/perl -w
#
# Copyright 2009-2011 The VOTCA Development Team (http://www.votca.org)
#
# Licensed under the Apache License, Version 2.0 (the "License");
# you may not use this file except in compliance with the License.
# You may obtain a copy of the License at
#
# http://www.apache.org/licenses/LICENSE-2.0
#
# Unless required by applicable law or agreed to in writing, software
# distributed under the License is distributed on an "AS IS" BASIS,
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
# See the License for the specific language governing permissions and
# limitations under the License.
#
use strict;
$_=$0;
s#^.*/##;
my $progname=$_;
my $usage="Usage: $progname [OPTIONS] <file> <a> <b>";
#Defaults
my $withflag=undef;
while ((defined ($ARGV[0])) and ($ARGV[0] =~ /^-./))
{
if (($ARGV[0] !~ /^--/) and (length($ARGV[0])>2)){
$_=shift(@ARGV);
#short opt having agruments examples fo
if ( $_ =~ /^-[fo]/ ) {
unshift(@ARGV,substr($_,0,2),substr($_,2));
}
else{
unshift(@ARGV,substr($_,0,2),"-".substr($_,2));
}
}
if (($ARGV[0] eq "-h") or ($ARGV[0] eq "--help"))
{
print <<END;
$progname, version %version%
crop the potential update at poorly sampled ends
$usage
Allowed options:
-h, --help Show this help message
Examples:
* $progname tmp.dpot.cur tmp.dpot.new
END
exit;
}
elsif ($ARGV[0] eq "--withflag")
{
shift(@ARGV);
die "nothing given for --withflag" unless $#ARGV > -1;
$withflag = $ARGV[0];
}
else
{
die "Unknown option '".$ARGV[0]."' !\n";
}
shift(@ARGV);
}
#Print usage
die "missing parameters\n$usage\n" unless $#ARGV >= 1;
use CsgFunctions;
my $infile="$ARGV[0]";
my $outfile="$ARGV[1]";
my @r;
my @val;
my @flag;
(readin_table($infile,@r,@val,@flag)) || die "$progname: error at readin_table\n";
# find last u/o
my $i_first;
# TODO: look for at least 3 successive points with i
for($i_first=0; ($i_first<$#r) && ($flag[$i_first] =~ /[uo]/); $i_first++) {}
my $ncrop=0;
while($i_first + $ncrop<=$#r-3) {
my $i = $i_first + $ncrop;
my $delta_1 = $val[$i] - $val[$i + 1];
my $delta_2 = $val[$i + 1 ] - $val[$i + 2];
# do both deltas have the same sign?
if($delta_1 * $delta_2 > 0) {
last;
} elsif (abs($val[$i]) < 0.5 && abs($val[$i+1]) < 0.5) {
last;
}
$flag[$i]='o';
$ncrop++;
if($ncrop > 3) {
print "error: need to crop more than 3 points in $infile. think about sampleing/grid interval.";
exit 1;
}
}
if($ncrop > 0) {
print "warnng, I cropped $ncrop points at the beginning\n";
}
saveto_table($outfile,@r,@val,@flag) || die "$progname: error at save table\n";
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment