/* @funcstatic seqwords_TermsRead *********************************************
**
** Read the next Terms object from a file in embl-like format. The search
** terms are modified with a leading and trailing space.
**
** @param [r] inf [AjPFile] Input file stream
** @param [w] thys [AjPTerms*] Terms object
**
** @return [AjBool] True on succcess
** @@
*****************************************************************************/
static AjBool seqwords_TermsRead(AjPFile inf,
AjPTerms *thys)
{
AjPStr line =NULL; /* Line of text. */
AjPStr temp =NULL;
AjPList list_terms =NULL; /* List of keywords for a scop node*/
AjBool ok =ajFalse;
AjPStr type = NULL;
/* Memory management */
(*thys)=seqwords_TermsNew();
list_terms = ajListstrNew();
line = ajStrNew();
type = ajStrNew();
/* Read first line. */
ok = ajReadlineTrim(inf,&line);
while(ok && !ajStrPrefixC(line,"//"))
{
if(ajStrPrefixC(line,"XX"))
{
ok = ajReadlineTrim(inf,&line);
continue;
}
else if(ajStrPrefixC(line,"TY"))
{
ajFmtScanS(line, "%*s %S", &type);
if(ajStrMatchC(type, "SCOP"))
(*thys)->Type = ajSCOP;
else if(ajStrMatchC(type, "CATH"))
(*thys)->Type = ajCATH;
}
else if(ajStrPrefixC(line,"CL"))
{
ajStrAssignC(&(*thys)->Class,ajStrGetPtr(line)+3);
ajStrRemoveWhiteExcess(&(*thys)->Class);
}
else if(ajStrPrefixC(line,"AR"))
{
ajStrAssignC(&(*thys)->Architecture,ajStrGetPtr(line)+3);
ajStrRemoveWhiteExcess(&(*thys)->Architecture);
}
else if(ajStrPrefixC(line,"TP"))
{
ajStrAssignC(&(*thys)->Topology,ajStrGetPtr(line)+3);
ajStrRemoveWhiteExcess(&(*thys)->Topology);
}
else if(ajStrPrefixC(line,"FO"))
{
ajStrAssignC(&(*thys)->Fold,ajStrGetPtr(line)+3);
while(ajReadlineTrim(inf,&line))
{
if(ajStrPrefixC(line,"XX"))
break;
ajStrAppendC(&(*thys)->Fold,ajStrGetPtr(line)+3);
}
ajStrRemoveWhiteExcess(&(*thys)->Fold);
}
else if(ajStrPrefixC(line,"SF"))
{
ajStrAssignC(&(*thys)->Superfamily,ajStrGetPtr(line)+3);
while(ajReadlineTrim(inf,&line))
{
if(ajStrPrefixC(line,"XX"))
break;
ajStrAppendC(&(*thys)->Superfamily,ajStrGetPtr(line)+3);
}
ajStrRemoveWhiteExcess(&(*thys)->Superfamily);
}
else if(ajStrPrefixC(line,"FA"))
{
ajStrAssignC(&(*thys)->Family,ajStrGetPtr(line)+3);
while(ajReadlineTrim(inf,&line))
{
if(ajStrPrefixC(line,"XX"))
break;
ajStrAppendC(&(*thys)->Family,ajStrGetPtr(line)+3);
}
ajStrRemoveWhiteExcess(&(*thys)->Family);
}
else if(ajStrPrefixC(line,"TE"))
{
/* Copy and clean up term. */
temp = ajStrNew();
ajStrAssignC(&temp,ajStrGetPtr(line)+3);
ajStrRemoveWhiteExcess(&temp);
/* Append a leading and trailing space to search term*/
ajStrAppendK(&temp, ' ');
ajStrInsertC(&temp, 0, " ");
/* Add the current term to the list. */
ajListstrPush(list_terms,temp);
}
ok = ajReadlineTrim(inf,&line);
}
if(!ok)
{
/* Clean up. */
ajListstrFree(&list_terms);
ajStrDel(&line);
/* Return. */
return ajFalse;
}
/* Convert the AjPList of terms to array of AjPSeq's. */
if(!((*thys)->N=ajListstrToarray((AjPList)list_terms,&(*thys)->Keywords)))
ajWarn("Zero sized list of terms passed into seqwords_TermsRead");
/* Clean up. Free the list (not the nodes!). */
ajListstrFree(&list_terms);
ajStrDel(&line);
ajStrDel(&type);
return ajTrue;
}
static void acdrelations_readdatfile
(AjPFile inf,
PEdam *P)
{
AjPStr line = NULL;
const AjPStr tok = NULL;
const AjPStr subtok = NULL;
AjPStr strtmp = NULL;
AjPList strlist = NULL;
AjPStr acdtype = NULL;
AjPStr relations = NULL;
PEdamdat dattmp = NULL;
AjPList datlist = NULL;
if(!P)
ajFatal("Null arg error 1 in acdrelations_readdatfile");
if(!inf)
ajFatal("Null arg error 3 in acdrelations_readdatfile");
/* Allocate memory */
line = ajStrNew();
acdtype = ajStrNew();
relations = ajStrNew();
datlist = ajListNew();
/* Read data from file */
while(ajReadline(inf,&line))
{
/* Discard comment lines */
if(ajStrPrefixC(line,"#"))
continue;
/* Tokenise line, delimited by '|'.
Parse first token (ACD datatype ) */
ajStrAssignS(&acdtype, ajStrParseC(line, "|"));
/* Parse second token (EDAM relations: value ) */
ajStrAssignS(&relations, ajStrParseC(NULL, "|"));
/* Parse third token (attribute:value strings block) */
tok = ajStrParseC(NULL, "|");
/* Create new string list */
strlist = ajListstrNew();
/* Tokenise third token itself into tokens delimited by ' ' (space)
Parse tokens (individual attribute:value strings)*/
if((subtok=ajStrParseC(tok, ";")))
{
strtmp = ajStrNew();
ajStrAssignS(&strtmp, subtok);
ajStrRemoveWhite(&strtmp);
ajListstrPushAppend(strlist, strtmp);
while((subtok=ajStrParseC(NULL, ";")))
{
strtmp = ajStrNew();
ajStrAssignS(&strtmp, subtok);
ajStrRemoveWhite(&strtmp);
ajListstrPushAppend(strlist, strtmp);
}
}
/* Write PEdamdat structure & push onto list */
dattmp = ajEdamdatNew();
ajStrRemoveWhite(&acdtype);
ajStrAssignS(&dattmp->acdtype, acdtype);
ajStrAssignS(&dattmp->edam, relations);
dattmp->n = ajListstrToarray(strlist, &dattmp->acdattr);
ajListPushAppend(datlist, dattmp);
/* Clear nodes (but not strings) from string list */
ajListstrFree(&strlist);
}
/* Write PEdam structure */
((*P)->n) = ajListToarray(datlist, (void***) &((*P)->dat));
/* Free memory */
ajStrDel(&line);
ajStrDel(&acdtype);
ajStrDel(&relations);
ajListFree(&datlist);
return;
}
int main(int argc, char** argv)
{
AjPFile outf = NULL;
AjPFile cachef = NULL;
AjIList iterator = NULL;
AjPList aliases = NULL;
AjPList dbas = NULL;
AjPList species = NULL;
AjPStr alias = NULL;
AjPStr dbname = NULL;
AjPStr spname = NULL;
AjPStr svrname = NULL;
AjPStr svrurl = NULL;
AjPStr dbcurl = NULL;
AjPTime svrtime = NULL;
EnsEDatabaseadaptorGroup dbag = ensEDatabaseadaptorGroupNULL;
EnsPDatabaseadaptor dba = NULL;
EnsPDatabaseconnection dbc = NULL;
embInit("cacheensembl", argc, argv);
ensInit();
svrname = ajAcdGetString("servername");
outf = ajAcdGetOutfile("outfile");
cachef = ajAcdGetOutfile("cachefile");
dbcurl = ajStrNew();
svrurl = ajStrNew();
dbname = ajStrNew();
ajNamSvrGetUrl(svrname, &svrurl);
if(!svrurl)
ajFatal("Could not resolve server name '%S'.", svrname);
dbc = ensDatabaseconnectionNewUrl(svrurl);
ensRegistryLoadDatabaseconnection(dbc);
ensDatabaseconnectionDel(&dbc);
/* Write the server file header. */
svrtime = ajTimeNewTodayFmt("cachefile");
ajFmtPrintF(cachef, "# %S %D\n", ajFileGetNameS(cachef), svrtime);
ajTimeDel(&svrtime);
ajFmtPrintF(cachef,
"# Automatically generated by cacheensembl "
"for server '%S'.\n\n",
svrname);
/*
** Get all Ensembl Database Adaptor objects and write them as
** EMBOSS Database definitions.
*/
aliases = ajListstrNew();
dbas = ajListNew();
species = ajListstrNew();
ensRegistryRetrieveAllSpecies(species);
while(ajListstrPop(species, &spname))
{
ensRegistryGetAllDatabaseadaptors(ensEDatabaseadaptorGroupNULL,
spname,
dbas);
while(ajListPop(dbas, (void**) &dba))
{
dbag = ensDatabaseadaptorGetGroup(dba);
if(dbag == ensEDatabaseadaptorGroupNULL)
{
ajDebug("cacheensembl main got unexpected "
"Ensembl Database Adaptor Group %d.\n",
dbag);
continue;
}
ajStrAssignS(&dbname, ensDatabaseadaptorGetSpecies(dba));
if(dbag != ensEDatabaseadaptorGroupCore)
{
ajStrAppendC(&dbname, "_");
ajStrAppendC(&dbname, ensDatabaseadaptorGroupToChar(dbag));
}
dbc = ensDatabaseadaptorGetDatabaseconnection(dba);
ensDatabaseconnectionFetchUrl(dbc, &dbcurl);
if(outf)
ajFmtPrintF(outf, "%S\n", dbname);
ajFmtPrintF(cachef, "DBNAME %S [\n", dbname);
ajFmtPrintF(cachef, " release: \"%s\"\n", ensSoftwareGetVersion());
ajFmtPrintF(cachef, " server: \"%S\"\n", svrname);
ajFmtPrintF(cachef, " url: \"%S\"\n", dbcurl);
ajFmtPrintF(cachef, "]\n");
ajFmtPrintF(cachef, "\n");
if(dbag != ensEDatabaseadaptorGroupCore)
continue;
ensRegistryAliasFetchAllbySpecies(
ensDatabaseadaptorGetSpecies(dba),
aliases);
/*
** Format all aliases to lower case,
** sort them alphabetically and remove duplicates.
*/
iterator = ajListIterNew(aliases);
while(!ajListIterDone(iterator))
{
alias = ajListstrIterGet(iterator);
ajStrFmtLower(&alias);
}
ajListIterDel(&iterator);
ajListSortUnique(aliases,
cacheensembl_stringcompare,
cacheensembl_stringdelete);
alias = NULL;
if(ajListGetLength(aliases) > 0)
{
while(ajListstrPop(aliases, &alias))
{
/*
** Reject any aliases with other than alpha-numeric
** characters like white space.
*/
if(ajStrIsAlnum(alias))
ajFmtPrintF(cachef,
"ALIAS %S %S\n",
alias,
ensDatabaseadaptorGetSpecies(dba));
ajStrDel(&alias);
}
ajFmtPrintF(cachef, "\n");
}
/* Ensembl Database Adaptor objects *must not* be deleted. */
}
ajStrDel(&spname);
}
ajListstrFree(&aliases);
ajListFree(&dbas);
ajStrDel(&dbcurl);
ajStrDel(&svrurl);
ajStrDel(&dbname);
ajStrDel(&svrname);
ajFileClose(&outf);
ajFileClose(&cachef);
embExit();
return EXIT_SUCCESS;
}
int main(int argc, char **argv)
{
/* Variable declarations */
AjPStr dbname = NULL;
AjPStr svrname = NULL;
AjPFile outfile = NULL;
AjPResource resource = NULL;
AjPResourcein resourcein = NULL;
AjPStr resourceqry = NULL;
AjPStr type = NULL;
AjBool id;
AjBool qry;
AjBool all;
AjBool verbose;
AjPStr methods = NULL;
AjPStr release = NULL;
AjPStr comment = NULL;
AjPStr defined = NULL;
AjPList list = NULL;
AjPTagval tagval = NULL;
AjIList iter = NULL;
ajuint space = 0;
AjPList aliaslist = NULL;
AjIList aliter = NULL;
const AjPStr alias = NULL;
ajuint maxlen;
AjPStr truedbname = NULL;
AjBool isalias = ajFalse;
/* ACD processing */
embInit("dbtell", argc, argv);
dbname = ajAcdGetString("database");
svrname = ajAcdGetString("server");
verbose = ajAcdGetBoolean("full");
outfile = ajAcdGetOutfile("outfile");
ajStrAssignS(&truedbname, dbname);
ajNamAliasDatabase(&truedbname);
ajStrFmtLower(&truedbname);
if(!ajStrMatchS(dbname, truedbname))
isalias = ajTrue;
/* Application logic */
/* Check EMBOSS database information.
Write output file */
if(ajNamDbDetailsSvr(truedbname, svrname, &type, &id, &qry, &all,
&comment, &release, &methods, &defined))
{
if(isalias)
ajFmtPrintF(outfile, "# %S is an alias for %S defined in %S\n",
dbname, truedbname, defined);
else
ajFmtPrintF(outfile, "# %S is defined in %S\n",
truedbname, defined);
ajFmtPrintF(outfile, "# access levels id: %B query: %B all: %B\n\n",
id, qry, all);
ajFmtPrintF(outfile, "DBNAME %S [\n", truedbname);
if(ajStrGetLen(svrname))
list = ajNamDbGetAttrlistSvr(truedbname, svrname);
else
list = ajNamDbGetAttrlist(truedbname);
iter = ajListIterNewread(list);
while(!ajListIterDone(iter))
{
tagval = ajListIterGet(iter);
space = 15 - ajStrGetLen(ajTagvalGetTag(tagval));
ajFmtPrintF(outfile, " %S:%.*s\"%S\"\n",
ajTagvalGetTag(tagval),
space, " ",
ajTagvalGetValue(tagval));
ajTagvalDel(&tagval);
}
ajListIterDel(&iter);
ajListFree(&list);
ajFmtPrintF(outfile, "]\n");
if(verbose)
{
aliaslist = ajListNew();
ajNamListFindAliases(truedbname, aliaslist);
if(ajListGetLength(aliaslist))
{
ajFmtPrintF(outfile, "\n");
aliter = ajListIterNewread(aliaslist);
maxlen = 1;
while(!ajListIterDone(aliter))
{
alias = ajListIterGet(aliter);
if(MAJSTRGETLEN(alias) > maxlen)
maxlen = MAJSTRGETLEN(alias);
}
ajListIterDel(&aliter);
aliter = ajListIterNewread(aliaslist);
while(!ajListIterDone(aliter))
{
alias = ajListIterGet(aliter);
if(ajStrFindK(alias, ':') < 0)
ajFmtPrintF(outfile, "ALIAS %-*S %S\n",
maxlen, alias, truedbname);
}
ajListIterDel(&aliter);
}
ajListstrFree(&aliaslist);
}
ajStrDel(&type);
ajStrDel(&methods);
ajStrDel(&release);
ajStrDel(&comment);
ajStrDel(&defined);
}
else
{
/* try looking in DRCAT */
resourcein = ajResourceinNew();
resource = ajResourceNew();
ajFmtPrintS(&resourceqry, "drcat:%S", dbname);
ajResourceinQryS(resourcein, resourceqry);
if(ajResourceinRead(resourcein, resource))
{
ajFmtPrintF(outfile, "DBNAME %S [\n", dbname);
ajFmtPrintF(outfile, " comment: \"defined in DRCAT\"\n");
ajFmtPrintF(outfile, "]\n");
}
ajResourceinDel(&resourcein);
ajResourceDel(&resource);
ajStrDel(&resourceqry);
}
/* Memory clean-up and exit */
ajFileClose(&outfile);
ajStrDel(&truedbname);
ajStrDel(&dbname);
ajStrDel(&svrname);
embExit();
return 0;
}
