Exemplo n.º 1
0
$wsize = mysql_real_escape_string($_POST['wsize']);
$minqHSP = mysql_real_escape_string($_POST['minqHSP']);
$minqHSPc = mysql_real_escape_string($_POST['minqHSPc']);
$minESTc = mysql_real_escape_string($_POST['minESTc']);
$maxnest = mysql_real_escape_string($_POST['maxnest']);
$gsq_job_name = mysql_real_escape_string($_POST['gsq_job_name']);
$job_name_submit = job_name_strip($gsq_job_name);
// jobs_functions.inc.php; strip any user-added illegal chars
// calculate file sizes
$trans_array_est = create_input_list($input_data_path, "est", $dbpass);
# calculate total file size - est (jobs_functions.inc.php)
$trans_array_cdna = create_input_list($input_data_path, "cdna", $dbpass);
# calculate total file size - cdna (jobs_functions.inc.php)
$trans_array_tsa = create_input_list($input_data_path, "tsa", $dbpass);
# calculate total file size - tsa (jobs_functions.inc.php)
$gdna_array = create_input_list($input_data_path, "gdna", $dbpass);
# calculate total file size - genome
$genome_file_size = $gdna_array[2];
// for jobs database as total genome file size
// Get processor info from database based on selected app_id (note: 'nodes' is not required, the app is pre-configured for node count, but we sent it anyway for informational purposes)
$query = "SELECT nodes, proc_per_node, memory_per_node, platform from {$global_DB}.apps where app_id ='{$app_id}'";
$result = mysql_query($query);
$row = mysql_fetch_array($result);
$nodes = $row['nodes'];
$HPC_name = $row['platform'];
$proc_per_node = $row['proc_per_node'];
// should be max
$memory_per_node = $row['memory_per_node'];
// default is '2' currently
$transcript_type = mysql_real_escape_string($_POST['transcript_type']);
//used to assign est, cdna, tsa. See below:
Exemplo n.º 2
0
 # this is the default
 $editable = "";
 ## GenomeThreader parameters
 $input_g = "{$user_input_base}SCFDIR/{$DBid}gdna.fa";
 # temp input path for Genome File e.g. /username/hpcGDB002/SCFDIR/GDB002gdna.fa
 $input_p = "{$user_input_base}Protein/{$DBid}prot.fa";
 # temp input path for Transcript File e.g. /username/hpcGDB002/Protein/GDB002prot.fa
 ##  Set validation style for input files
 $valid_est_array = create_input_list($gsq_input_path, "transcript", $dbpass);
 // generic; covers est, cdna, tsa
 $valid_est = $valid_est_array[0];
 $valid_prot_array = create_input_list($gsq_input_path, "protein", $dbpass);
 $valid_prot = $valid_prot_array[0];
 $valid_gdna_array = create_input_list($gsq_input_path, "gdna", $dbpass);
 $valid_gdna = $valid_gdna_array[0];
 $valid_out_array = create_input_list($output_data_path, "gsq." . $DBid, $dbpass);
 # e.g. 	gsq.GDB001cdna, test for type in function.
 $valid_out = $valid_out_array[0];
 ## Set return job estimates based on size distribution of gdna:
 $scaffolds = calculate_scaffolds($gsq_input_path, $gsq_proc, $gth_proc);
 # jobs_functions.inc.php
 // return array($large_scaffold_count, $small_scaffold_count, $chunks, $remainder_size, $gsq_split, $gth_split, $scaffold_size_display, $gsq_split_display, $gth_split_display, $gsq_time_display, $gth_time_display);
 $large_scaffold_count = $scaffolds[0];
 $small_scaffold_count = $scaffolds[1];
 $total_scaffold_count = $large_scaffold_count + $small_scaffold_count;
 $break_point_count = $scaffolds[2];
 $remainder_size = $scaffolds[3];
 $gsq_split = $scaffolds[4];
 $gth_split = $scaffolds[5];
 $scaffold_sizes_display = $scaffolds[6];
 $gsq_split_display = $scaffolds[7];