+++ /dev/null
-
-#include <stdio.h>
-#include <stdlib.h>
-#include <math.h>
-#include <stdarg.h>
-
-
-#include "io_lib_header.h"
-#include "util_lib_header.h"
-#include "define_header.h"
-
-HaschT * hcreate ( int n_elements,struct Hasch_data * declare_data(struct Hasch_entry *), struct Hasch_data *free_data(struct Hasch_data *) )
- {
- HaschT *T;
- int a;
-
- n_elements=n_elements*2+1;
-
- T=vcalloc ( 1, sizeof (HaschT));
- T->ne=n_elements;
- T->p=vcalloc (n_elements,sizeof ( Hasch_entry*));
- for ( a=0; a<n_elements; a++)
- {
- T->p[a]=allocate_hasch_entry(NULL,DECLARE,declare_data, free_data);
- }
- return T;
- }
-HaschT *hdestroy (HaschT *T,struct Hasch_data * declare_data(struct Hasch_entry *), struct Hasch_data *free_data(struct Hasch_data *) )
-
-
- {
- int a;
- Hasch_entry *p, *pp;
-
- if ( T==NULL)return NULL;
-
- for (a=0; a< T->ne; a++)
- {
- p=T->p[a];
- while (p)
- {
- pp=p;
- p=p->n;
- allocate_hasch_entry(pp,FREE, declare_data, free_data);
- }
- }
- vfree (T->p);
- vfree ( T);
- return NULL;
- }
-
-
-Hasch_entry* hsearch (HaschT *T, int k, int action, struct Hasch_data * declare_data(struct Hasch_entry *), struct Hasch_data *free_data(struct Hasch_data *) )
-
-
- {
- /*action: FIND,ADD, REMOVE*/
- Hasch_entry *p, *pi;
- int h;
-
-
-
- /* find the key: k->h*/
-
- h=k%T->ne;
-
-
- if ( action==ADD || action==FIND)
- {
- p=pi=T->p[h];
- while (p && p->k!=k){p=p->n;}
- if (action==ADD && !p)
- {
- p=insert_hasch_entry_in_list (pi, NULL, NULL, declare_data, free_data);
- p->k=k;
- }
- else if (action==FIND && !p)p=NULL;
- return p;
- }
- else if ( action==REMOVE)
- {
- allocate_hasch_entry(hsearch ( T, k, FIND, declare_data, free_data), FREE, declare_data, free_data);
- return NULL;
- }
- return NULL;
- }
-
-
-Hasch_entry * extract_hasch_entry_from_list (Hasch_entry *e, struct Hasch_data * declare_data(struct Hasch_entry *), struct Hasch_data *free_data(struct Hasch_data *) )
-
-
- {
- /*extracts entry e and returns p, or next if is NULL*/
- Hasch_entry *p=NULL, *n=NULL;
-
- if (!e);
- else
- {
- p=e->p;
- n=e->n;
-
- if (p)p->n=n;
- if (n)n->p=p;
- e->p=e->n=NULL;
- }
- return e;
- }
-
-Hasch_entry * insert_hasch_entry_in_list (Hasch_entry *p, Hasch_entry *e, Hasch_entry *n, struct Hasch_data * declare_data(struct Hasch_entry *), struct Hasch_data *free_data(struct Hasch_data *) )
-
-
-{
- /*inserts entry e between entry p and entry n and returns e*/
-
- if (!e)e=allocate_hasch_entry (NULL,DECLARE, declare_data, free_data);
-
-
-
- if (!p && !n);
- else if ( !p)p=n->p;
- else if ( !n)n=p->n;
-
- e->p=p;
- if (p)p->n=e;
-
- e->n=n;
- if (n)n->p=e;
-
- return e;
- }
-
-Hasch_entry * allocate_hasch_entry (Hasch_entry *e, int action,struct Hasch_data * declare_data(struct Hasch_entry *), struct Hasch_data *free_data(struct Hasch_data *) )
-
-
-{
- static Hasch_entry *s;
- Hasch_entry *ns;
-
- if ( !s)s=vcalloc ( 1, sizeof (Hasch_entry));
-
- if ( action==DECLARE)
- {
- ns=s->p;
- e=extract_hasch_entry_from_list (s, declare_data, free_data);
- if ( e->free_data)(e->free_data)(e->data);
- e->declare_data=declare_data;
- e->free_data=free_data;
- e->declare_data (e);
- e->k=UNDEFINED;
- s=ns;
- }
- else if ( action==FREE)
- {
- extract_hasch_entry_from_list (e,declare_data, free_data );
- e->k=UNDEFINED;
- if ( e->free_data)e->data=(e->free_data)(e->data);
- e->free_data=NULL;
- e->declare_data=NULL;
- s=insert_hasch_entry_in_list (s, e, NULL, declare_data, free_data);
-
- }
- else if ( action==FREE_STACK)
- {
- while (s)
- {
- e=s->p;
- allocate_hasch_entry (s, FREE, declare_data,free_data);
- vfree (s);
- s=e;
- }
- }
- else crash ("Unknown MODE for allocate_hasch_entry\n");
- return e;
- }
-
-/*********************************************************************/
-/* */
-/* Get string key */
-/* */
-/* */
-/*********************************************************************/
-
-
-int string2key (char *s, Char_node *n)
-{
- static Char_node *root;
-
- if ( !root)root=declare_char_node (DECLARE);
-
- if ( n==NULL && s==NULL)
- {
- declare_char_node (FREE_STACK);
- }
- else if (n==NULL)
- {
- return string2key(s, root);
- }
- else if ( s[0]=='\0')
- {
- return n->key;
- }
- else
- {
- return string2key(s+1, (n->c[(int)s[0]])?(n->c[(int)s[0]]):(n->c[(int)s[0]]=declare_char_node (DECLARE)));
- }
- return 0;
-}
-
-Char_node * declare_char_node (int action)
-{
-static struct Char_node **heap;
-static int heap_size, free_heap, a;
-static int key;
- if ( action==DECLARE)
- {
- if ( free_heap==0)
- {
- free_heap=100;
-
- heap=vrealloc (heap,(heap_size+free_heap)*sizeof (struct Char_node *));
- for ( a=heap_size; a<heap_size+free_heap; a++)
- {
- (heap[a])=vcalloc ( 1, sizeof ( struct Char_node));
- (heap[a])->c=vcalloc ( 256, sizeof (Char_node*));
- (heap[a])->key=key++;
- }
- heap_size+=free_heap;
- }
- return heap[heap_size-(free_heap--)];
- }
- else if ( action==FREE_STACK)
- {
- for (a=0; a< heap_size; a++)
- {
- heap[a]->key=key++;
- vfree ( heap[a]->c);
- (heap[a])->c=vcalloc ( 256, sizeof (Char_node*));
- }
- free_heap=heap_size;
- return NULL;
- }
- return NULL;
-}
-
-/* old declare_char_node (too hungry)
-Char_node * declare_char_node (int action)
-{
- static int key;
- Char_node *cn;
- static Char_node *root;
-
- if ( action==DECLARE)
- {
- cn=vcalloc (1, sizeof (Char_node));
- cn->key=++key;
- cn->c=vcalloc (256, sizeof (Char_node *));
-
- }
- return cn;
-}
-*/
-/*********************************COPYRIGHT NOTICE**********************************/
-/*© Centro de Regulacio Genomica */
-/*and */
-/*Cedric Notredame */
-/*Tue Oct 27 10:12:26 WEST 2009. */
-/*All rights reserved.*/
-/*This file is part of T-COFFEE.*/
-/**/
-/* T-COFFEE 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.*/
-/**/
-/* T-COFFEE 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 Foobar; if not, write to the Free Software*/
-/* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA*/
-/*............................................... |*/
-/* If you need some more information*/
-/* cedric.notredame@europe.com*/
-/*............................................... |*/
-/**/
-/**/
-/* */
-/*********************************COPYRIGHT NOTICE**********************************/