File: promotercluster.c

package info (click to toggle)
wise 2.4.1-21
  • links: PTS, VCS
  • area: main
  • in suites: bullseye, buster, sid
  • size: 27,140 kB
  • sloc: ansic: 276,365; makefile: 1,003; perl: 886; lex: 93; yacc: 81; sh: 24
file content (68 lines) | stat: -rw-r--r-- 1,646 bytes parent folder | download
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
#include "dnaprofileengine.h"
#include "version.h"


char * program_name = "promotercluster";

void show_version(FILE * ofp)
{
  fprintf(ofp,"%s\nVersion: %s\nReleased: %s\nCompiled: %s\n",program_name,VERSION_NUMBER,RELEASE_DAY,COMPILE_DATE);
  fprintf(ofp,"\nThis program is freely distributed under a GNU General Public License\n");
  fprintf(ofp,"The source code is copyright (c) EMBL and others\n");
  fprintf(ofp,"There is no warranty, implied or otherwise on the performance of this program\n");
  fprintf(ofp,"For more information read the GNULICENSE file in the distribution\n\n");
  fprintf(ofp,"Credits: Ewan Birney <birney@ebi.ac.uk>\n");
  exit(63);   
}

void show_help(FILE * ofp)
{
  fprintf(ofp,"%s fasta-file-for-promoter-set\n",program_name);

  show_help_DnaProfileEnginePara(ofp);

  show_standard_options(ofp);

}


int main(int argc,char ** argv)
{
  DnaProfileEnginePara * dpep;
  DnaProfileNode * root;
  DnaProfileSet * set;
  FILE * ifp;

  boolean is_four = FALSE;

  dpep = new_DnaProfileEnginePara_from_argv(&argc,argv);

  is_four = strip_out_boolean_argument(&argc,argv,"bfour");

  strip_out_standard_options(&argc,argv,show_help,show_version);


  if( argc != 2 ) {
    show_help(stdout);
    exit(12);
  }

  ifp = openfile(argv[1],"r");
  if( ifp == NULL ) {
    fatal("could not open file %s",argv[1]);
  }

  if( is_four ) {
    root = balanced_4_Sequence_fasta_stream(ifp);
  } else {
    root = simple_cascade_Sequence_fasta_stream(ifp);
  }

  populate_DnaProfileNode_from_root(root,dpep);

  /*set = filter_DnaProfileSet(root->set,0,0);*/

  show_DnaProfileSet(root->set,dpep->rm,stdout);

  return 0;
}