#ifndef DYNAMITEgenewise4HEADERFILE #define DYNAMITEgenewise4HEADERFILE #ifdef _cplusplus extern "C" { #endif #include "dyna.h" #include "geneparser4.h" #include "genewisemodel.h" #include "genewisemodeldb.h" struct Wise2_GeneWise4 { int dynamite_hard_link; BaseMatrix * basematrix; int leni; int lenj; GeneWiseScore* query; ComplexSequence* target; GeneParser4Score * gp; } ; /* GeneWise4 defined */ #ifndef DYNAMITE_DEFINED_GeneWise4 typedef struct Wise2_GeneWise4 Wise2_GeneWise4; #define GeneWise4 Wise2_GeneWise4 #define DYNAMITE_DEFINED_GeneWise4 #endif #ifdef PTHREAD struct thread_pool_holder_GeneWise4 { GeneWiseScore* query; /* Query object placeholder */ GeneWiseDB* querydb;/* Query database object */ boolean query_init; ComplexSequence* target;/* Target object placeholder */ GenomicDB* targetdb;/* Target database object */ boolean target_init; GeneParser4Score * gp; pthread_mutex_t input_lock; pthread_mutex_t output_lock; Hscore * out; pthread_t * pool; int number_of_threads; boolean search_has_ended; DBSearchImpl * dbsi; } ; #endif /* PTHREAD */ /***************************************************/ /* Callable functions */ /* These are the functions you are expected to use */ /***************************************************/ /* Function: search_GeneWise4(dbsi,out,querydb,targetdb,gp) * * Descrip: This function makes a database search of GeneWise4 * It uses the dbsi structure to choose which implementation * to use of the database searching. This way at run time you * can switch between single threaded/multi-threaded or hardware * * * Arg: dbsi [UNKN ] Undocumented argument [DBSearchImpl *] * Arg: out [UNKN ] Undocumented argument [Hscore *] * Arg: querydb [UNKN ] Undocumented argument [GeneWiseDB*] * Arg: targetdb [UNKN ] Undocumented argument [GenomicDB*] * Arg: gp [UNKN ] Undocumented argument [GeneParser4Score *] * * Return [UNKN ] Undocumented return value [Search_Return_Type] * */ Search_Return_Type Wise2_search_GeneWise4(DBSearchImpl * dbsi,Hscore * out,GeneWiseDB* querydb,GenomicDB* targetdb ,GeneParser4Score * gp); #define search_GeneWise4 Wise2_search_GeneWise4 /* Function: serial_search_GeneWise4(out,querydb,targetdb,gp) * * Descrip: This function makes a database search of GeneWise4 * It is a single processor implementation * * * Arg: out [UNKN ] Undocumented argument [Hscore *] * Arg: querydb [UNKN ] Undocumented argument [GeneWiseDB*] * Arg: targetdb [UNKN ] Undocumented argument [GenomicDB*] * Arg: gp [UNKN ] Undocumented argument [GeneParser4Score *] * * Return [UNKN ] Undocumented return value [Search_Return_Type] * */ Search_Return_Type Wise2_serial_search_GeneWise4(Hscore * out,GeneWiseDB* querydb,GenomicDB* targetdb ,GeneParser4Score * gp); #define serial_search_GeneWise4 Wise2_serial_search_GeneWise4 /* Function: PackAln_bestmemory_GeneWise4(query,target,gp,dpenv,dpri) * * Descrip: This function chooses the best memory set-up for the alignment * using calls to basematrix, and then implements either a large * or small memory model. * * It is the best function to use if you just want an alignment * * If you want a label alignment, you will need * /convert_PackAln_to_AlnBlock_GeneWise4 * * * Arg: query [UNKN ] query data structure [GeneWiseScore*] * Arg: target [UNKN ] target data structure [ComplexSequence*] * Arg: gp [UNKN ] Resource [GeneParser4Score *] * Arg: dpenv [UNKN ] Undocumented argument [DPEnvelope *] * Arg: dpri [UNKN ] Undocumented argument [DPRunImpl *] * * Return [UNKN ] Undocumented return value [PackAln *] * */ PackAln * Wise2_PackAln_bestmemory_GeneWise4(GeneWiseScore* query,ComplexSequence* target ,GeneParser4Score * gp,DPEnvelope * dpenv,DPRunImpl * dpri); #define PackAln_bestmemory_GeneWise4 Wise2_PackAln_bestmemory_GeneWise4 /* Function: allocate_Expl_GeneWise4(query,target,gp) * * Descrip: This function allocates the GeneWise4 structure * and the basematrix area for explicit memory implementations * It calls /allocate_GeneWise4_only * * * Arg: query [UNKN ] query data structure [GeneWiseScore*] * Arg: target [UNKN ] target data structure [ComplexSequence*] * Arg: gp [UNKN ] Resource [GeneParser4Score *] * * Return [UNKN ] Undocumented return value [GeneWise4 *] * */ GeneWise4 * Wise2_allocate_Expl_GeneWise4(GeneWiseScore* query,ComplexSequence* target ,GeneParser4Score * gp); #define allocate_Expl_GeneWise4 Wise2_allocate_Expl_GeneWise4 /* Function: recalculate_PackAln_GeneWise4(pal,mat) * * Descrip: This function recalculates the PackAln structure produced by GeneWise4 * For example, in linear space methods this is used to score them * * * Arg: pal [UNKN ] Undocumented argument [PackAln *] * Arg: mat [UNKN ] Undocumented argument [GeneWise4 *] * */ void Wise2_recalculate_PackAln_GeneWise4(PackAln * pal,GeneWise4 * mat); #define recalculate_PackAln_GeneWise4 Wise2_recalculate_PackAln_GeneWise4 /* Function: allocate_Small_GeneWise4(query,target,gp) * * Descrip: This function allocates the GeneWise4 structure * and the basematrix area for a small memory implementations * It calls /allocate_GeneWise4_only * * * Arg: query [UNKN ] query data structure [GeneWiseScore*] * Arg: target [UNKN ] target data structure [ComplexSequence*] * Arg: gp [UNKN ] Resource [GeneParser4Score *] * * Return [UNKN ] Undocumented return value [GeneWise4 *] * */ GeneWise4 * Wise2_allocate_Small_GeneWise4(GeneWiseScore* query,ComplexSequence* target ,GeneParser4Score * gp); #define allocate_Small_GeneWise4 Wise2_allocate_Small_GeneWise4 /* Function: PackAln_calculate_Small_GeneWise4(mat,dpenv) * * Descrip: This function calculates an alignment for GeneWise4 structure in linear space * If you want only the start/end points * use /AlnRangeSet_calculate_Small_GeneWise4 * * The function basically * finds start/end points * foreach start/end point * calls /full_dc_GeneWise4 * * * Arg: mat [UNKN ] Undocumented argument [GeneWise4 *] * Arg: dpenv [UNKN ] Undocumented argument [DPEnvelope *] * * Return [UNKN ] Undocumented return value [PackAln *] * */ PackAln * Wise2_PackAln_calculate_Small_GeneWise4(GeneWise4 * mat,DPEnvelope * dpenv); #define PackAln_calculate_Small_GeneWise4 Wise2_PackAln_calculate_Small_GeneWise4 /* Function: AlnRangeSet_calculate_Small_GeneWise4(mat) * * Descrip: This function calculates an alignment for GeneWise4 structure in linear space * If you want the full alignment, use /PackAln_calculate_Small_GeneWise4 * If you have already got the full alignment, but want the range set, use /AlnRangeSet_from_PackAln_GeneWise4 * If you have got the small matrix but not the alignment, use /AlnRangeSet_from_GeneWise4 * * * Arg: mat [UNKN ] Undocumented argument [GeneWise4 *] * * Return [UNKN ] Undocumented return value [AlnRangeSet *] * */ AlnRangeSet * Wise2_AlnRangeSet_calculate_Small_GeneWise4(GeneWise4 * mat); #define AlnRangeSet_calculate_Small_GeneWise4 Wise2_AlnRangeSet_calculate_Small_GeneWise4 /* Function: AlnRangeSet_from_GeneWise4(mat) * * Descrip: This function reads off a start/end structure * for GeneWise4 structure in linear space * If you want the full alignment use * /PackAln_calculate_Small_GeneWise4 * If you have not calculated the matrix use * /AlnRange_calculate_Small_GeneWise4 * * * Arg: mat [UNKN ] Undocumented argument [GeneWise4 *] * * Return [UNKN ] Undocumented return value [AlnRangeSet *] * */ AlnRangeSet * Wise2_AlnRangeSet_from_GeneWise4(GeneWise4 * mat); #define AlnRangeSet_from_GeneWise4 Wise2_AlnRangeSet_from_GeneWise4 /* Function: convert_PackAln_to_AlnBlock_GeneWise4(pal) * * Descrip: Converts a path alignment to a label alignment * The label alignment is probably much more useful than the path * * * Arg: pal [UNKN ] Undocumented argument [PackAln *] * * Return [UNKN ] Undocumented return value [AlnBlock *] * */ AlnBlock * Wise2_convert_PackAln_to_AlnBlock_GeneWise4(PackAln * pal); #define convert_PackAln_to_AlnBlock_GeneWise4 Wise2_convert_PackAln_to_AlnBlock_GeneWise4 /* Function: PackAln_read_Expl_GeneWise4(mat) * * Descrip: Reads off PackAln from explicit matrix structure * * * Arg: mat [UNKN ] Undocumented argument [GeneWise4 *] * * Return [UNKN ] Undocumented return value [PackAln *] * */ PackAln * Wise2_PackAln_read_Expl_GeneWise4(GeneWise4 * mat); #define PackAln_read_Expl_GeneWise4 Wise2_PackAln_read_Expl_GeneWise4 /* Function: calculate_GeneWise4(mat) * * Descrip: This function calculates the GeneWise4 matrix when in explicit mode * To allocate the matrix use /allocate_Expl_GeneWise4 * * * Arg: mat [UNKN ] GeneWise4 which contains explicit basematrix memory [GeneWise4 *] * * Return [UNKN ] Undocumented return value [boolean] * */ boolean Wise2_calculate_GeneWise4(GeneWise4 * mat); #define calculate_GeneWise4 Wise2_calculate_GeneWise4 /* Function: GeneWise4_alloc(void) * * Descrip: Allocates structure: assigns defaults if given * * * * Return [UNKN ] Undocumented return value [GeneWise4 *] * */ GeneWise4 * Wise2_GeneWise4_alloc(void); #define GeneWise4_alloc Wise2_GeneWise4_alloc /* Function: free_GeneWise4(obj) * * Descrip: Free Function: removes the memory held by obj * Will chain up to owned members and clear all lists * * * Arg: obj [UNKN ] Object that is free'd [GeneWise4 *] * * Return [UNKN ] Undocumented return value [GeneWise4 *] * */ GeneWise4 * Wise2_free_GeneWise4(GeneWise4 * obj); #define free_GeneWise4 Wise2_free_GeneWise4 /* Unplaced functions */ /* There has been no indication of the use of these functions */ /***************************************************/ /* Internal functions */ /* you are not expected to have to call these */ /***************************************************/ void * Wise2_thread_loop_GeneWise4(void * ptr); #define thread_loop_GeneWise4 Wise2_thread_loop_GeneWise4 int Wise2_score_only_GeneWise4(GeneWiseScore* query,ComplexSequence* target ,GeneParser4Score * gp); #define score_only_GeneWise4 Wise2_score_only_GeneWise4 GeneWise4 * Wise2_allocate_GeneWise4_only(GeneWiseScore* query,ComplexSequence* target ,GeneParser4Score * gp); #define allocate_GeneWise4_only Wise2_allocate_GeneWise4_only void Wise2_init_GeneWise4(GeneWise4 * mat); #define init_GeneWise4 Wise2_init_GeneWise4 AlnRange * Wise2_AlnRange_build_GeneWise4(GeneWise4 * mat,int stopj,int stopspecstate,int * startj,int * startspecstate); #define AlnRange_build_GeneWise4 Wise2_AlnRange_build_GeneWise4 boolean Wise2_read_hidden_GeneWise4(GeneWise4 * mat,int starti,int startj,int startstate,int stopi,int stopj,int stopstate,PackAln * out); #define read_hidden_GeneWise4 Wise2_read_hidden_GeneWise4 int Wise2_max_hidden_GeneWise4(GeneWise4 * mat,int hiddenj,int i,int j,int state,boolean isspecial,int * reti,int * retj,int * retstate,boolean * retspecial,int * cellscore); #define max_hidden_GeneWise4 Wise2_max_hidden_GeneWise4 boolean Wise2_read_special_strip_GeneWise4(GeneWise4 * mat,int stopi,int stopj,int stopstate,int * startj,int * startstate,PackAln * out); #define read_special_strip_GeneWise4 Wise2_read_special_strip_GeneWise4 int Wise2_max_special_strip_GeneWise4(GeneWise4 * mat,int i,int j,int state,boolean isspecial,int * reti,int * retj,int * retstate,boolean * retspecial,int * cellscore); #define max_special_strip_GeneWise4 Wise2_max_special_strip_GeneWise4 int Wise2_max_matrix_to_special_GeneWise4(GeneWise4 * mat,int i,int j,int state,int cscore,int * reti,int * retj,int * retstate,boolean * retspecial,int * cellscore); #define max_matrix_to_special_GeneWise4 Wise2_max_matrix_to_special_GeneWise4 void Wise2_calculate_hidden_GeneWise4(GeneWise4 * mat,int starti,int startj,int startstate,int stopi,int stopj,int stopstate,DPEnvelope * dpenv); #define calculate_hidden_GeneWise4 Wise2_calculate_hidden_GeneWise4 void Wise2_init_hidden_GeneWise4(GeneWise4 * mat,int starti,int startj,int stopi,int stopj); #define init_hidden_GeneWise4 Wise2_init_hidden_GeneWise4 boolean Wise2_full_dc_GeneWise4(GeneWise4 * mat,int starti,int startj,int startstate,int stopi,int stopj,int stopstate,PackAln * out,int * donej,int totalj,DPEnvelope * dpenv); #define full_dc_GeneWise4 Wise2_full_dc_GeneWise4 boolean Wise2_do_dc_single_pass_GeneWise4(GeneWise4 * mat,int starti,int startj,int startstate,int stopi,int stopj,int stopstate,DPEnvelope * dpenv,int perc_done); #define do_dc_single_pass_GeneWise4 Wise2_do_dc_single_pass_GeneWise4 void Wise2_push_dc_at_merge_GeneWise4(GeneWise4 * mat,int starti,int stopi,int startj,int * stopj,DPEnvelope * dpenv); #define push_dc_at_merge_GeneWise4 Wise2_push_dc_at_merge_GeneWise4 void Wise2_follow_on_dc_GeneWise4(GeneWise4 * mat,int starti,int stopi,int startj,int stopj,DPEnvelope * dpenv,int perc_done); #define follow_on_dc_GeneWise4 Wise2_follow_on_dc_GeneWise4 void Wise2_run_up_dc_GeneWise4(GeneWise4 * mat,int starti,int stopi,int startj,int stopj,DPEnvelope * dpenv,int perc_done); #define run_up_dc_GeneWise4 Wise2_run_up_dc_GeneWise4 void Wise2_init_dc_GeneWise4(GeneWise4 * mat); #define init_dc_GeneWise4 Wise2_init_dc_GeneWise4 boolean Wise2_dc_start_end_calculate_GeneWise4(GeneWise4 * mat,DPEnvelope * dpenv); #define dc_start_end_calculate_GeneWise4 Wise2_dc_start_end_calculate_GeneWise4 int Wise2_start_end_find_end_GeneWise4(GeneWise4 * mat,int * endj); #define start_end_find_end_GeneWise4 Wise2_start_end_find_end_GeneWise4 void Wise2_init_start_end_linear_GeneWise4(GeneWise4 * mat); #define init_start_end_linear_GeneWise4 Wise2_init_start_end_linear_GeneWise4 AlnConvertSet * Wise2_AlnConvertSet_GeneWise4(void); #define AlnConvertSet_GeneWise4 Wise2_AlnConvertSet_GeneWise4 int Wise2_find_end_GeneWise4(GeneWise4 * mat,int * ri,int * rj,int * state,boolean * isspecial); #define find_end_GeneWise4 Wise2_find_end_GeneWise4 void Wise2_GeneWise4_debug_show_matrix(GeneWise4 * mat,int starti,int stopi,int startj,int stopj,FILE * ofp); #define GeneWise4_debug_show_matrix Wise2_GeneWise4_debug_show_matrix int Wise2_max_calc_GeneWise4(GeneWise4 * mat,int i,int j,int state,boolean isspecial,int * reti,int * retj,int * retstate,boolean * retspecial,int * cellscore); #define max_calc_GeneWise4 Wise2_max_calc_GeneWise4 int Wise2_max_calc_special_GeneWise4(GeneWise4 * mat,int i,int j,int state,boolean isspecial,int * reti,int * retj,int * retstate,boolean * retspecial,int * cellscore); #define max_calc_special_GeneWise4 Wise2_max_calc_special_GeneWise4 #ifdef _cplusplus } #endif #endif