lib/nonmem_subs.pm

   1 # {{{ include statements
   2
   3 start include statements
   4 if( $0 =~ /nonmem.pm$/ ) {
   5   use FindBin qw($Bin);
   6   use lib "$Bin";
   7 }
   8 use File::Copy 'cp';
   9 #use ext::Config::Tiny;
  10 use ext::IPC::Run3;
  11 use Config;
  12 if( $0 =~ /nonmem.pm$/ ) {
  13   require PsN;
  14   my $nonmem = nonmem -> new( 'modelfile'       => $ARGV[0],
  15                               'outputfile'      => $ARGV[1],
  16                               'nice'            => $ARGV[2],
  17                               'version'         => $ARGV[3],
  18                               'fsubs'           => [split( /,/ , $ARGV[6] )],
  19                               'nm_directory'    => $ARGV[7],
  20                               'show_version'    => 0);
  21
  22   my $no_remote_compile = $ARGV[4];
  23   my $no_remote_execution = $ARGV[5];
  24
  25   unless( $no_remote_compile ){
  26     unless( $nonmem -> compile() ){
  27       debug -> die( message => $nonmem -> error_message );
  28     }
  29   }
  30
  31   unless( $no_remote_execution ){
  32     $nonmem -> execute;
  33   }
  34 }
  35 end include statements
  36
  37 # }}}
  38
  39 # {{{ new
  40
  41 start new
  42     {
  43       unless( defined $this -> {'outputfile'} ){
  44         $this -> {'outputfile'} = $this -> {'modelfile'};
  45         $this -> {'outputfile'} =~ s/\.mod$/\.lst/;
  46       }
  47
  48       unless( defined $this -> {'nm_directory'} ){
  49         my $nmdir = $PsN::config -> { 'nm_versions' } -> { $this -> {'version'} };
  50
  51         unless( defined $nmdir ){
  52           print "Unknown NONMEM version ",$this -> {'version'}," specified.\n";
  53           print "Configuration file is ",$PsN::config_file,"\n";
  54           my @nmkeys = keys %{$PsN::config -> { 'nm_versions' }};
  55           if ( $#nmkeys == 0 and
  56                defined $nmkeys[0] and
  57                defined $PsN::config -> { 'nm_versions' } -> {$nmkeys[0]} ) {
  58             print "Using the only defined version $nmkeys[0] in ",$PsN::config -> { 'nm_versions' } -> {$nmkeys[0]}," instead\n";
  59             $nmdir = $PsN::config -> { 'nm_versions' } -> {$nmkeys[0]};
  60           }
  61         }
  62
  63         $this -> {'nm_directory'} = $nmdir;
  64
  65         unless( defined $PsN::config -> {'compiler'} and defined $PsN::config -> {'compiler'} -> {'name'} ){
  66           debug -> warn( level => 1,
  67                          message => "No compiler defined, assuming g77" );
  68           $this -> {'compiler'} = 'g77';
  69         } else {
  70           $this -> {'compiler'} = $PsN::config -> {'compiler'} -> {'name'};
  71           $this -> {'compiler_options'} = $PsN::config -> {'compiler'} -> {'options'};
  72         }
  73       }
  74     }
  75 end new
  76
  77 # }}}
  78
  79 # {{{ compile
  80
  81 start compile
  82     {
  83       my $version = $self -> {'version'};
  84       my $nmdir = $self -> {'nm_directory'};
  85       my $includes;
  86       my $nmlink;
  87       my $adaptive = $self -> {'adaptive'} ? '_adaptive' : '';
  88       if( -e "$nmdir/util/nmlink5.exe" ){
  89         $nmlink = "$nmdir/util/nmlink5.exe";
  90         $self -> {'version'} = 5;
  91         if( $self -> {'show_version'} ){
  92           $version = $self -> {'version'};
  93         } else {
  94           $version = '';
  95         }
  96       } elsif ( -e "$nmdir/util/nmlink6.exe" ){
  97         $self -> {'version'} = 6;
  98         if( $self -> {'show_version'} ){
  99           $version = $self -> {'version'};
 100         } else {
 101           $version = '';
 102         }
 103         $nmlink = "$nmdir/util/nmlink6.exe";
 104         if( $Config{osname} eq 'MSWin32' ){
 105           $includes = "-I$nmdir\\sizes";
 106         } else {
 107           $includes = "-I$nmdir/sizes";
 108         }
 109       } else {
 110         my $err_version = ( defined $nmdir and $nmdir ne '' ) ? $nmdir : '[not configured]';
 111         debug -> die( message => "Unable to find a supported version of NONMEM\n".
 112                       "The NONMEM installation directory is $err_version for version [".
 113                       $self -> {'version'}."] according to psn.conf." );
 114       }
 115
 116       # first clean up from old compile
 117       unlink( 'FLIB','FCON', 'FDATA', 'FREPORT','FSUBS', 'FSUBS.f','LINK.LNK','FSTREAM', 'PRDERR', 'nonmem.exe', 'nonmem', 'nonmem5', 'nonmem6', 'nonmem5_adaptive', 'nonmem6_adaptive' );
 118
 119       my $nm;
 120
 121       if( $Config{osname} eq 'MSWin32' ){
 122         $nm="$nmdir\\nm";
 123       } else {
 124         $nm="$nmdir/nm";
 125       }
 126
 127       my $modelfile = $self -> {'modelfile'};
 128
 129       run3( "$nmdir/tr/nmtran.exe", $modelfile, \$self -> {'nmtran_message'}, \$self -> {'nmtran_message'} );
 130
 131       #$self -> {'nmtran_message'} = `$nmdir/tr/nmtran.exe < $modelfile 2>&1`;
 132
 133       open( NMMSG, '>compilation_output.txt' );
 134       print( NMMSG $self -> {'nmtran_message'}, "\n" );
 135
 136       unless(-e 'FREPORT'){
 137         $self -> {'error_message'} = "NMtran failed: \n" . $self -> {'nmtran_message'} ;
 138         return 0;
 139       }
 140
 141       my $nmlink_message;
 142       run3( "$nmlink", undef, \$nmlink_message, \$nmlink_message );
 143
 144       print( NMMSG $nmlink_message, "\n" );
 145
 146       my $fsub;
 147
 148       if(-e 'FSUBS'){
 149         if( $self -> {'compiler'} eq 'g77' ){
 150           rename("FSUBS", "FSUBS.f");
 151           $fsub='FSUBS.f';
 152         } else {
 153           rename("FSUBS", "FSUBS.for");
 154           $fsub='FSUBS.for';
 155         }
 156       }
 157
 158       if( defined $self -> {'fsubs'} ){
 159         foreach my $sub ( @{$self -> {'fsubs'}} ){
 160           $fsub .= " $sub";
 161         }
 162       }
 163
 164       open( FH, "<", 'LINK.LNK' );
 165       my @link;
 166       while( <FH> ){
 167         my $tmp = $_;
 168         $tmp =~ s/^\s+//;
 169         $tmp =~ s/\s+$//;
 170         push(@link, $tmp);
 171       }
 172       close( FH );
 173
 174       my $compile_message;
 175       my $compile_command;
 176
 177       if( $self -> {'compiler'} eq 'g77' ){
 178
 179         my @nmlib = ("$nm/NONMEM.o", "$nm/BLKDAT.o", "$nm/nonmem.a");
 180         $compile_command = "g77 " . $self -> {'compiler_options'} . " -ononmem$version$adaptive $includes $fsub @link @nmlib 2>&1";
 181
 182       } elsif( $self -> {'compiler'} eq 'df' or $self -> {'compiler'} eq 'fl32' ){
 183
 184         my @nmlib = ("$nm\\nonmem.obj", "$nm\\blkdat.obj", "$nm\\nonmem.lib");
 185         $compile_command = $self -> {'compiler'}." " . $self -> {'compiler_options'} . " /Fenonmem $fsub @nmlib @link";
 186
 187       } elsif( $self -> {'compiler'} eq 'ifort' ){
 188         my @nmlib = ("$nm/NONMEM.o", "$nm/BLKDAT.o", "$nm/nonmem.a");
 189         $compile_command = "ifort        " . $self -> {'compiler_options'} . " -ononmem$version$adaptive $includes $fsub @link @nmlib 2>&1";
 190       }
 191
 192
 193       run3( "$compile_command", undef, \$compile_message, \$compile_message );
 194
 195       print( NMMSG $compile_message, "\n" );
 196
 197       close( NMMSG );
 198
 199 #      if( $self -> {'compiler'} eq 'df' or $self -> {'compiler'} eq 'fl32' ){
 200         unless(-e "nonmem.exe" or -e "nonmem$version$adaptive" ) {
 201           $self -> {'error_message'} = "Fortran Compilation failed: \n" . $compile_message ;
 202           return 0;
 203         }
 204 #      } else {
 205 #       unless(-e "nonmem$version") {
 206 #         $self -> {'error_message'} = "Fortran Compilation failed: \n" . $compile_message ;
 207 #         return 0;
 208 #       }
 209 #      }
 210     }
 211 end compile
 212
 213 # }}}
 214
 215 # {{{ execute
 216
 217 start execute
 218     {
 219       my $version = $self -> {'version'};
 220       my $nmdir = $self -> {'nm_directory'};
 221       my $outputfile = $self -> {'outputfile'};
 222       my $adaptive = $self -> {'adaptive'} ? '_adaptive' : '';
 223
 224       my ( $start_time, $fin_time );
 225       if( $Config{osname} eq 'MSWin32' ){
 226         $start_time = `date /T`;
 227         chomp($start_time);
 228         $start_time = $start_time.' '.`time /T`;
 229         run3( "nonmem.exe", undef, \undef );
 230         $fin_time = `date /T`;
 231         chomp($fin_time);
 232         $fin_time = $fin_time.' '.`time /T`;
 233         if ( -e "output" ) {
 234           cp( "output", $outputfile );
 235           unlink( "output" );
 236         } elsif ( -e "OUTPUT" ) {
 237           cp( "OUTPUT", $outputfile );
 238           unlink( "OUTPUT" );
 239         }
 240       } else {
 241         unless( $self -> {'show_version'} ){
 242           $version = '';
 243         }
 244         $start_time = `date`;
 245         run3( "nice -n " . $self -> {'nice'} . " ./nonmem$version$adaptive", "FCON", $outputfile );
 246         $fin_time = `date`;
 247       }
 248       open( OUT, '>>', $outputfile );
 249       print( OUT  "This file was created using the NONMEM version in directory $nmdir\n" );
 250       print( OUT "Started  $start_time" );
 251       print( OUT "Finished $fin_time" );
 252       close( OUT );
 253     }
 254 end execute
 255
 256 # }}}