/***************************************************************************** Major portions of this software are copyrighted by the Medical College of Wisconsin, 1994-2000, and are released under the Gnu General Public License, Version 2. See the file README.Copyright for details. ******************************************************************************/ #ifndef _MCW_CS_STUFF_HEADER_ #define _MCW_CS_STUFF_HEADER_ /*----- Miscellaneous "computer science" stuff [legacy of when I was a professor of Computer Science -- RWCox] -----*/ #include #include #include #include #include "mcw_malloc.h" #ifdef __cplusplus extern "C" { #endif /***** Quicksort routines in various flavonoids *****/ extern void qsort_floatint ( int , float * , int * ) ; extern void qsort_doubleint ( int , double * , int * ) ; extern void qsort_intint ( int , int * , int * ) ; extern void qsort_floatfloat( int , float * , float * ) ; /* 30 Jan 2000 */ extern void qsort_floatstuff( int , float * , void ** ) ; /* 06 Feb 2000 */ extern void qsort_doublestuff( int, double * , void ** ) ; /* 18 Dec 2006 */ extern void qsort_intstuff ( int, int * , void ** ) ; /* 25 Mar 2010 */ extern void qsort_double ( int , double * ) ; /* 06 Jul 2010 */ /***** Quicksort-ish median *****/ extern float qmed_float ( int , float * ) ; extern float qmean_float ( int , float * ) ; /* 26 Mar 2013 */ extern void qmedmad_float ( int,float *,float *,float * ) ; /* 08 Mar 2001 */ extern void meansigma_float( int,float *,float *,float * ) ; /* 07 Dec 2006 */ extern void qmedmadbmv_float( int, float *, float *, float *, float * ) ; extern float centromean_float( int n , float *ar ) ; /* 01 Nov 2010 */ /***** Eigensolutions *****/ extern void symeig_double( int , double * , double * ) ; extern void symeigval_double( int , double * , double * ) ; extern void svd_double( int, int, double *, double *, double *, double * ) ; extern void svd_float ( int, int, float *, float *, float *, float * ) ; extern void set_svd_sort(int) ; extern int svd_desingularize( int m , int n , double *aa ) ; /* 11 Mar 2010 */ extern void symeig_3( double *, double *, int ) ; /* 30 Sep 2005 */ extern void symeig_2( double *, double *, int ) ; extern void symeig_forbid_23( int ) ; extern int symeig_irange( int n, double *a, double *e, int bb, int tt, int novec ) ; extern int first_principal_vectors( int n, int m, float *xx, int nvec, float *sval, float *uvec ) ; extern float principal_vector( int n , int m , int xtyp , void *xp , float *uvec , float *tvec , float *ws , unsigned short xran[] ) ; extern float mean_vector( int n , int m , int xtyp , void *xp , float *uvec ) ; extern void * pv_get_workspace( int n , int m ) ; #ifndef TYPEDEF_float_pair #define TYPEDEF_float_pair typedef struct { float a,b ; } float_pair ; #endif extern float_pair principal_vector_pair( int n, int m, int xtyp, void *xp, float *uvec, float *vvec, float *tvec, float *ws , unsigned short xran[] ) ; /***** Argument list mangling *****/ extern void addto_args( int , char ** , int * , char *** ) ; extern void append_string_to_args ( char *, int, char **, int *, char *** ) ; extern void prepend_string_to_args( char *, int, char **, int *, char *** ) ; /***** Misc stuff *****/ extern void get_laguerre_table( int , double ** , double ** ) ; /* 12 Mar 2000 */ extern int qhull_wrap( int , float * , int ** ) ; /* 07 Jun 2001 */ extern int sphere_voronoi_angles( int , float *, float *, float ** ) ; extern int sphere_voronoi_vectors( int , float *, float ** ) ; extern float cl1_solve ( int, int, float *, float **, float *,int ) ; /* cl1.c */ extern float cl1_solve_res( int, int, float *, float **, float *,int, float*,int ) ; /* cl1.c */ extern float cl2_solve ( int, int, float *, float **, float *,int ) ; /* cl2.c */ extern int powell_newuoa( int ndim , double *x , double rstart , double rend , int maxcall , double (*ufunc)(int,double *) ) ; extern void powell_set_mfac( float mm , float aa ) ; extern void powell_set_verbose( int ) ; extern int powell_newuoa_con( int ndim , double *x , double *xb, double *xt , int nrand, double rstart , double rend , int maxcall , double (*ufunc)(int,double *) ) ; extern int powell_newuoa_constrained( int ndim, double *x, double *cost , double *xbot, double *xtop , int nrand, int nkeep, int ntry , double rstart , double rend , int maxcall , double (*ufunc)(int,double *) ) ; extern char * approximate_number_string( double ) ; /* 16 Jan 2004 */ extern char * commaized_integer_string( long long ); /* 18 Mar 2010 */ extern char * string_substitute( char *src , char *targ , char *repl ) ; extern char * nice_time_string( int ) ; /* 18 Dec 2012 */ extern int is_a_number(char *) ; /* 03 Apr 2012 */ extern int is_an_int (char *) ; extern int strcmp_aboot( char * , char * ) ; /* 12 Mar 2007 */ extern char * afni_fgets( char *buf , int nbuf , FILE *fp ) ; /* 20 Dec 2011 */ extern void afni_fgets_setskip(int) ; #ifndef DEBLANK #define DEBLANK(cc) do{ if( (cc) != NULL ){ \ char *qc ; \ for( qc=(cc) ; *qc != '\0' ; qc++ ) \ if( isspace(*qc) ) *qc = '_' ; \ }} while(0) #endif typedef struct { double x; int Index; } Z_QSORT_DOUBLE; typedef struct { float x; int Index; } Z_QSORT_FLOAT; typedef struct { char * x; int Index; } Z_QSORT_STRING; typedef struct { int x; int Index; } Z_QSORT_INT; extern int compare_Z_IQSORT_STRING (Z_QSORT_STRING *a, Z_QSORT_STRING *b ); extern int compare_Z_IQSORT_DOUBLE (Z_QSORT_DOUBLE *a, Z_QSORT_DOUBLE *b ); extern int compare_Z_IQSORT_FLOAT (Z_QSORT_FLOAT *a, Z_QSORT_FLOAT *b ); extern int compare_Z_IQSORT_INT (Z_QSORT_INT *a, Z_QSORT_INT *b ); extern int compare_string (const void *a, const void *b ); extern int compare_double (double *a, double *b ); extern int compare_float (float *a, float *b ); extern int compare_int (int *a, int *b ); extern int compare_short (short *a, short *b ); extern int compare_char (char *a, char *b ); /***** little things to format values for nice printing *****/ typedef enum { CCALC_NOT_SET = 0, CCALC_DOUBLE = 1, CCALC_NICE, CCALC_INT, CCALC_FINT, CCALC_CINT, CCALC_CUSTOM } FORMAT_VALUE; #define AFNI_EOL '\n' extern char *format_value_4print (double value, int oform, char *formatstr); /***** pca calc *****/ double covariance(float *data_mat, double *cov_mat, unsigned char * row_mask, int num_rows, int num_cols, int norm, int remove_mean, int be_quiet); void pca (float *data_mat, unsigned char * row_mask, int num_rows, int num_cols, int be_quiet); double pca_fast3 (float *data_mat, int num_rows, int be_quiet, double *pca_vec, double *pca_eig); #ifdef __cplusplus } #endif #endif