makemesh

   1 #!/usr/bin/perl -w
   2
   3 #=======================================================================
   4 # $Id$
   5 # Create random 3D surface for use in plotting programs.
   6 #
   7 # Character set: UTF-8
   8 # ©opyleft 2007– Øyvind A. Holm <sunny@sunbase.org>
   9 # License: GNU General Public License version 2 or later, see end of
  10 # file for legal stuff.
  11 #=======================================================================
  12
  13 BEGIN {
  14     our @version_array;
  15 }
  16
  17 use strict;
  18 use Getopt::Long;
  19
  20 $| = 1;
  21
  22 our $Debug = 0;
  23
  24 our %Opt = (
  25     'debug' => 0,
  26     'help' => 0,
  27     'max' => 1000,
  28     'min' => 0,
  29     'verbose' => 0,
  30     'version' => 0,
  31 );
  32
  33 our $progname = $0;
  34 $progname =~ s/^.*\/(.*?)$/$1/;
  35
  36 my $rcs_id = '$Id$';
  37 my $id_date = $rcs_id;
  38 $id_date =~ s/^.*?\d+ (\d\d\d\d-.*?\d\d:\d\d:\d\d\S+).*/$1/;
  39
  40 push(@main::version_array, $rcs_id);
  41
  42 Getopt::Long::Configure("bundling");
  43 GetOptions(
  44     "debug" => \$Opt{'debug'},
  45     "help|h" => \$Opt{'help'},
  46     "max=f" => \$Opt{'max'},
  47     "min=f" => \$Opt{'min'},
  48     "verbose|v+" => \$Opt{'verbose'},
  49     "version" => \$Opt{'version'},
  50 ) || die("$progname: Option error. Use -h for help.\n");
  51
  52 $Opt{'debug'} && ($Debug = 1);
  53 $Opt{'help'} && usage(0);
  54 $Opt{'version'} && print_version();
  55
  56 my ($Lat, $Lon) =
  57     (  0,    0);
  58
  59 my @Ele = ();
  60
  61 my $Size = 20;
  62
  63 while (1) {
  64     D("Lon = '$Lon'");
  65     for $Lon (0..$Size) {
  66         my $Left = defined($Ele[$Lat][$Lon-1]) ? $Ele[$Lat][$Lon-1]  : rand(1);
  67         my $Right = defined($Ele[$Lat][$Lon+1]) ? $Ele[$Lat][$Lon+1] : rand(1);
  68         my $Over = defined($Ele[$Lat+1][$Lon]) ? $Ele[$Lat+1][$Lon]  : rand(1);
  69         my $Under = defined($Ele[$Lat-1][$Lon]) ? $Ele[$Lat-1][$Lon] : rand(1);
  70         my $Curr = ($Left+$Right+$Over+$Under)/4;
  71         $Curr += rand(1) >= 0.5 ? rand(1) : 0-rand(1);
  72         # until ($Curr >= $Opt{'min'}) {
  73         #     D("Curr >= Min");
  74         #     $Curr += rand(1) * 1.0;
  75         # }
  76         # until ($Curr <= $Opt{'max'}) {
  77         #     D("Curr <= Max");
  78         #     $Curr -= rand(1) * 1.0;
  79         # }
  80         if ($Lat > 0 && $Lat < $Size && $Lon > 0 && $Lon < $Size) {
  81             # print("$Lat\t$Lon\t$Curr\n");
  82             # print($Lat   . "\t" . $Lon-1 . "\t" . $Left  . "\n");
  83             print($Lat   . "\t" . $Lon   . "\t" . $Curr  . "\n");
  84             # print("\n");
  85             print($Lat+1 . "\t" . $Lon   . "\t" . $Over  . "\n");
  86             print($Lat   . "\t" . $Lon   . "\t" . $Curr  . "\n");
  87             # print("\n");
  88             # print($Lat   . "\t" . $Lon+1 . "\t" . $Right . "\n");
  89             print($Lat   . "\t" . $Lon   . "\t" . $Curr  . "\n");
  90             # print("\n");
  91             print($Lat-1 . "\t" . $Lon   . "\t" . $Under . "\n");
  92             print($Lat   . "\t" . $Lon   . "\t" . $Curr  . "\n");
  93             # print("\n");
  94         }
  95         $Ele[$Lat][$Lon] = $Curr;
  96         print("\n");
  97     }
  98     $Lat += 1;
  99     last if ($Lat > $Size);
 100 }
 101
 102 sub print_version {
 103     # Print program version {{{
 104     for (@main::version_array) {
 105         print("$_\n");
 106     }
 107     exit(0);
 108     # }}}
 109 } # print_version()
 110
 111 sub usage {
 112     # Send the help message to stdout {{{
 113     my $Retval = shift;
 114
 115     print(<<END);
 116
 117 $rcs_id
 118
 119 Usage: $progname [options] [file [files [...]]]
 120
 121 Options:
 122
 123   -h, --help
 124     Show this help.
 125   -v, --verbose
 126     Increase level of verbosity. Can be repeated.
 127   --version
 128     Print version information.
 129   --debug
 130     Print debugging messages.
 131
 132 END
 133     exit($Retval);
 134     # }}}
 135 } # usage()
 136
 137 sub msg {
 138     # Print a status message to stderr based on verbosity level {{{
 139     my ($verbose_level, $Txt) = @_;
 140
 141     if ($Opt{'verbose'} >= $verbose_level) {
 142         print(STDERR "$progname: $Txt\n");
 143     }
 144     # }}}
 145 } # msg()
 146
 147 sub D {
 148     # Print a debugging message {{{
 149     $Debug || return;
 150     my @call_info = caller;
 151     chomp(my $Txt = shift);
 152     my $File = $call_info[1];
 153     $File =~ s#\\#/#g;
 154     $File =~ s#^.*/(.*?)$#$1#;
 155     print(STDERR "$File:$call_info[2] $$ $Txt\n");
 156     return("");
 157     # }}}
 158 } # D()
 159
 160 __END__
 161
 162 # Plain Old Documentation (POD) {{{
 163
 164 =pod
 165
 166 =head1 NAME
 167
 168
 169
 170 =head1 REVISION
 171
 172 $Id$
 173
 174 =head1 SYNOPSIS
 175
 176  [options] [file [files [...]]]
 177
 178 =head1 DESCRIPTION
 179
 180
 181
 182 =head1 OPTIONS
 183
 184 =over 4
 185
 186 =item B<-h>, B<--help>
 187
 188 Print a brief help summary.
 189
 190 =item B<-v>, B<--verbose>
 191
 192 Increase level of verbosity. Can be repeated.
 193
 194 =item B<--version>
 195
 196 Print version information.
 197
 198 =item B<--debug>
 199
 200 Print debugging messages.
 201
 202 =back
 203
 204 =head1 BUGS
 205
 206
 207
 208 =head1 AUTHOR
 209
 210 Made by Øyvind A. Holm S<E<lt>sunny@sunbase.orgE<gt>>.
 211
 212 =head1 COPYRIGHT
 213
 214 Copyleft © Øyvind A. Holm E<lt>sunny@sunbase.orgE<gt>
 215 This is free software; see the file F<COPYING> for legalese stuff.
 216
 217 =head1 LICENCE
 218
 219 This program is free software; you can redistribute it and/or modify it
 220 under the terms of the GNU General Public License as published by the
 221 Free Software Foundation; either version 2 of the License, or (at your
 222 option) any later version.
 223
 224 This program is distributed in the hope that it will be useful, but
 225 WITHOUT ANY WARRANTY; without even the implied warranty of
 226 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
 227 See the GNU General Public License for more details.
 228
 229 You should have received a copy of the GNU General Public License along
 230 with this program; if not, write to the Free Software Foundation, Inc.,
 231 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 232
 233 =head1 SEE ALSO
 234
 235 =cut
 236
 237 # }}}
 238
 239 # vim: set fenc=UTF-8 ft=perl fdm=marker ts=4 sw=4 sts=4 et fo+=w :
 240 # End of file $Id$