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