X-Git-Url: http://git.kpe.io/?a=blobdiff_plain;f=cgi-bin%2Fctsim.cgi.in;h=660ac60b504445ad9182d9db2427424b2b6030c4;hb=ba39c5a65a2aeeb695b1434b1d7cb781d110c1b2;hp=c4d9e531b39dadb4f4a560c7b5d1119af9b97bda;hpb=13838bda88b63a2535b5baaf7197006767de4b8e;p=ctsim.git diff --git a/cgi-bin/ctsim.cgi.in b/cgi-bin/ctsim.cgi.in index c4d9e53..660ac60 100755 --- a/cgi-bin/ctsim.cgi.in +++ b/cgi-bin/ctsim.cgi.in @@ -6,9 +6,10 @@ use CGI; use File::Basename; use Fcntl ':flock'; +require './ctsim.conf'; + my $opt_d = 0; -my $bindir = "@prefix@/bin"; -$ENV{'PATH'} = $bindir; +$ENV{'PATH'} = "$::bindir:$::lamrundir"; my $fromhost = $ENV{'REMOTE_HOST'}; @@ -55,59 +56,56 @@ my $Disp_Min = "auto"; my $Disp_Max = "auto"; $Disp_Min = FilterToNumber($in{'Disp_Min'}) if ($in{'Disp_Min'} ne "auto" && $in{'Disp_Min'} ne ""); $Disp_Max = FilterToNumber($in{'Disp_Max'}) if ($in{'Disp_Max'} ne "auto" && $in{'Disp_Max'} ne ""); -if ($Disp_Min ne 'auto' && ! ($Disp_Min =~ /^[\d\.]+$/)) { - $error .= "Display Minimum must be 'auto' or numeric (received '$Disp_Min')
"; +if ($Disp_Min ne 'auto' && ! ($Disp_Min =~ /^[\d\.\-]+$/)) { + $error .= "Display Minimum must be 'auto' or numeric (received $Disp_Min)
"; } -if ($Disp_Max ne 'auto' && ! ($Disp_Max =~ /^[\d\.]+$/)) { - $error .= "Display Maximum must be 'auto' or numeric (received '$Disp_Max')
"; +if ($Disp_Max ne 'auto' && ! ($Disp_Max =~ /^[\d\.\-]+$/)) { + $error .= "Display Maximum must be 'auto' or numeric (received $Disp_Max)
"; } my $MPI_Str = FilterMetaChars($in{'MPI'}); my $MPI = 0; -$MPI = 1 if ($MPI_Str eq "yes"); +$MPI = 1 if ($MPI_Str eq "yes" && $::mpi_enable ne ""); $error .= "IR Nx and Ny must be between 5 and 1024
" if ($IR_Nx < 5 || $IR_Nx > 1024 || $IR_Ny < 5 || $IR_Ny > 1024); $error .= "IR Filter Parameter must be between 0 and 1
" if ($IR_Filter_Param < 0 || $IR_Filter_Param > 1); -my $jobdir = "@webdatadir@"; my $tmpid = $$; my $auto_window_img = "std0.1"; my $auto_window_diff = "std1"; my $auto_window_rs = "full"; -my $datadir = "@webdatadir@"; -my $url_datadir = "@webdataurl@"; -my $logfile = "$jobdir/ctsim.log"; - -my $result_fname = "$datadir/result-$tmpid.html"; -my $phantom_fname = "$datadir/phantom-$tmpid.sdf"; -my $rs_fname = "$datadir/rs-$tmpid.rs"; -my $ir_fname = "$datadir/ir-$tmpid.sdf"; -my $rs_sdf_fname = "$datadir/rs-$tmpid.sdf"; -my $diff_fname = "$datadir/diff-$tmpid.sdf"; -my $phantom_png = "$datadir/phantom-$tmpid.png"; -my $ir_png = "$datadir/ir-$tmpid.png"; -my $rs_png = "$datadir/rs-$tmpid.png"; -my $diff_png = "$datadir/diff-$tmpid.png"; - -my $result_url = "$url_datadir/result-$tmpid.html"; -my $phantom_png_url = "$url_datadir/phantom-$tmpid.png"; -my $ir_png_url = "$url_datadir/ir-$tmpid.png"; -my $rs_png_url = "$url_datadir/rs-$tmpid.png"; -my $diff_png_url = "$url_datadir/diff-$tmpid.png"; - -my $ir_ver = "$bindir/ir"; -my $rs_ver = "$bindir/rs"; -my $gp_ver = "$bindir/gp"; -my $diff_ver = "$bindir/sdf-2"; -$ir_ver = "/opt/lam/bin/mpirun N N $bindir/ir-mpi" if $MPI; -$rs_ver = "/opt/lam/bin/mpirun N N $bindir/rs-mpi" if $MPI; -$gp_ver = "/opt/lam/bin/mpirun N N $bindir/gp-mpi" if $MPI; - -my $gp_cmd = "$gp_ver $phantom_fname $Phantom_Nx $Phantom_Ny --phantom $Phantom_Name --nsample $Phantom_NSample"; -my $rs_cmd = "$rs_ver $rs_fname $RS_NDet $RS_NRot --phantom $Phantom_Name --nray $RS_NRay --rotangle $RS_RotAngle"; -my $rs_sdf_cmd = "$bindir/rs2sdf $rs_fname $rs_sdf_fname"; -my $ir_cmd = "$ir_ver $rs_fname $ir_fname $IR_Nx $IR_Ny --filter $IR_Filter --filter-param $IR_Filter_Param --interp $IR_Interp --backproj $IR_Backproj"; +my $logfile = "$::jobdir/ctsim.log"; + +my $result_fname = "$::datadir/result-$tmpid.html"; +my $phantom_fname = "$::datadir/phantom-$tmpid.sdf"; +my $rs_fname = "$::datadir/rs-$tmpid.rs"; +my $ir_fname = "$::datadir/ir-$tmpid.sdf"; +my $rs_sdf_fname = "$::datadir/rs-$tmpid.sdf"; +my $diff_fname = "$::datadir/diff-$tmpid.sdf"; +my $phantom_png = "$::datadir/phantom-$tmpid.png"; +my $ir_png = "$::datadir/ir-$tmpid.png"; +my $rs_png = "$::datadir/rs-$tmpid.png"; +my $diff_png = "$::datadir/diff-$tmpid.png"; + +my $result_url = "$::url_datadir/result-$tmpid.html"; +my $phantom_png_url = "$::url_datadir/phantom-$tmpid.png"; +my $ir_png_url = "$::url_datadir/ir-$tmpid.png"; +my $rs_png_url = "$::url_datadir/rs-$tmpid.png"; +my $diff_png_url = "$::url_datadir/diff-$tmpid.png"; + +my $ctrec_ver = "$::bindir/ctrec"; +my $phm2rs_ver = "$::bindir/phm2rs"; +my $phm2sdf_ver = "$::bindir/phm2sdf"; +my $diff_ver = "$::bindir/sdf-2"; +$ctrec_ver = "/opt/lam/bin/mpirun N N $::lamrundir/ctrec-lam" if $MPI; +$phm2rs_ver = "/opt/lam/bin/mpirun N N $::lamrundir/phm2rs-lam" if $MPI; +$phm2sdf_ver = "/opt/lam/bin/mpirun N N $::lamrundir/phm2sdf-lam" if $MPI; + +my $gp_cmd = "$phm2sdf_ver $phantom_fname $Phantom_Nx $Phantom_Ny --phantom $Phantom_Name --nsample $Phantom_NSample"; +my $rs_cmd = "$phm2rs_ver $rs_fname $RS_NDet $RS_NRot --phantom $Phantom_Name --nray $RS_NRay --rotangle $RS_RotAngle"; +my $rs_sdf_cmd = "$::bindir/rs2sdf $rs_fname $rs_sdf_fname"; +my $ir_cmd = "$ctrec_ver $rs_fname $ir_fname $IR_Nx $IR_Ny --filter $IR_Filter --filter-param $IR_Filter_Param --interp $IR_Interp --backproj $IR_Backproj"; my $diff_cmd = "$diff_ver $phantom_fname $ir_fname $diff_fname --comp"; my $window_options = "--auto $auto_window_img"; @@ -118,10 +116,10 @@ if ($Disp_Max ne 'auto') { $window_options .= " --max $Disp_Max"; } -my $png1_cmd = "$bindir/sdf2img $phantom_fname $phantom_png $window_options --stats --format png"; -my $png2_cmd = "$bindir/sdf2img $ir_fname $ir_png $window_options --stats --format png"; -my $png3_cmd = "$bindir/sdf2img $rs_sdf_fname $rs_png --auto $auto_window_rs --stats --format png"; -my $png4_cmd = "$bindir/sdf2img $diff_fname $diff_png --auto $auto_window_diff --stats --format png"; +my $png1_cmd = "$::bindir/sdf2img $phantom_fname $phantom_png $window_options --stats --format png"; +my $png2_cmd = "$::bindir/sdf2img $ir_fname $ir_png $window_options --stats --format png"; +my $png3_cmd = "$::bindir/sdf2img $rs_sdf_fname $rs_png --auto $auto_window_rs --stats --format png"; +my $png4_cmd = "$::bindir/sdf2img $diff_fname $diff_png --auto $auto_window_diff --stats --format png"; my $title = "CT Simulation Results"; @@ -142,68 +140,62 @@ if ($error ne "") { $out .= "Please correct them and try submitting your request again.
\n"; $out .= "$error"; } else { - my $gp_out; - my $rs_out; - my $rs_sdf_out; - my $ir_out; - my $diff_out; - my $png_gp_out; - my $png_ir_out; - my $png_rs_out; - my $png_diff_out; - if (1) { - $gp_out = `$gp_cmd`; - $rs_out .= `$rs_cmd`; - $rs_sdf_out .= `$rs_sdf_cmd`; - $ir_out .= `$ir_cmd`; - $diff_out .= `$diff_cmd`; - $png_gp_out .= `$png1_cmd`; + my $gp_out; + my $rs_out; + my $rs_sdf_out; + my $ir_out; + my $diff_out; + my $png_gp_out; + my $png_ir_out; + my $png_rs_out; + my $png_diff_out; + $gp_out = `$gp_cmd`; + if (-s $phantom_fname) { + $rs_out .= `$rs_cmd`; + $png_gp_out .= `$png1_cmd`; + if (-s $rs_fname) { + $rs_sdf_out .= `$rs_sdf_cmd`; + $png_rs_out .= `$png3_cmd`; + $ir_out .= `$ir_cmd`; + if (-s $ir_fname) { $png_ir_out .= `$png2_cmd`; - $png_rs_out .= `$png3_cmd`; + $diff_out .= `$diff_cmd`; $png_diff_out .= `$png4_cmd`; - } else { - $gp_out = `$gp_cmd 2>&1`; - $rs_out .= `$rs_cmd 2>&1`; - $rs_sdf_out .= `$rs_sdf_cmd 2>&1`; - $ir_out .= `$ir_cmd 2>&1`; - $diff_out .= `$diff_cmd 2>&1`; - $png_gp_out .= `$png1_cmd 2>&1`; - $png_ir_out .= `$png2_cmd 2>&1`; - $png_rs_out .= `$png3_cmd 2>&1`; - $png_diff_out .= `$png4_cmd 2>&1`; - } - $cmdout = "$gp_cmd\n $gp_out $rs_cmd\n $rs_out $rs_sdf_cmd\n $rs_sdf_out $ir_cmd\n $ir_out $diff_cmd\n $diff_out $png1_cmd\n $png_gp_out $png2_cmd\n $png_ir_out $png3_cmd\n $png_rs_out $png4_cmd\n $png_diff_out"; - if (open(LOGFILE,">> $logfile")) { - flock(LOGFILE,LOCK_EX); - seek(LOGFILE, 0, 2); - print LOGFILE "Job $tmpid\n"; - print LOGFILE $cmdout; - print LOGFILE "----------------------------------------------------\n"; - flock(LOGFILE,LOCK_UN); - close(LOGFILE); - } - $gmt_end = time(); - $execution_time = $gmt_end - $gmt_start; - if ($opt_d) - { - $out .= "

Command Output

$cmdout
\n"; + } } - my $png_gp_out_html = $png_gp_out; - my $png_ir_out_html = $png_ir_out; - my $png_rs_out_html = $png_rs_out; - my $png_diff_out_html = $png_diff_out; - $png_gp_out_html =~ s/\n/
/gms; - $png_ir_out_html =~ s/\n/
/gms; - $png_rs_out_html =~ s/\n/
/gms; - $png_diff_out_html =~ s/\n/
/gms; - $out .= "\n"; - $out .= "\n"; - $out .= "\n"; - $out .= "\n"; - $out .= "\n"; - $out .= "\n"; - $out .= "
Phantom ImageReconstructed Image

$png_gp_out
$png_ir_out
Raysum SinusoidPhantom/Reconst Error

$png_rs_out
$diff_out
$png_diff_out
"; - $out .= "Execution time: $execution_time seconds\n"; + } + + $cmdout = "$gp_cmd\n $gp_out $rs_cmd\n $rs_out $rs_sdf_cmd\n $rs_sdf_out $ir_cmd\n $ir_out $diff_cmd\n $diff_out $png1_cmd\n $png_gp_out $png2_cmd\n $png_ir_out $png3_cmd\n $png_rs_out $png4_cmd\n $png_diff_out"; + if (open(LOGFILE,">> $logfile")) { + flock(LOGFILE,LOCK_EX); + seek(LOGFILE, 0, 2); + print LOGFILE "Job $tmpid\n"; + print LOGFILE $cmdout; + print LOGFILE "----------------------------------------------------\n"; + flock(LOGFILE,LOCK_UN); + close(LOGFILE); + } + $gmt_end = time(); + $execution_time = $gmt_end - $gmt_start; + if ($opt_d) { + $out .= "

Command Output

$cmdout
\n"; + } + my $png_gp_out_html = $png_gp_out; + my $png_ir_out_html = $png_ir_out; + my $png_rs_out_html = $png_rs_out; + my $png_diff_out_html = $png_diff_out; + $png_gp_out_html =~ s/\n/
/gms; + $png_ir_out_html =~ s/\n/
/gms; + $png_rs_out_html =~ s/\n/
/gms; + $png_diff_out_html =~ s/\n/
/gms; + $out .= "\n"; + $out .= "\n"; + $out .= "\n"; + $out .= "\n"; + $out .= "\n"; + $out .= "\n"; + $out .= "
Phantom ImageReconstructed Image

$png_gp_out
$png_ir_out
Raysum SinusoidPhantom/Reconst Error

$png_rs_out
$diff_out
$png_diff_out
"; + $out .= "Execution time: $execution_time seconds\n"; } $out .= "
\n"; @@ -225,7 +217,7 @@ else print "Content-type: text/plain\n\n"; print "The simulator was unable to create an result file.\n"; } -if (open(JOBFILES,"> $jobdir/$tmpid")) +if (open(JOBFILES,"> $::jobdir/$tmpid")) { flock(JOBFILES,LOCK_EX); print JOBFILES "gmt_start=$gmt_start\n"; @@ -284,6 +276,6 @@ sub FilterMetaChars sub FilterToNumber { my $var = pop(@_); - $var =~ /^(\d*\.*\d*)$/; + $var =~ /^(\-*\d*\.*\d*)$/; $1; }