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
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
|
/* @source martfilters application
**
** Return biomart filters given mart dataset info
**
** @author Copyright (C) Alan Bleasby
** @@
**
** This program is free software; you can redistribute it and/or
** modify it under the terms of the GNU General Public License
** as published by the Free Software Foundation; either version 2
** of the License, or (at your option) any later version.
**
** This program is distributed in the hope that it will be useful,
** but WITHOUT ANY WARRANTY; without even the implied warranty of
** MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
** GNU General Public License for more details.
**
** You should have received a copy of the GNU General Public License
** along with this program; if not, write to the Free Software
** Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
******************************************************************************/
#include "emboss.h"
/* @prog martfilters *********************************************************
**
** Reads and writes (returns) sequences
**
******************************************************************************/
int main(int argc, char **argv)
{
AjPSeqin seqin = NULL;
AjPMartquery mq = NULL;
AjPMartquery pmq = NULL;
AjPMartFilter filt = NULL;
AjPStr dataset = NULL;
AjPStr host = NULL;
AjPStr path = NULL;
AjPFile outf = NULL;
ajint iport = 0;
ajuint i;
AjPStr key1 = NULL;
AjPStr key2 = NULL;
const AjPStr value1 = NULL;
const AjPStr value2 = NULL;
embInit("martfilters", argc, argv);
host = ajAcdGetString("host");
path = ajAcdGetString("path");
iport = ajAcdGetInt("port");
dataset = ajAcdGetString("dataset");
outf = ajAcdGetOutfile("outfile");
mq = ajMartqueryNew();
seqin = ajSeqinNew();
ajMartAttachMartquery(seqin,mq);
ajMartSetMarthostS(seqin,host);
ajMartSetMartpathS(seqin,path);
ajMartSetMartport(seqin,iport);
/*
**
** This section retrieves the filters given a mart dataset name
** from a mart host.
*/
ajMartGetFilters(seqin, dataset);
ajMartfiltersParse(seqin);
key1 = ajStrNewC("name");
key2 = ajStrNewC("displayName");
pmq = ajMartGetMartqueryPtr(seqin);
filt = pmq->Filters;
for(i=0; i < filt->Nfilters; ++i)
{
value1 = ajTableFetchS(filt->Filters[i], key1);
value2 = ajTableFetchS(filt->Filters[i], key2);
ajFmtPrintF(outf,"%-40S %S\n",value1,value2);
}
ajStrDel(&key1);
ajStrDel(&key2);
ajStrDel(&host);
ajStrDel(&path);
ajStrDel(&dataset);
ajMartquerySeqinFree(seqin);
ajSeqinDel(&seqin);
ajFileClose(&outf);
embExit();
return 0;
}
|
