int main(int argc, char **argv)
{
/* Variable declarations */
AjPOutfile outfile = NULL;
AjBool obsolete = ajTrue;
AjPObo obo = NULL;
AjPObo oboanc = NULL;
AjPObo oboparent = NULL;
AjPOboall oboall = NULL;
AjPOboin oboin = NULL;
AjPStr oboqryup = NULL;
AjPStr oboqryanc = NULL;
AjPTable alltable = NULL;
AjPTable newtable = NULL;
AjPStr up = NULL;
AjPList uplist = NULL;
ajuint iterms = 0;
AjPStr topid = NULL;
AjPStr obodb = NULL;
AjBool saved = ajFalse;
/* ACD processing */
embInit("ontogetcommon", argc, argv);
oboall = ajAcdGetOboall("oboterms");
outfile = ajAcdGetOutobo("outfile");
obsolete = ajAcdGetBoolean("obsolete");
oboin = ajOboinNew();
oboparent = ajOboNew();
oboanc = ajOboNew();
uplist = ajListNew();
alltable = ajTablestrNew(600);
newtable = ajTablestrNew(600);
while(ajOboallNext(oboall, &obo))
{
saved = ajFalse;
if(!obsolete && ajOboIsObsolete(obo))
continue;
if(!iterms)
{
ajDebug("store id '%S'\n", ajOboGetId(obo));
ajStrAssignS(&obodb, ajOboGetDb(obo));
ajTablePut(alltable, ajStrNewS(ajOboGetId(obo)), NULL);
saved = ajTrue;
}
else
{
ajDebug("test id '%S'\n", ajOboGetId(obo));
if(ajTableMatchS(alltable, ajOboGetId(obo)))
{
ajDebug("keep id '%S'\n", ajOboGetId(obo));
ajTablePut(newtable, ajStrNewS(ajOboGetId(obo)), NULL);
saved = ajTrue;
}
}
if(saved)
ajStrAssignS(&topid, ajOboGetId(obo));
if(ajOboGetParents(obo, uplist)) /* that was the root */
{
while(ajListstrPop(uplist, &up))
{
ajDebug("up: '%S'\n", up);
ajFmtPrintS(&oboqryup, "%S-id:%S", ajOboGetDb(obo), up);
ajOboinQryS(oboin, oboqryup);
while(ajOboinRead(oboin, oboparent))
{
if(!obsolete && ajOboIsObsolete(oboparent))
continue;
if(!iterms)
{
ajDebug("store parent '%S'\n", ajOboGetId(oboparent));
ajTablePut(alltable, ajStrNewS(ajOboGetId(oboparent)),
NULL);
}
else
{
ajDebug("test parent '%S'\n", ajOboGetId(oboparent));
if(ajTableMatchS(alltable, ajOboGetId(oboparent)))
{
ajDebug("keep parent '%S'\n",
ajOboGetId(oboparent));
ajTablePut(newtable,
ajStrNewS(ajOboGetId(oboparent)),
NULL);
if(!saved)
{
ajStrAssignS(&topid, ajOboGetId(oboparent));
saved = ajTrue;
}
}
}
if(!ajOboGetParents(oboparent, uplist))
continue;
}
ajStrDel(&up);
}
}
if(iterms)
ajTableMergeAnd(alltable, newtable);
ajDebug("id: '%S' saved %u\n",
ajOboGetId(obo),
ajTableGetLength(alltable));
if(!ajTableGetLength(alltable))
ajDie("Query '%S' no matching ancestor found for obo term '%S:%S'",
ajOboallGetQryS(oboall), ajOboGetDb(obo), ajOboGetId(obo));
iterms++;
}
ajFmtPrintS(&oboqryanc, "%S-id:%S", obodb, topid);
ajOboinQryS(oboin, oboqryanc);
while(ajOboinRead(oboin, oboanc))
{
ajObooutWrite(outfile, oboanc);
}
/* Memory clean-up and exit */
ajOboallDel(&oboall);
ajOboinDel(&oboin);
ajOboDel(&oboanc);
ajOboDel(&oboparent);
ajListFree(&uplist);
ajTablestrFreeKey(&alltable);
ajTablestrFreeKey(&newtable);
ajStrDel(&oboqryup);
ajStrDel(&oboqryanc);
ajStrDel(&obodb);
ajStrDel(&topid);
ajOutfileClose(&outfile);
embExit();
return 0;
}
void ajResourceClear(AjPResource resource)
{
AjPReslink lnk = NULL;
AjPResquery qry = NULL;
AjPResterm resterm = NULL;
AjPStr ptr = NULL;
if(MAJSTRGETLEN(resource->Id))
ajStrSetClear(&resource->Id);
if(ajListGetLength(resource->Idalt))
while(ajListstrPop(resource->Idalt,&ptr))
ajStrDel(&ptr);
if(MAJSTRGETLEN(resource->Acc))
ajStrSetClear(&resource->Acc);
if(MAJSTRGETLEN(resource->Name))
ajStrSetClear(&resource->Name);
if(MAJSTRGETLEN(resource->Desc))
ajStrSetClear(&resource->Desc);
if(MAJSTRGETLEN(resource->Url))
ajStrSetClear(&resource->Url);
if(MAJSTRGETLEN(resource->Urllink))
ajStrSetClear(&resource->Urllink);
if(MAJSTRGETLEN(resource->Urlrest))
ajStrSetClear(&resource->Urlrest);
if(MAJSTRGETLEN(resource->Urlsoap))
ajStrSetClear(&resource->Urlsoap);
if(ajListGetLength(resource->Cat))
while(ajListstrPop(resource->Cat,&ptr))
ajStrDel(&ptr);
if(ajListGetLength(resource->Taxon))
while(ajListPop(resource->Taxon,(void**)&resterm))
ajRestermDel(&resterm);
if(ajListGetLength(resource->Edamdat))
while(ajListPop(resource->Edamdat,(void**)&resterm))
ajRestermDel(&resterm);
if(ajListGetLength(resource->Edamfmt))
while(ajListPop(resource->Edamfmt,(void**)&resterm))
ajRestermDel(&resterm);
if(ajListGetLength(resource->Edamid))
while(ajListPop(resource->Edamid,(void**)&resterm))
ajRestermDel(&resterm);
if(ajListGetLength(resource->Edamtpc))
while(ajListPop(resource->Edamtpc,(void**)&resterm))
ajRestermDel(&resterm);
if(ajListGetLength(resource->Xref))
while(ajListPop(resource->Xref,(void**)&lnk))
ajReslinkDel(&lnk);
if(ajListGetLength(resource->Query))
while(ajListPop(resource->Query,(void**)&qry))
ajResqueryDel(&qry);
if(ajListGetLength(resource->Example))
while(ajListstrPop(resource->Example,&ptr))
ajStrDel(&ptr);
if(MAJSTRGETLEN(resource->Db))
ajStrSetClear(&resource->Db);
if(MAJSTRGETLEN(resource->Setdb))
ajStrSetClear(&resource->Setdb);
if(MAJSTRGETLEN(resource->Full))
ajStrSetClear(&resource->Full);
if(MAJSTRGETLEN(resource->Qry))
ajStrSetClear(&resource->Qry);
if(MAJSTRGETLEN(resource->Formatstr))
ajStrSetClear(&resource->Formatstr);
if(MAJSTRGETLEN(resource->Filename))
ajStrSetClear(&resource->Filename);
ajStrDel(&resource->TextPtr);
resource->Count = 0;
resource->Fpos = 0L;
resource->Format = 0;
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;
}
