lowlevel.h

Go to the documentation of this file.

/* 
 *This file is part of biomcmc-lib, a low-level library for phylogenomic analysis.
 * Copyright (C) 2019-today  Leonardo de Oliveira Martins [ leomrtns at gmail.com;  http://www.leomartins.org ]
 *
 * biomcmc 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 (file "COPYING" or http://www.gnu.org/copyleft/gpl.html).
 */

#ifndef _biomcmc_lowlevel_h_
#define _biomcmc_lowlevel_h_

#include "config.h"

#include <stdio.h>      /* [ANSI C C89] */
#include <stdlib.h>     /* random number, searching, sorting, EXIT_SUCCESS [ANSI C C89] */
#include <string.h>     /* string manipulation [ANSI C C89] */
#include <stdarg.h>     /* Access to va_arg, va_list [ANSI C C89] */
#include <stdint.h>     /* standard integer types (int32_t typedef etc.) [C99]*/
#include <ctype.h>      /* char operation functions (e.g. isspace() ), case convertion [ANSI C C89] */
#include <math.h>       /* standard math functions (e.g. exp() ) [ANSI C C89] */
#include <float.h>      /* DBL_MAX_EXP, DBL_EPSILON constants (to avoid underflow etc) */
#include <time.h>       /* speed profiling(e.g. clock(), clock_gettime(), struct timespec ) [ANSI C C89] */
#include <unistd.h>     /* system values checking at runtime (e.g. sysconf() ) [POSIX C] */
#include <sys/time.h>   /* random seed (e.g. gettimeofday(), struct timeval) [POSIX C] */
#include <sys/times.h>  /* speed profiling in clock ticks (e.g. times() ) [POSIX.1 (or GNU extension?)] */ 
#include <sys/types.h>  /* pid_t for process ID, used together with unistd.h (may not be necessary) */
#include <unistd.h>     /* getpid() function, used together with sys/types.h (may not be necessary) */
#include <assert.h>    
//#include <sys/resource.h> // suggested by goptics (gpu), but don't seem needed
//#include <sys/stat.h>     // suggested by goptics (gpu), but don't seem needed


#include <libgen.h> /* standard XPG basename() - the one provided by string.h is a GNU extension, fails on macOSX*/

#ifdef _OPENMP
#include <omp.h>         /* OpenMP parallel threading library when available */
#endif

/* Global constants */

#define EXP_1 2.71828182845904523536028747135266 /* Euler's number */

#define true  1U 
#define false 0U 
#if defined(__GNUC__) && __GNUC__ >= 7
 #define attribute_FALLTHROUGH __attribute__ ((fallthrough));
#else
 #define attribute_FALLTHROUGH ((void)0);
#endif /* __GNUC__ >= 7 */

#define MIN(x,y) (((x)<(y)) ? (x) : (y))
#define MAX(x,y) (((x)>(y)) ? (x) : (y))
#define MOD(a)   (((a)>0)   ? (a) :(-a))


typedef unsigned char bool;

typedef struct hungarian_struct* hungarian;

struct hungarian_struct
{
  int **cost; 
  int size,  
      initial_cost, 
      final_cost;   
  int *col_mate, *unchosen_row, *slack_row, *row_mate, *parent_row;  
  double **dcost, initial_dcost, final_dcost; 
  double *row_dec_d, *col_inc_d, *slack_d; 
  int    *row_dec,   *col_inc,   *slack; /* aux vectors */
  bool is_double;
};

void *biomcmc_malloc (size_t size);

void *biomcmc_realloc (void *ptr, size_t size);

FILE *biomcmc_fopen (const char *path, const char *mode);

void biomcmc_error (const char *template, ...);

int compare_int_increasing (const void *a, const void *b);
int compare_int_decreasing (const void *a, const void *b);
int compare_uint64_increasing (const void *a, const void *b);
int compare_uint64_decreasing (const void *a, const void *b);
int compare_double_increasing (const void *a, const void *b);
int compare_double_decreasing (const void *a, const void *b);

int biomcmc_getline (char **lineptr, size_t *n, FILE *stream);

uint32_t biomcmc_levenshtein_distance (const char *s1, uint32_t n1, const char *s2, uint32_t n2, uint32_t cost_sub, uint32_t cost_indel, bool skip_borders);

/* Hungarian method for bipartite matching (assignment) */
hungarian new_hungarian (int size, bool is_double);
void hungarian_reset (hungarian p);
void hungarian_update_cost (hungarian p, int row, int col, void *cost); /* pointer since we decide if int/double later */
void del_hungarian (hungarian p);
void hungarian_solve (hungarian p, int this_size);

#endif