Header file for libtpcmisc. More...

#include "tpcclibConfig.h"
#include <errno.h>
#include <stdio.h>
#include <math.h>
#include <time.h>
#include <stdlib.h>
#include <string.h>
#include <strings.h>
#include <ctype.h>
#include <stdint.h>
#include <unistd.h>

Data Structures
struct	IFT_KEY_AND_VALUE

struct	IFT

struct	INT_list

struct	INTEGER_LIST

struct	STR_TOKEN_LIST

Macros
#define	BRANCHING_O 0.999

#define	BRANCHING_C 0.998

#define	BRANCHING_Cu64 0.174

#define	BRANCHING_N 0.998

#define	BRANCHING_F 0.967

#define	BRANCHING_Ge 0.891

#define	BRANCHING_Ga 0.891

#define	BRANCHING_Rb 0.950

#define	HL_O15 2.05 /* 123 s */

#define	HL_N13 10.0

#define	HL_C11 20.4

#define	HL_F18 109.8

#define	HL_Ge68 396000.0 /* 275 d */

#define	HL_Ga68 68.0

#define	HL_Br75 98.0

#define	HL_Br76 978.33 /* 58700 s */

#define	HL_Cu62 9.7 /* 582 s */

#define	HL_Cu64 762.018 /* 12.7003 h */

#define	HL_Fe52 4980.0

#define	HL_Na22 1368000.0

#define	HL_O14 1.1818

#define	HL_Rb82 1.25 /* 75 s */

#define	HL_Zn62 558.0

#define	HL_I124 6013.44 /* 4.176 */

#define	M_LN2 0.69314718055994530942

#define	MAX_UNITS_LEN 31

#define	MAX_REGIONNAME_LEN 20

#define	MAX_REGIONSUBNAME_LEN 6

#define	MAX_STUDYNR_LEN 255

Enumerations
enum	{ TPCISOT_BR_75 , TPCISOT_BR_76 , TPCISOT_CU_62 , TPCISOT_CU_64 , TPCISOT_FE_52 , TPCISOT_GA_68 , TPCISOT_GE_68 , TPCISOT_NA_22 , TPCISOT_RB_82 , TPCISOT_ZN_62 , TPCISOT_F_18 , TPCISOT_C_11 , TPCISOT_N_13 , TPCISOT_O_15 , TPCISOT_O_14 , TPCISOT_I_124 , TPCISOT_UNKNOWN }

enum	{ CUNIT_UNKNOWN , CUNIT_CPS , CUNIT_COUNTS , CUNIT_KBQ_PER_ML , CUNIT_SEC_KBQ_PER_ML , CUNIT_PER_SEC , CUNIT_PER_MIN , CUNIT_ML_PER_ML , CUNIT_ML_PER_DL , CUNIT_ML_PER_ML_PER_MIN , CUNIT_ML_PER_DL_PER_MIN , CUNIT_UNITLESS , CUNIT_NCI_PER_ML , CUNIT_MBQ_PER_ML , CUNIT_BQ_PER_ML , CUNIT_UCI_PER_ML , CUNIT_UMOL_PER_MIN_PER_100G , CUNIT_MG_PER_MIN_PER_100G , CUNIT_UMOL_PER_MIN_PER_DL , CUNIT_MG_PER_MIN_PER_DL , CUNIT_PERCENTAGE , CUNIT_KCPS , CUNIT_MIN_KBQ_PER_ML , CUNIT_BQ , CUNIT_KBQ , CUNIT_MBQ , CUNIT_GBQ , CUNIT_NCI , CUNIT_UCI , CUNIT_MCI , CUNIT_PID , CUNIT_PIDM , CUNIT_PIDV , CUNIT_G_PER_ML , CUNIT_ML_PER_G , CUNIT_HU }

enum	{ TUNIT_UNKNOWN , TUNIT_SEC , TUNIT_MIN , TUNIT_UM , TUNIT_MM , TUNIT_CM , TUNIT_M , TUNIT_HOUR , TUNIT_MONTH , TUNIT_YEAR , TUNIT_MSEC }

enum	{ IFT_OK , IFT_FAULT , IFT_NOMEMORY , IFT_CANNOTREAD , IFT_CANNOTWRITE , IFT_UNKNOWNFORMAT , IFT_KEYNOTFOUND , IFT_NODATA , IFT_VALUENOTFOUND }

Functions
int	backupExistingFile (char filename, char backup_ext, char *status)

int	fileCopy (char filename1, char filename2, char *status)

float	branchingFraction (int isotope)

struct tm *	gmtime_r (const time_t t, struct tm tm)
	Convert time_t to GMT struct tm.

struct tm *	localtime_r (const time_t t, struct tm tm)
	Convert time_t to local time in struct tm.

time_t	timegm (struct tm *tm)
	Inverse of gmtime, converting struct tm to time_t.

char *	ctime_r_int (const time_t t, char buf)
	Convert calendard time t into a null-terminated string of the form YYYY-MM-DD hh:mm:ss, with length of 19 characters and the null.

int	isdate (char *str)

int	isdate2 (char str, char intdate)

int	isdate3 (char str, char intdate)

int	isdate4 (int dateint, int year, int month, int *day)

int	istime (char *str)

int	isdatetime (char str, char intdate)

int	get_datetime (char str, struct tm date, int verbose)

int	get_date (char str, struct tm date)

long int	math_div (long int a, long int b)

int	isleapyear (long int year)

long int	leaps_between (long int year1, long int year2)

void	time_to_tm (time_t totalsecs, int offset, struct tm *result)

double	tmDifference (struct tm tm1, struct tm tm0)

void	tmAdd (int s, struct tm *d)

int	dec_comma_is (char *str)

int	dec_separator (char *str)

void	dec_separator_change (char *str, int decsep)

double	atof_dpi (char *str)

int	dec_nr (char *str)

int	atof_with_check (char double_as_string, double result_value)

char *	strPtrToNextValue (char str, char *nxtp)

int	atoi_with_check (const char int_as_string, int result_value)

void	filenameRmPath (char *s)

int	filenameRmExtension (char *s)

void	filenameRmExtensions (char *s)

int	fnmatch (const char fname, const char key)

int	fncasematch (const char fname, const char key)

char *	filenameGetExtension (char *s)
	Get the last extension of a filename.

char *	filenameGetExtensions (char *s)
	Get all extensions of a filename.

char *	hlIsotopeCode (int isotope)

double	hlFromIsotope (char *isocode)

double	hl2lambda (double halflife)

double	hlLambda2factor (double lambda, double frametime, double framedur)

float	hlLambda2factor_float (float lambda, float frametime, float framedur)

char *	hlCorrectIsotopeCode (char *isocode)

int	hlIsotopeFromHalflife (double halflife)

const char *	libpet_idcrypt_version (void)

int	id_crypt (const char string, const char key, char *out, int decrypt)

void	iftSetStatus (IFT *ift, int status)

void	iftInit (IFT *ift)

void	iftEmpty (IFT *ift)

int	iftPut (IFT ift, char key, char value, char cmt_type, int verbose)

int	iftPutDouble (IFT ift, char key, double value, char *cmt_type, int verbose)

int	iftRead (IFT ift, char filename, int is_key_required, int verbose)

char *	iftReadValue (char filename, char keystr, int verbose)

int	iftWriteItem (IFT ift, int item, FILE fp, int verbose)

int	iftWrite (IFT ift, char filename, int verbose)

int	defRead (IFT ift, char filename, int verbose)

int	iftGet (IFT ift, char key, int verbose)

int	iftGetNth (IFT ift, char key, int n, int verbose)

int	iftFindNthKey (IFT ift, char str, int n, int verbose)

int	iftFindNthValue (IFT ift, char str, int n, int verbose)

int	iftGetKeyNr (IFT ift, const char key, int verbose)

int	iftGetFrom (IFT ift, int si, const char key, int verbose)

int	iftGetFullmatchFrom (IFT ift, int si, const char key, const char *value, int verbose)

int	iftGetFloatValue (IFT ift, int si, const char key, float *value, int verbose)

int	iftGetDoubleValue (IFT ift, int si, const char key, double *value, int verbose)

int	iftGetIntValue (IFT ift, int si, const char key, int *value, int verbose)

int	iftDeleteItem (IFT *ift, int item, int verbose)

int	iftReplaceNthValue (IFT ift, int item, char value, int verbose)

int	iftdup (IFT ift1, IFT ift2, int verbose)

void	intInit (INT_list *list)

void	intEmpty (INT_list *list)

int	intExpand (char text, INT_list list)

INT_list	intMerge (INT_list list1, INT_list list2)

int	_intexadd (INT_list *list, int a)

int	integerListInit (INTEGER_LIST *l)

int	integerListEmpty (INTEGER_LIST *l)

int	integerListAdd (INTEGER_LIST *l, int v, int ifnew)

int	integerListSort (INTEGER_LIST *l)

int	integerListAddFromString (const char s1, const char s2, INTEGER_LIST *l, const int ifnew)

int	integerListExpandFromString (const char s1, const char s2, INTEGER_LIST *l, const int ifnew)

int	temp_roundf (float e)

int	petCunitId (const char *unit)

int	petTunitId (const char *timeunit)

char *	petCunit (int cunit)

char *	petTunit (int tunit)

int	cunitFromFilename (char *fname)

int	tpcProcessStdOptions (const char s, int print_usage, int print_version, int verbose_level)

void	tpcProgramName (const char program, int version, int copyright, char prname, int n)

void	tpcPrintUsage (const char program, char text[], FILE *fp)

int	tpcHtmlUsage (const char program, char text[], const char *path)

void	tpcPrintBuild (const char program, FILE fp)

char *	strstr_noquotation (const char str1, const char str2)

int	strnCopyClean (char str1, const char str2, int maxlen)

void	str_token_list_init (STR_TOKEN_LIST *lst)

void	str_token_list_empty (STR_TOKEN_LIST *lst)

int	str_token_list_add (STR_TOKEN_LIST lst, char new_item)

int	str_token_list_del (STR_TOKEN_LIST *lst, int item)

int	str_token_list_read (const char filename, STR_TOKEN_LIST lst)

int	textfileReadLines (const char filename, STR_TOKEN_LIST lst)

int	readStrtokens (const char filename, char **toklist)

int	asciiCommentLine (const char line, int cont)

int	rnameSplit (char rname, char name1, char name2, char name3, int max_name_len)

int	rnameMatch (char rname, int rnr, char test_str)

int	rnameRmDots (char rname1, char rname2)

int	rnameCatenate (char rname, int max_rname_len, char name1, char name2, char name3, char space)

int	roinameExists (char *roiname)
	Verifies whether TAC name exists or not.

int	strTokenNr (const char str1, const char str2)

int	strTokenNCpy (const char str1, const char str2, int i, char *str3, int count)

char *	strTokenDup (const char s1, const char s2, int *next)

int	strChrCount (const char str1, const char str2)

int	strUppercaseCount (const char *s)

void	strReplaceChar (char *str, char c1, char c2)

size_t	strnlen (const char *s, size_t n)

size_t	strlcat (char dst, const char src, size_t dstsize)

size_t	strlcpy (char dst, const char src, size_t dstsize)

char *	strcasestr (const char haystack, const char needle)

int	strncpyCleanSpaces (char s1, const char s2, int maxlen)

int	strCleanSpaces (char *s)

char *	strEncodeForXML (const char *s)

void	strCleanForXML (char *s)

int	studynr_in_fname (char fname, char studynr)

int	studynr_from_fname (char fname, char studynr)

int	studynr_from_fname2 (char fname, char studynr, int force)

int	studynr_match (char studynr1, char studynr2)

int	studynr_validity_check2 (char *studynr, int zero_ok)

int	studynr_validity_check (char *studynr)

int	studynr_rm_zeroes (char *studynr)

int	studynr_to_lowercase (char *studynr)

int	little_endian ()

void	swap (void orig, void new, int size)

void	swabip (void *buf, long long int size)

void	swawbip (void *buf, long long int size)

void	swawip (void *buf, long long int size)

void	printf32bits (void *buf)

int	doubleMatch (const double v1, const double v2, const double lim)

int	doubleMatchRel (const double v1, const double v2, const double lim)

double	doubleMachEps ()

void	doubleCopy (double t, double s, const unsigned int n)

unsigned int	doubleMaxIndex (double *a, const unsigned int n)

double	doubleSum (double *a, const unsigned int n)

double	doubleMean (double *a, const unsigned int n)

int	doubleSpanPositives (double *a, const int n)

int	doubleCSpanPositives (double *a, const int n)

double	inverfc (double x)

void	statSortDouble (double *data, unsigned int n, int order)

void	statSortFloat (float *data, unsigned int n, int order)

Detailed Description

Header file for libtpcmisc.

Author: Vesa Oikonen

Definition in file libtpcmisc.h.

Macro Definition Documentation

◆ BRANCHING_C

#define BRANCHING_C 0.998

Isotope branching ratio

Definition at line 28 of file libtpcmisc.h.

Referenced by branchingFraction().

◆ BRANCHING_Cu64

#define BRANCHING_Cu64 0.174

Isotope branching ratio

Definition at line 30 of file libtpcmisc.h.

Referenced by branchingFraction().

◆ BRANCHING_F

#define BRANCHING_F 0.967

Isotope branching ratio

Definition at line 34 of file libtpcmisc.h.

Referenced by branchingFraction().

◆ BRANCHING_Ga

#define BRANCHING_Ga 0.891

Isotope branching ratio

Definition at line 38 of file libtpcmisc.h.

Referenced by branchingFraction().

◆ BRANCHING_Ge

#define BRANCHING_Ge 0.891

Isotope branching ratio

Definition at line 36 of file libtpcmisc.h.

Referenced by branchingFraction().

◆ BRANCHING_N

#define BRANCHING_N 0.998

Isotope branching ratio

Definition at line 32 of file libtpcmisc.h.

Referenced by branchingFraction().

◆ BRANCHING_O

#define BRANCHING_O 0.999

Isotope branching ratio

Definition at line 26 of file libtpcmisc.h.

Referenced by branchingFraction().

◆ BRANCHING_Rb

#define BRANCHING_Rb 0.950

Isotope branching ratio

Definition at line 40 of file libtpcmisc.h.

Referenced by branchingFraction().

◆ HL_Br75

#define HL_Br75 98.0

Isotope halflife in minutes; not verified from the reference

Definition at line 58 of file libtpcmisc.h.

◆ HL_Br76

#define HL_Br76 978.33 /* 58700 s */

Isotope halflife in minutes; not verified from the reference

Definition at line 60 of file libtpcmisc.h.

◆ HL_C11

#define HL_C11 20.4

Isotope halflife in minutes

Definition at line 47 of file libtpcmisc.h.

◆ HL_Cu62

#define HL_Cu62 9.7 /* 582 s */

Isotope halflife in minutes; not verified from the reference

Definition at line 62 of file libtpcmisc.h.

◆ HL_Cu64

#define HL_Cu64 762.018 /* 12.7003 h */

Isotope halflife in minutes; not verified from the reference

Definition at line 64 of file libtpcmisc.h.

◆ HL_F18

#define HL_F18 109.8

Isotope halflife in minutes

Definition at line 49 of file libtpcmisc.h.

◆ HL_Fe52

#define HL_Fe52 4980.0

Isotope halflife in minutes; not verified from the reference

Definition at line 66 of file libtpcmisc.h.

◆ HL_Ga68

#define HL_Ga68 68.0

Isotope halflife in minutes

Definition at line 53 of file libtpcmisc.h.

◆ HL_Ge68

#define HL_Ge68 396000.0 /* 275 d */

Isotope halflife in minutes

Definition at line 51 of file libtpcmisc.h.

◆ HL_I124

#define HL_I124 6013.44 /* 4.176 */

Isotope halflife in minutes; not verified from the reference

Definition at line 76 of file libtpcmisc.h.

◆ HL_N13

#define HL_N13 10.0

Isotope halflife in minutes

Definition at line 45 of file libtpcmisc.h.

◆ HL_Na22

#define HL_Na22 1368000.0

Isotope halflife in minutes; not verified from the reference

Definition at line 68 of file libtpcmisc.h.

◆ HL_O14

#define HL_O14 1.1818

Isotope halflife in minutes; not verified from the reference

Definition at line 70 of file libtpcmisc.h.

◆ HL_O15

#define HL_O15 2.05 /* 123 s */

Isotope halflife in minutes

Definition at line 43 of file libtpcmisc.h.

◆ HL_Rb82

#define HL_Rb82 1.25 /* 75 s */

Isotope halflife in minutes; not verified from the reference

Definition at line 72 of file libtpcmisc.h.

◆ HL_Zn62

#define HL_Zn62 558.0

Isotope halflife in minutes; not verified from the reference

Definition at line 74 of file libtpcmisc.h.

◆ M_LN2

#define M_LN2 0.69314718055994530942

ln(2)

Definition at line 90 of file libtpcmisc.h.

Referenced by hl2lambda(), and plasma_to_blood_ratio().

◆ MAX_REGIONNAME_LEN

#define MAX_REGIONNAME_LEN 20

Max length of Region name (+1)

Definition at line 154 of file libtpcmisc.h.

Referenced by cptReadOne(), cptrnameSplit(), csv2dft_a(), csv2dft_b(), csv2dft_linkset(), csv2dft_mat(), dft_fill_hdr_from_IFT(), dftGetPmodTitle(), dftRead(), dftRNameSimplify(), fitRead(), resRead(), rnameMatch(), rnameSplit(), roikbqRead(), and tsvRead().

◆ MAX_REGIONSUBNAME_LEN

#define MAX_REGIONSUBNAME_LEN 6

Max length of Region name subfield (+1)

Definition at line 158 of file libtpcmisc.h.

Referenced by clusterTACs(), cptReadOne(), csv2dft(), csv2dft_a(), csv2dft_b(), csv2dft_linkset(), csv2dft_mat(), dft_fill_hdr_from_IFT(), dftGetPmodTitle(), dftRead(), fitRead(), fitWrite(), idwcRead(), resRead(), resWrite(), tsvRead(), and xelRead().

◆ MAX_STUDYNR_LEN

#define MAX_STUDYNR_LEN 255

Max length of Study number (+1)

Definition at line 163 of file libtpcmisc.h.

Referenced by csv2dft_b(), csv2dft_mat(), dft_fill_hdr_from_IFT(), dftAllocateWithIMG(), dftRead(), ecat63ReadAllToImg(), ecat63ReadPlaneToImg(), fitRead(), img2sif(), imgGetAnalyzeHeader(), imgGetEcat63MHeader(), imgGetEcat7MHeader(), imgInfo(), imgReadAnalyze(), imgReadAnalyzeHeader(), imgReadNiftiHeader(), resRead(), sif2img(), studynr_from_fname2(), studynr_in_fname(), studynr_match(), and studynr_validity_check2().

◆ MAX_UNITS_LEN

#define MAX_UNITS_LEN 31

Max length of units string (+1), based on ECAT7 format

Definition at line 95 of file libtpcmisc.h.

Referenced by cptReadOne(), dftAllocateWithIMG(), and dftGetPmodTitle().

Enumeration Type Documentation

◆ anonymous enum

anonymous enum

isotope_code

Definition at line 79 of file libtpcmisc.h.

     {
  TPCISOT_BR_75, TPCISOT_BR_76, TPCISOT_CU_62, TPCISOT_CU_64, 
  TPCISOT_FE_52, TPCISOT_GA_68, TPCISOT_GE_68, TPCISOT_NA_22,
  TPCISOT_RB_82, TPCISOT_ZN_62, TPCISOT_F_18, TPCISOT_C_11,
  TPCISOT_N_13, TPCISOT_O_15, TPCISOT_O_14, TPCISOT_I_124,
  TPCISOT_UNKNOWN
};

◆ anonymous enum

anonymous enum

Data y units

Definition at line 99 of file libtpcmisc.h.

     {
  /*  0 */ CUNIT_UNKNOWN,
  /*  1 */ CUNIT_CPS,
  /*  2 */ CUNIT_COUNTS,
  /*  3 */ CUNIT_KBQ_PER_ML,
  /*  4 */ CUNIT_SEC_KBQ_PER_ML,
  /*  5 */ CUNIT_PER_SEC,
  /*  6 */ CUNIT_PER_MIN,
  /*  7 */ CUNIT_ML_PER_ML,
  /*  8 */ CUNIT_ML_PER_DL,
  /*  9 */ CUNIT_ML_PER_ML_PER_MIN,
  /* 10 */ CUNIT_ML_PER_DL_PER_MIN,
  /* 11 */ CUNIT_UNITLESS,
  /* 12 */ CUNIT_NCI_PER_ML,
  /* 13 */ CUNIT_MBQ_PER_ML,
  /* 14 */ CUNIT_BQ_PER_ML,
  /* 15 */ CUNIT_UCI_PER_ML,
  /* 16 */ CUNIT_UMOL_PER_MIN_PER_100G,
  /* 17 */ CUNIT_MG_PER_MIN_PER_100G,
  /* 18 */ CUNIT_UMOL_PER_MIN_PER_DL,
  /* 19 */ CUNIT_MG_PER_MIN_PER_DL,
  /* 20 */ CUNIT_PERCENTAGE,
  /* 21 */ CUNIT_KCPS,
  /* 22 */ CUNIT_MIN_KBQ_PER_ML,
  /* 23 */ CUNIT_BQ,
  /* 24 */ CUNIT_KBQ,
  /* 25 */ CUNIT_MBQ,
  /* 26 */ CUNIT_GBQ,
  /* 27 */ CUNIT_NCI,
  /* 28 */ CUNIT_UCI,
  /* 29 */ CUNIT_MCI,
  /* 30 */ CUNIT_PID,
  /* 31 */ CUNIT_PIDM,
  /* 32 */ CUNIT_PIDV,
  /* 33 */ CUNIT_G_PER_ML, // SUV unit
  /* 34 */ CUNIT_ML_PER_G, // SUV unit
  /* 35 */ CUNIT_HU // Hounsfield unit
};

◆ anonymous enum

anonymous enum

Data x units

Definition at line 138 of file libtpcmisc.h.

     {
  /*  0 */ TUNIT_UNKNOWN,
  /*  1 */ TUNIT_SEC,
  /*  2 */ TUNIT_MIN,
  /*  3 */ TUNIT_UM,
  /*  4 */ TUNIT_MM,
  /*  5 */ TUNIT_CM,
  /*  6 */ TUNIT_M,
  /*  7 */ TUNIT_HOUR,
  /*  8 */ TUNIT_MONTH,
  /*  9 */ TUNIT_YEAR,
  /* 10 */ TUNIT_MSEC
};

◆ anonymous enum

anonymous enum

Definitions for ift status message

Definition at line 249 of file libtpcmisc.h.

249 {IFT_OK, IFT_FAULT, IFT_NOMEMORY, IFT_CANNOTREAD, IFT_CANNOTWRITE,

250 IFT_UNKNOWNFORMAT, IFT_KEYNOTFOUND, IFT_NODATA, IFT_VALUENOTFOUND};

Function Documentation

◆ _intexadd()

int _intexadd	(	INT_list *	list,
		int	a )

extern

int _intexadd(int *list, int a) ; local function. Deprecated.

Parameters

list
a

Definition at line 92 of file intex.c.

                                     {
  int i, j, n;
 
  /* Check if list is yet empty */
  if(list->nr==0) {
    list->i=(int*)malloc(sizeof(int)); if(list->i==NULL) return(-1);
    /* Put the first integer to list and return */
    list->nr=1; list->i[0]=a; return(1);
  }
  n=list->nr;
  /* Check through the existing list */
  for(i=0; i<n; i++) {
    /* if it already is listed, just return */
    if(list->i[i]==a) return(0);
    /* make room for this integer */
    list->i=(int*)realloc(list->i, (n+1)*sizeof(int)); if(list==NULL) return(-1);
    if(list->i[i]>a) {for(j=n-1; j>=i; j--) list->i[j+1]=list->i[j]; break;}
  }
  list->i[i]=a; list->nr=n+1;
  return(1);
}

Referenced by intExpand().

◆ asciiCommentLine()

int asciiCommentLine	(	const char *	line,
		int *	cont )

extern

Check if ASCII text line starts with comment character '#'. Comment character is searched from the first non-space character (space characters here include spaces and tabs).

Returns: 1 if this is comment line and 0 if not.

Author: Vesa Oikonen

Parameters

line	Pointer to string containing one line of ASCII text file
cont	Optional pointer which is set to the index of line where the first non-space character after the comment character starts. If line does not start with comment character, then this will point to the first non-space character of the line. Enter NULL if not needed.

Definition at line 246 of file readfile.c.

  {
  if(cont!=NULL) *cont=0; 
  if(line==NULL) return 0;
  char *cptr=(char*)line;
  int i=strspn(cptr, " \t"); cptr+=i; if(cont!=NULL) *cont=i;
  if(*cptr!='#') return 0;
  if(cont==NULL) return 1;
  cptr++; i=strspn(cptr, " \t"); *cont+=(i+1);
  return 1;
}

Referenced by sifRead().

◆ atof_dpi()

double atof_dpi ( char * str )

extern

Replacement of atof(), which works whether string contains decimal dots or decimal commas. Possible commas are replaced by dots in the argument string.

Returns: Returns the double float.

Parameters

str	Pointer to string (not modified).

Definition at line 59 of file decpoint.c.

  {
  char *cptr;
  double f;
 
  if(str==NULL) return(nan(""));
  /* If string contains a dot, then use atof directly */
  if(strchr(str, '.')!=NULL) return(atof(str));
  /* Otherwise, convert all commas to dots */
  //while((cptr=strchr(str, ','))!=NULL) *cptr='.';
  cptr=strchr(str, ','); if(cptr!=NULL) *cptr='.';
  f=atof(str); if(cptr!=NULL) *cptr=',';
  return(f);
}

Referenced by atof_with_check(), csv2dft_a(), csv2dft_b(), dft_fill_hdr_from_IFT(), dftRead(), fitRead(), and resRead().

◆ atof_with_check()

int atof_with_check	(	char *	double_as_string,
		double *	result_value )

extern

Converts a string to float using atof(), but if its return value is zero this function checks that argument string actually contains a number. Result value is set to NaN if string was not valid value. Both decimal point and comma are accepted.

Returns: Returns 0 if successful, and 1 in case of an error.

Parameters

double_as_string	String which is converted to a double; not modified
result_value	Pointer to the double float; enter NULL, if not needed

Definition at line 107 of file decpoint.c.

  {
  char* cptr;
  double f;
 
  if(result_value!=NULL) *result_value=nan("");
  if(double_as_string==NULL) return(1);
  f=atof_dpi(double_as_string);
  if(f!=0.0) {if(result_value!=NULL) *result_value=f; return(0);}
  cptr=double_as_string;
  while(*cptr!=0 && (*cptr=='+' || *cptr=='-' || *cptr==' ')) cptr++;
  if(*cptr=='0') {if(result_value!=NULL) *result_value=f; return(0);}
  return(1);
}

Referenced by csv2dft_linkset(), csv2dft_mat(), dftRead(), and readEcat931Calibrationfile().

◆ atoi_with_check()

int atoi_with_check	(	const char *	int_as_string,
		int *	result_value )

extern

Converts a string to integer (int) using atoi, but this function verifies that argument string actually contains an integer number. String must end in NULL character. Exponentials are not accepted. Result value is set to 0 if string was not valid value.

Returns: Returns 0 if successful, and 1 in case of an error.

Parameters

int_as_string	String which is converted to a int; not modified
result_value	Pointer to the int; enter NULL, if not needed

Definition at line 238 of file decpoint.c.

  {
  int i, len;
 
  if(result_value!=NULL) *result_value=0;
  if(int_as_string==NULL) return(1);
  len=strlen(int_as_string); if(len<1) return(1);
  /* First character can be + or - */
  i=strspn(int_as_string, "+-");
  if(i>1) return(1);
  if(i==1 && len==1) return(1); 
  /* Check that (rest of) characters are digits */
  for( ; i<len; i++) if(!isdigit(int_as_string[i])) return(1);
  /* Convert to int */
  if(result_value!=NULL) *result_value=atoi(int_as_string);
  return(0);
}

Referenced by integerListAddFromString(), and pxlRead().

◆ backupExistingFile()

int backupExistingFile	(	char *	filename,
		char *	backup_ext,
		char *	status )

extern

Check if specified file exists; rename existing file to a backup file. If also backup file exists, then remove that.

Returns: Returns 0, if successful, and <>0 in case of an error.

Parameters

filename	Name of file which, if it exists, is renamed to a backup file
backup_ext	Extension for backup file; NULL will set the default ".bak" extension.
status	Pointer to a string (allocated for at least 64 chars) where error message or other execution status will be written; enter NULL, if not needed

Definition at line 14 of file backup.c.

  {
  char bakfile[FILENAME_MAX];
  int ret;
 
  // Check the input
  if(filename==NULL || strlen(filename)<1) {
    if(status!=NULL) sprintf(status, "invalid filename");
    return 1;
  }
 
  // Check if file exists; if not then no need to make any backup
  if(access(filename, 0) == -1) {
    if(status!=NULL) sprintf(status, "file does not pre-exist");
    return 0;
  }
  // Create filename for the backup file
  strlcpy(bakfile, filename, FILENAME_MAX);
  if(backup_ext==NULL) strlcat(bakfile, ".bak", FILENAME_MAX);
  else strlcat(bakfile, backup_ext, FILENAME_MAX);
  // If also backup file exists, then just delete it 
  if(access(bakfile, 0) != -1) {
    ret=remove(bakfile);
    if(ret!=0) {
      if(status!=NULL) sprintf(status, "cannot delete previous backup file");
      return 3;
    }
  }
  // Rename file
  ret=rename(filename, bakfile);
  if(ret!=0) {
    if(status!=NULL) sprintf(status, "cannot rename file as backup");
    return 5;
  }
  if(status!=NULL) sprintf(status, "file renamed as backup");
  return 0;
}  

Referenced by dftWrite(), fitWrite(), plot_fit_svg(), plot_fitrange_svg(), plot_svg(), plotdata(), and resWrite().

◆ branchingFraction()

float branchingFraction ( int isotope )

Branching fraction for specified isotope.

Returns: Returns the branching factor, or 0 in case branching fraction is unknown.

Parameters

isotope Isotope code; see hlIsotopeFromHalflife()

Definition at line 14 of file branch.c.

  {
  float bf=0.0;
  switch(isotope) {
    case TPCISOT_CU_64: bf=BRANCHING_Cu64; break;
    case TPCISOT_GA_68: bf=BRANCHING_Ga; break;
    case TPCISOT_GE_68: bf=BRANCHING_Ge; break;
    case TPCISOT_RB_82: bf=BRANCHING_Rb; break;
    case TPCISOT_F_18:  bf=BRANCHING_F; break;
    case TPCISOT_C_11:  bf=BRANCHING_C; break;
    case TPCISOT_N_13:  bf=BRANCHING_N; break;
    case TPCISOT_O_15:  bf=BRANCHING_O; break;
    case TPCISOT_BR_75:
    case TPCISOT_BR_76:
    case TPCISOT_CU_62:
    case TPCISOT_FE_52:
    case TPCISOT_NA_22:
    case TPCISOT_O_14:
    case TPCISOT_I_124:
    case TPCISOT_ZN_62:
    case TPCISOT_UNKNOWN:
    default:            bf=0.0;
  }
  return(bf);
}

Referenced by hrrtMakeCalHdr(), imgBranchingCorrection(), and imgGetMicropetHeader().

◆ ctime_r_int()

char * ctime_r_int	(	const time_t *	t,
		char *	buf )

extern

Convert calendard time t into a null-terminated string of the form YYYY-MM-DD hh:mm:ss, with length of 19 characters and the null.

This is a replacement of the thread-safe ctime_r function which converts to date and time in english format.

Author: Vesa Oikonen

Returns: Returns pointer to the string, or null in case of an error.

Parameters

t	Pointer to calendar time ; do not give pointer to int here, like &e7mhdr.scan_start_time
buf	Pointer to string where the date and time will be written. It must be pre-allocated for at least 20 characters.

Definition at line 110 of file datetime.c.

  {
  if(buf==NULL) return(NULL);
  buf[0]=(char)0;
  struct tm tm;
  if(!gmtime_r(t, &tm)) return(NULL);
  if(!strftime(buf, 20, "%Y-%m-%d %H:%M:%S", &tm)) return(NULL);
  return(buf);
}

Referenced by atnMake(), ecat7MHeaderToIFT(), ecat7PrintMainheader(), fitWrite(), imgGetEcat63MHeader(), imgInfo(), imgSetEcat63MHeader(), imgWriteEcat63Frame(), imgWriteFrame(), plotdata(), resWrite(), resWriteHTML_table(), and sifPrint().

◆ cunitFromFilename()

int cunitFromFilename ( char * fname )

extern

Tries to find calibration unit from filename.

Returns: Returns CUNIT, which is CUNIT_UNKNOWN if not successful.

See also: petCunitId, petTunitId, petCunit, petTunit

Parameters

fname Pointer to file name, where calibration unit is tried to be found.

Definition at line 242 of file petunits.c.

  {
  char *cptr;
 
  if(fname==NULL || strlen(fname)<3) return CUNIT_UNKNOWN;
  for(int i=0; i<2; i++) {
    if(i==0) /* First, look in the extension */
      {cptr=strrchr(fname, '.'); if(cptr==NULL) {cptr=fname; i++;}}
    else /* Then, look into whole filename */
      cptr=fname;
    if(strcasestr(cptr, "KBQ")!=NULL) return CUNIT_KBQ_PER_ML;
    if(strcasestr(cptr, "MBQ")!=NULL) return CUNIT_MBQ_PER_ML;
    if(strcasestr(cptr, "BQ")!=NULL) return CUNIT_BQ_PER_ML;
    if(strcasestr(cptr, "NCI")!=NULL) return CUNIT_NCI_PER_ML;
    if(strcasestr(cptr, "KCPS")!=NULL) return CUNIT_KCPS;
    if(strcasestr(cptr, "CPS")!=NULL) return CUNIT_CPS;
  }
  return CUNIT_UNKNOWN;
}

◆ dec_comma_is()

int dec_comma_is ( char * str )

extern

Checks whether argument string contains a decimal comma instead of dot.

Returns: Returns 1 if decimal comma is found and 0 if not found.

Parameters

str	Pointer to string (not modified).

Definition at line 14 of file decpoint.c.

  {
  if(strchr(str, '.')!=NULL) return(0);
  if(strchr(str, ',')!=NULL) return(1);
  return(0);
}

◆ dec_nr()

int dec_nr ( char * str )

extern

Returns the number of decimal places in the argument string, representing a floating point value. String can contain either decimal dots or commas.

Definition at line 81 of file decpoint.c.

{
  char *cptr;
  int n;
 
  if(str==NULL) return 0;
  /* Find the first dot or comma, but stop with E */
  cptr=str;
  while(*cptr!='.' && *cptr!=',') {
    if(*cptr==(char)0 || *cptr=='E' || *cptr=='e') return 0;
    cptr++;
  }
  if(*cptr==(char)0 || (*cptr!='.' && *cptr!=',')) return 0;
  /* Calculate the number of digits that follows */
  cptr++; n=0; while(cptr!=NULL && isdigit(*cptr)) {n++; cptr++;}
  return n;
}

Referenced by dftRead().

◆ dec_separator()

int dec_separator ( char * str )

extern

Checks whether argument string contains a decimal comma or dot, or neither.

Returns: Returns 0, if neither is found, 1 if dot, and 2 if comma is found.

Parameters

str	Pointer to string (not modified).

Definition at line 28 of file decpoint.c.

  {
  if(strchr(str, '.')!=NULL) return(1);
  if(strchr(str, ',')!=NULL) return(2);
  return(0);
}

◆ dec_separator_change()

void dec_separator_change	(	char *	str,
		int	decsep )

extern

Convert the first decimal separator to comma or dot, as required.

Parameters

str	Pointer to string (modified when necessary).
decsep	Requested decimal separator: 0=dot, 1=comma.

Definition at line 40 of file decpoint.c.

  {
  char *cptr;
  cptr=strchr(str, '.'); if(cptr!=NULL && decsep==1) {*cptr=','; return;}
  cptr=strchr(str, ','); if(cptr!=NULL && decsep==0) {*cptr='.'; return;}
  return;
}

◆ defRead()

int defRead	(	IFT *	ift,
		char *	filename,
		int	verbose )

extern

Read definition file, for example microPET header file, into IFT structure.

Returns: Returns 0 if ok. Sets ift->status.

See also: iftWrite, iftRead

Parameters

ift	Pointer to initiated but empty IFT
filename	Input filename
verbose	Verbose level; if zero, then only warnings are printed into stderr

Definition at line 321 of file iftfile.c.

  {
  int i, ret, nr=0, line=0, initial_key_nr=0, nonprintable=0;
  char *allfile, *cptr, *key_ptr, *value_ptr, *cmt_ptr;
  FILE *fp;
 
  /* Check function input */
  if(verbose>0) printf("%s(*ift, %s)\n", __func__, filename);
  if(ift==NULL) return(1);
  if(filename==NULL || strlen(filename)<1) {
    iftSetStatus(ift, IFT_FAULT); return(1);
  }
  if(ift->keyNr>0) initial_key_nr=ift->keyNr;
 
  /* Open file */
  if(strcasecmp(filename, "stdin")==0) {
    fp=stdin;
  } else {
    fp=fopen(filename, "r");
    if(fp==NULL) {iftSetStatus(ift, IFT_CANNOTREAD); return(2);}
  }
 
  /* Get file size */
  nr=nonprintable=0; while((ret=fgetc(fp))!=EOF) {
    if(iscntrl(ret) && ret!=13 && ret!=10 && ret!=9) {
      nonprintable=1; break;}
    nr++;
  }
  if(nr<2) {
    if(strcasecmp(filename, "stdin")!=0) fclose(fp);
    if(nonprintable>0) {
      /* File contains non-printable characters; maybe binary file */
      iftSetStatus(ift, IFT_UNKNOWNFORMAT);
    } else {
      /* File just din't have any content */
      iftSetStatus(ift, IFT_NODATA);
    }
    return(3);
  }
  if(verbose>1) printf("  the size of file is %d bytes\n", nr);
  if(nr>5000000) {
    if(strcasecmp(filename, "stdin")!=0) fclose(fp);
    iftSetStatus(ift, IFT_UNKNOWNFORMAT); return(3);
  }
  rewind(fp);
 
  /* Allocate memory for file contents */
  allfile=(char*)malloc((nr+1)*sizeof(char));
  if(allfile==NULL) {
    if(strcasecmp(filename, "stdin")!=0) fclose(fp);
    iftSetStatus(ift, IFT_NOMEMORY); return(4);
  }
 
  /* Read file contents and close the file */
  i=0; while((ret=fgetc(fp))!=EOF && i<nr) allfile[i++]=(char)ret;
  allfile[i]=(char)0;
  if(strcasecmp(filename, "stdin")!=0) fclose(fp);
 
  /* and then fill the list */
  /* separate the first line */
  cptr=strtok(allfile, "\n\r"); line=0;
  do {
    if(verbose>10) printf("line %d: '%s'\n", line, cptr);
    /* Remove initial spaces and tabs */
    i=strspn(cptr, " \t"); cptr+=i;
    if(strlen(cptr)<1) {cptr=strtok(NULL, "\n\r"); continue;}
    /* Check if line starts with a comment character */
    if((cmt_ptr=strchr("#!;%", cptr[0]))!=NULL) {
      /* save the whole line as a key */
      key_ptr=cptr; value_ptr=NULL;
      ret=iftPut(ift, key_ptr, value_ptr, " ", verbose-2);
      if(ret) {
        free(allfile); iftEmpty(ift); iftSetStatus(ift, IFT_FAULT);
        return(10+ret);
      }
      /* separate the next line */
      cptr=strtok(NULL, "\n\r"); line++;
      continue;
    }
    if(verbose>11) printf("  line %d: '%s'\n", line, cptr);
    /* Get the first line string representing key name */
    key_ptr=cptr; cptr=strchr(cptr, ' ');
    if(cptr!=NULL) {*cptr=(char)0; cptr++;}
    /* Rest of the line is the key value */
    if(cptr==NULL) value_ptr=NULL; else value_ptr=cptr;
    ret=iftPut(ift, key_ptr, value_ptr, " ", verbose-2);
    if(ret) {
      free(allfile); iftEmpty(ift); iftSetStatus(ift, IFT_FAULT);
      return(10+ret);
    }
    /* separate the next line */
    cptr=strtok(NULL, "\n\r"); line++;
  } while(cptr!=NULL);
  free(allfile);
 
  /* Did we actually get any data? */
  if(ift->keyNr<=initial_key_nr) {iftSetStatus(ift, IFT_NODATA); return(7);}
 
  /* Set parameter file type */
  ift->type=4;
 
  iftSetStatus(ift, IFT_OK);
  return(0);
}

Referenced by imgReadMicropetFrame(), and imgReadMicropetHeader().

◆ doubleCopy()

void doubleCopy	(	double *	t,
		double *	s,
		const unsigned int	n )

extern

Copy double values from the 2nd array to the first.

Parameters

t	Target array
s	Source array
n	Length of arrays

Definition at line 93 of file doubleutil.c.

  {
  unsigned int i;
  if(t==NULL || s==NULL || n<1) return;
  for(i=0; i<n; i++) t[i]=s[i];
}

◆ doubleCSpanPositives()

int doubleCSpanPositives	(	double *	a,
		const int	n )

extern

Returns the length of array consisting of non-positive (<=0 and NaN) values.

Returns: Returns the index of first positive value in the given array.

See also: doubleSpanPositives

Parameters

a	Pointer to the array.
n	Length of the array.

Definition at line 186 of file doubleutil.c.

  {
  if(a==NULL || n<1) return(0);
  int i=0;
  for(i=0; i<n; i++) if(a[i]>0.0) break;
  return(i);
}

◆ doubleMachEps()

double doubleMachEps ( )

extern

Estimates the machine epsilon, the upper bound on the relative error due to rounding in floating point arithmetic, within one order of magnitude of the true machine epsilon.

Standard C library should also have DBL_EPSILON in float.h.

Returns: Estimate of machine epsilon.

Author: Vesa Oikonen

See also: doubleMatchRel, doubleMatch

Definition at line 83 of file doubleutil.c.

{
  double macheps=1.0;
  do {macheps/=2.0;} while((1.0+macheps/2.0)!=1.0);
  return(macheps);
}

Referenced by householder_transform().

◆ doubleMatch()

int doubleMatch	(	const double	v1,
		const double	v2,
		const double	lim )

extern

Verifies that given two doubles have the same value inside given limits.

Values are considered to match also if both are NaNs.

Returns: 1 if match is found, and 0 if not.

Author: Vesa Oikonen

See also: doubleMatchRel, doubleMachEps

Parameters

v1	First value
v2	Second value
lim	Limit for absolute difference (if <0 then test will fail every time)

Definition at line 28 of file doubleutil.c.

  {
  if(isnan(v1) && isnan(v2)) return 1;
  if(isnan(v1) || isnan(v2)) return 0;
  if(v1==v2) return 1;
  if(isnan(lim) || lim<0.0) return 0;
  if(fabs(v1-v2)<=lim) return 1;
  return 0;
}

◆ doubleMatchRel()

int doubleMatchRel	(	const double	v1,
		const double	v2,
		const double	lim )

extern

Verifies that given two doubles have the same value inside given relative limit |2*(v1-v2)/(v1+v2)|.

Values are considered to match also if both are NaNs, but not if mean is zero (unless both are exactly zero).

Returns: 1 if match is found, and 0 if not.

Author: Vesa Oikonen

See also: doubleMatch, doubleMachEps

Parameters

v1	First value
v2	Second value
lim	Limit for relative difference (if <0 then test will fail every time)

Definition at line 55 of file doubleutil.c.

  {
  if(isnan(v1) && isnan(v2)) return 1;
  if(isnan(v1) || isnan(v2)) return 0;
  if(v1==v2) return 1;
  if(isnan(lim)) return 0;
  double mean;
  mean=0.5*(v1+v2); if(!isnormal(mean)) return 0;
  if(fabs((v1-v2)/mean)<=lim) return 1;
  return 0;
}

◆ doubleMaxIndex()

unsigned int doubleMaxIndex	(	double *	a,
		const unsigned int	n )

extern

Find the maximum value in given double array.

Returns: The index [0..n-1] of maximum value in array; 0 is returned also in case of errors.

Parameters

a	Pointer to double array.
n	Length of array.

Definition at line 111 of file doubleutil.c.

  {
  if(a==NULL || n<1) return(0);
  unsigned int i, mi=0;
  double mv=nan("");
  for(i=0; i<n; i++) if(isnan(mv) || a[i]>mv) {mv=a[i]; mi=i;}
  return(mi);
}

◆ doubleMean()

double doubleMean	(	double *	a,
		const unsigned int	n )

extern

Calculate the mean of values in given double array.

Returns: The mean of array values, or NaN in case of an error.

See also: doubleSum

Parameters

a	Pointer to double array.
n	Length of array.

Definition at line 148 of file doubleutil.c.

  {
  if(a==NULL || n<1) return(nan(""));
  double s=0.0;
  unsigned int i, ci=0;
  for(i=0; i<n; i++) if(!isnan(a[i])) {ci++; s+=a[i];}
  if(ci<1) return(nan(""));
  return(s/(double)ci);
}

◆ doubleSpanPositives()

int doubleSpanPositives	(	double *	a,
		const int	n )

extern

Returns the length of array consisting of only positive (>0 and not NaN) values.

Returns: Returns the index of first nonpositive value in the given array.

See also: doubleCSpanPositives

Parameters

a	Pointer to the array.
n	Length of the array.

Definition at line 168 of file doubleutil.c.

  {
  if(a==NULL || n<1) return(0);
  int i=0;
  for(i=0; i<n; i++) if(!(a[i]>0.0)) break;
  return(i);
}

◆ doubleSum()

double doubleSum	(	double *	a,
		const unsigned int	n )

extern

Calculate the sum of values in given double array.

Returns: The sum of array values.

See also: doubleMean

Parameters

a	Pointer to double array.
n	Length of array.

Definition at line 130 of file doubleutil.c.

  {
  double s=0.0;
  if(a==NULL || n<1) return(s);
  for(unsigned int i=0; i<n; i++) if(!isnan(a[i])) s+=a[i];
  return(s);
}

◆ fileCopy()

int fileCopy	(	char *	filename1,
		char *	filename2,
		char *	status )

extern

Copy file contents to another file. Existing file will be overwritten, to prevent it call backupExistingFile() before calling this function.

Returns: Returns 0 if successfull, otherwise >0.

Parameters

filename1	Name of file to be copied
filename2	Name of new file
status	Pointer to a string (allocated for at least 64 chars) where error message or other execution status will be written; enter NULL, if not needed

Definition at line 65 of file backup.c.

  {
  FILE *from, *to;
  char c;
 
  // Check the input
  if(filename1==NULL || filename2==NULL) {
    if(status!=NULL) sprintf(status, "invalid filename");
    return 1;
  }
  // Open the file1 for reading
  if((from=fopen(filename1, "rb"))==NULL) {
    if(status!=NULL) sprintf(status, "cannot open file for read");
    return 2;
  }
  // Open file2 for writing
  if((to=fopen(filename2, "wb"))==NULL) {
    if(status!=NULL) sprintf(status, "cannot open file for write");
    fclose(from); return 3;
  }
  // Copy the file
  while(!feof(from)) {
    c=fgetc(from);
    if(ferror(from)) {
      if(status!=NULL) sprintf(status, "cannot read from file");
      fclose(from); fclose(to); (void)remove(filename2); return 4;
    }
    if(!feof(from)) fputc(c, to);
    if(ferror(to)) {
      if(status!=NULL) sprintf(status, "cannot write to file");
      fclose(from); fclose(to); (void)remove(filename2); return 6;
    }
  }
  // Close files
  if(fclose(from)==EOF) {
    if(status!=NULL) sprintf(status, "cannot close file");
    fclose(to); return 7;
  }
  if(fclose(to)==EOF) {
    if(status!=NULL) sprintf(status, "cannot close file");
    return 8;
  }
  return 0;
}

◆ filenameGetExtension()

char * filenameGetExtension ( char * s )

extern

Get the last extension of a filename.

Extension(s) in pathname are not searched for. Note that pointer points to the original string.

Returns: Pointer to the filename extension starting with '.', or NULL if no extension is found.

Parameters

s	Pointer to string; string is not edited here.

Definition at line 139 of file filename.c.

  {
  if(s==NULL || strlen(s)<1) return((char*)NULL);
  /* Identify path */
  char *pptr=strrchr(s, '/'); if(pptr==NULL) pptr=strrchr(s, '\\');
  if(pptr==NULL) pptr=s; else pptr++;
  /* Search for the last '.' after the path */
  char *cptr=strrchr(pptr, '.'); if(cptr==NULL) return(cptr);
  /* If filename starts with '.', that is not counted as extension */
  if(strlen(pptr)==strlen(cptr)) return((char*)NULL);
  return(cptr);
}

◆ filenameGetExtensions()

char * filenameGetExtensions ( char * s )

extern

Get all extensions of a filename.

Extension(s) in pathname are not searched for. Note that pointer points to the original string.

Returns: Pointer to the filename extension starting with '.', or NULL if no extension is found.

Parameters

s	Pointer to string; string is not edited here.

Definition at line 164 of file filename.c.

  {
  if(s==NULL || strlen(s)<1) return((char*)NULL);
  /* Identify path */
  char *pptr=strrchr(s, '/'); if(pptr==NULL) pptr=strrchr(s, '\\');
  if(pptr==NULL) pptr=s; else pptr++;
  /* Search for the last '.' after the path, ignoring first character
     in case filename starts with '.' */
  char *cptr=strchr(pptr+1, '.');
  return(cptr);
}

◆ filenameRmExtension()

int filenameRmExtension ( char * s )

extern

Remove the last extension from file name.

Author: Vesa Oikonen

Returns: 1 if extension was found (and removed), 0 if not.

Parameters

s	Pointer to string

Definition at line 34 of file filename.c.

  {
  if(s==NULL || strlen(s)<1) return(0);
  char *cptr;
  cptr=strrchr(s, '.'); if(cptr==NULL) return(0);
  if(cptr[1]=='/' || cptr[1]=='\\') return(0);
  *cptr=(char)0;
  return(1);
}

Referenced by filenameRmExtensions(), tpcHtmlUsage(), tpcPrintBuild(), tpcPrintUsage(), and tpcProgramName().

◆ filenameRmExtensions()

void filenameRmExtensions ( char * s )

extern

Remove all extensions from file name.

Author: Vesa Oikonen

Parameters

s	Pointer to string

Definition at line 51 of file filename.c.

  {
  if(s==NULL || strlen(s)<1) return;
  while(filenameRmExtension(s)) {}
  return;
}

◆ filenameRmPath()

void filenameRmPath ( char * s )

extern

Remove path from file name.

Author: Vesa Oikonen

Parameters

s	Pointer to string

Definition at line 13 of file filename.c.

  {
  if(s==NULL || strlen(s)<1) return;
  char *cptr;
  cptr=strrchr(s, '/'); if(cptr==NULL) cptr=strrchr(s, '\\');
  if(cptr==NULL) return;
  cptr++;
  int i, n=strlen(cptr);
  for(i=0; i<n; i++, cptr++) s[i]=*cptr;
  s[i]=(char)0;
  return;
}

Referenced by tpcHtmlUsage(), tpcPrintBuild(), tpcPrintUsage(), and tpcProgramName().

◆ fncasematch()

int fncasematch	(	const char *	fname,
		const char *	key )

extern

Case-independent check whether string fname matches string key, which may contain wildcards ? and *.

Returns: 1 if strings do match and 0 if not.

Parameters

fname	filename that is evaluated
key	key string which may contain wildcards '?' and '*'

Definition at line 100 of file filename.c.

  {
  if(fname==NULL || key==NULL) return(0);
  char *key_ptr=NULL, *fname_ptr=NULL;
 
  while((*key)&&(*fname)) {
    if((*key=='?')||(toupper((int)*key)==toupper((int)*fname))) {
      key++; fname++;
    } else if(*key=='*') {
      if(toupper((int)*(key+1))==toupper((int)*fname)) {
        key_ptr=(char*)key++; fname_ptr=(char*)fname+1;
      } else {
        fname++;
        if(*(key+1)=='?') {key_ptr=(char*)key++; fname_ptr=(char*)fname;}
      }
    } else if((key_ptr!=NULL) && (*fname_ptr)) {
      return(fnmatch(key_ptr, fname_ptr));
    } else {
      return(0);
    }
  }
  if((*fname)&&(key_ptr!=NULL)) {return(fnmatch(key_ptr, fname_ptr));}
  else {if(*key=='*') key++; return(toupper((int)*key)==toupper((int)*fname));}
}

Referenced by dftFormat(), dftWrite(), and rnameMatch().

◆ fnmatch()

int fnmatch	(	const char *	fname,
		const char *	key )

extern

Check if string fname matches string key, which may contain wildcards ? and *.

Returns: 1 if strings do match and 0 if not.

Parameters

fname	filename that is evaluated
key	key string which may contain wildcards '?' and '*'

Definition at line 66 of file filename.c.

  {
  if(fname==NULL || key==NULL) return(0);
  char *key_ptr=NULL, *fname_ptr=NULL;
 
  while((*key)&&(*fname)) {
    if((*key=='?')||(*key==*fname)) {
      key++; fname++;
    } else if(*key=='*') {
      if(*(key+1)==*fname) {key_ptr=(char*)key++; fname_ptr=(char*)fname+1;
      } else {
        fname++;
        if(*(key+1)=='?') {key_ptr=(char*)key++; fname_ptr=(char*)fname;}
      }
    } else if((key_ptr!=NULL) && (*fname_ptr)) {
      return(fnmatch(key_ptr, fname_ptr));
    } else {
      return(0);
    }
  }
  if((*fname)&&(key_ptr!=NULL)) {return(fnmatch(key_ptr, fname_ptr));}
  else {if(*key=='*') key++; return(*key==*fname);}
}

Referenced by fncasematch(), and fnmatch().

◆ get_date()

int get_date	(	char *	str,
		struct tm *	date )

extern

Reads date from a standard string representation of date.

Returns: Returns 0 when successful, <>0 in case of an error.

Parameters

str	Pointer to string that contains date in one of the formats YYYY-MM-DD, DD/MM/YYYY, DD.MM.YYYY, DD/MM/YY, or DD.MM.YYYY. This string is not modified.
date	Pointer to allocated struct tm where the date is written

Definition at line 377 of file datetime.c.

  {
  char buf[32];
  int ret, n, YYYY=0, MM=0, DD=0;
 
  if(strlen(str)<8) return 1;
  if((ret=isdate(str))<=0) {
    strncpy(buf, str, 10); buf[10]=(char)0; 
  } else {
    ret=isdate2(str, buf);
    if(ret>0) ret=isdate3(str, buf);
  }
  if(ret>0) return 2;
  n=sscanf(buf, "%d-%d-%d", &YYYY, &MM, &DD);
  if(n!=3) return(3);
  date->tm_year=YYYY-1900; date->tm_mday=DD; date->tm_mon=MM-1;
  date->tm_hour=0; date->tm_min=0; date->tm_sec=0;  date->tm_isdst=-1;
  ret=strftime(buf, 32, "%Y-%m-%d %H:%M:%S", date);
  if(ret==0) return(4);
  return(0);
}

Referenced by resRead().

◆ get_datetime()

int get_datetime	(	char *	str,
		struct tm *	date,
		int	verbose )

extern

Reads time and date from a standard string representation of date and time.

Returns: Returns 0 when successful, <>0 in case of an error.

Parameters

str	Pointer to string that contains date and time in one of the formats YYYY-MM-DD hh:mm:ss, DD.MM.YYYY hh:mm:ss, or DD/MM/YY hh:mm:ss This string is not modified.
date	Pointer to allocated struct tm where date and time is written
verbose	Verbose level; if <=0, then nothing is printed into stdout

Definition at line 322 of file datetime.c.

  {
  char buf[32];
  int ret, n, YYYY=0, MM=0, DD=0, hh=0, mm=0, ss=0;
 
  if(verbose>0) printf("get_datetime(%s)\n", str);
  if(date==NULL) return(1);
  ret=isdatetime(str, buf); if(ret!=0) return(ret);
  n=sscanf(buf, "%d-%d-%d %d:%d:%d", &YYYY, &MM, &DD, &hh, &mm, &ss);
  if(n!=6) return(40);
  date->tm_year=YYYY-1900; date->tm_mday=DD; date->tm_mon=MM-1;
  date->tm_hour=hh; date->tm_min=mm; date->tm_sec=ss; 
  date->tm_isdst=-1;
#ifdef HAVE_TM_GMTOFF
  date->tm_gmtoff=0L;
#endif
  if(verbose>2) {
    printf("tm_mday=%d\n", date->tm_mday);
    printf("tm_mon=%d\n", date->tm_mon);
    printf("tm_year=%d\n", date->tm_year);
    printf("tm_hour=%d\n", date->tm_hour);
    time_t t=timegm(date);
    printf(" tm_hour=%d\n", date->tm_hour);
    printf("time_t=%ld\n", (long)t);
    printf("time_t_int=%d\n", (int)t);
  }
//  if(mktime(date)==-1 && date->tm_year<70) {
  if(timegm(date)==-1 && date->tm_year<70) {
    if(verbose>2) printf("tm_mday=%d\n", date->tm_mday);
    if(verbose>1) printf("mktime() fails\n");
    /* Date and time format was previously verified, so lets put
       1970-01-01 01:00:00 instead */  
    date->tm_year=70; date->tm_mday=1; date->tm_mon=0;
    date->tm_hour=1; date->tm_min=0; date->tm_sec=0; 
    date->tm_yday=0;
    date->tm_isdst=-1;
    mktime(date);
  }
  if(verbose>1) printf(" tm_hour=%d\n", date->tm_hour);
  return(0);
}

Referenced by ecat7EditMHeader(), fitRead(), resRead(), and selectEcat931Calibrationfile().

◆ gmtime_r()

struct tm * gmtime_r	(	const time_t *	t,
		struct tm *	tm )

extern

Convert time_t to GMT struct tm.

This version of gmtime_r function is here for systems (at least Windows) where it is not defined. Uses gmtime, which is threadsafe in Windows.

Returns: Pointer to struct tm, or null in case of an error.

Parameters

t	Pointer to time_t; do not give pointer to int here, like &e7mhdr.scan_start_time
tm	Pointer to struct tm, to be filled here

Definition at line 22 of file datetime.c.

  {
  struct tm *ltm=gmtime(t);
  if(ltm==NULL || tm==NULL) return(NULL);
  *tm=*ltm; tm->tm_isdst=-1;
  return tm;
}

Referenced by ctime_r_int(), ecat63AddImg(), ecat63WriteAllImg(), ecat7EditMHeader(), ecat7ReadMainheader(), ecatCopy7to63mainheader(), imgSetAnalyzeHeader(), imgSetEcat63MHeader(), imgWriteAnalyze(), res2ift(), sifWrite(), time_to_tm(), and timegm().

◆ hl2lambda()

double hl2lambda ( double halflife )

extern

Calculates the isotope lambda from specified halflife.

Parameters

halflife halflife time value

Returns: A negative value is returned in case of error.

Definition at line 84 of file halflife.c.

                                  {
  if(halflife>0.0) return(M_LN2/halflife); else return(-1.0); // M_LN2=log(2.0)
}

Referenced by atnMake(), dftDecayCorrection(), img2sif(), imgDecayCorrection(), imgSetDecayCorrFactors(), and noiseSD4Simulation().

◆ hlCorrectIsotopeCode()

char * hlCorrectIsotopeCode ( char * isocode )

extern

Check that isotope code, e.g. F-18, is in valid format, containing '-' and in this order. Returns the correct isotope code.

Parameters

isocode Pointer to string "C-11", "11c" etc; contents of this string is not changed, and this is not returned in any case

Returns: pointer to correct isotope code, and NULL if it was not valid and could not be corrected.

Definition at line 141 of file halflife.c.

                                          {
  int i, ok=0, n, mass_nr=0, ic_mass_nr=0;
  char *cptr, atom[3], ic_atom[3];
 
  /* Check, if isocode can be found in the list */
  i=ok=0;
  while(isotope_code[i]!=0) {
    if(strcasecmp(isotope_code[i], isocode)==0) {
      ok=1;
  break;
    }
    i++;
  }
  if(ok==1) return(isotope_code[i]);
 
  /* Try to figure out what it is */
  /* Separate the atom name and mass number */
  n=strcspn(isocode, "-1234567890");
  if(n>2) { /* cannot be */
    return(NULL);
  } else if(n>0) { /* start with atom name */
    strncpy(atom, isocode, n); atom[n]=(char)0;
    mass_nr=atoi(isocode+n); if(mass_nr<0) mass_nr=-mass_nr;
  } else { /* starts with mass number */
    mass_nr=atoi(isocode);
    cptr=isocode; while(isdigit((int)cptr[0])) cptr++;
    if(strlen(cptr)>2) return(NULL);
    strcpy(atom, cptr); 
  }
  /* Check if it matches any of the listed isotopes */
  i=ok=0;
  while(isotope_code[i]!=0) {
    /* Separate the atom name and mass number from the listed isotope */
    n=strcspn(isotope_code[i], "-1234567890");
    strncpy(ic_atom, isotope_code[i], n); ic_atom[n]=(char)0;
    ic_mass_nr=atoi(isotope_code[i]+n+1);
    /* Check the atom name */
    if(strcasecmp(ic_atom, atom)!=0) {i++; continue;}
    /* Check the mass number, if given */
    if(mass_nr>0 && ic_mass_nr!=mass_nr) {i++; continue;}
    /* Match was found! */
    ok=1; break;
  }
  if(ok==1) return(isotope_code[i]); else return(NULL);
}

Referenced by hlFromIsotope(), and sif2dft().

◆ hlFromIsotope()

double hlFromIsotope ( char * isocode )

extern

Identify the isotope from the specified isotope code string and return the halflife (min). This function checks the validity of the isotope string using hlCorrectIsotopeCode(), but does not change it in any way.

Returns: A negative value is returned in case of error.

Parameters

isocode Pointer to string "C-11", "18f" etc. This argument is not changed.

Definition at line 55 of file halflife.c.

  {
  char *corrected_isocode;
  int i = 0;
  int ok = 0;
 
  /* Validate the isotope string */
  corrected_isocode=hlCorrectIsotopeCode(isocode);
  if(corrected_isocode==NULL) return(-1.0);
 
  /* Determine the isotope and return the half-life */
  while(isotope_code[i]!=0) {
    if(strcmp(isotope_code[i], corrected_isocode)==0) {ok=1; break;}
    i++;
  }
  if(ok==1) return(isotope_halflife[i]);
  else return(-2.0);
} 

Referenced by dftDecayCorrection(), imgReadAnalyzeHeader(), imgReadModelingData(), imgReadNiftiHeader(), noiseSD4SimulationFromDFT(), and sif2img().

◆ hlIsotopeCode()

char * hlIsotopeCode ( int isotope )

extern

Isotope code as a string, based on isotope list number.

Parameters

isotope index of PET isotope in the list in halflife.c

Returns: pointer to static string or "unknown".

Definition at line 36 of file halflife.c.

                                 {
  static char unknown_isotope[]="unknown";
  int n=0;
  
  while(isotope_code[n]!=0) n++;
  if(isotope<0 || isotope>n-1) return(unknown_isotope);
  else return(isotope_code[isotope]);
}

Referenced by dftDecayCorrection(), and imgIsotope().

◆ hlIsotopeFromHalflife()

int hlIsotopeFromHalflife ( double halflife )

extern

Identify the isotope based on its halflife (in minutes).

Parameters

halflife Half-life in minutes

Returns: the isotope list number, or negative value if not identified.

Definition at line 195 of file halflife.c.

                                           {
  int i=0, ok=0;
  if(halflife<=0.01) return(-1);
  while(isotope_halflife[i]>0.0) {
    if( fabs(halflife/isotope_halflife[i]-1.0)<0.05 ) {ok=1; break;}
    i++;
  }
  if(ok==1) return(i);
  else return(-2.0);
}

Referenced by dftDecayCorrection(), imgBranchingCorrection(), and imgIsotope().

◆ hlLambda2factor()

double hlLambda2factor	(	double	lambda,
		double	frametime,
		double	framedur )

extern

Calculate the decay correction factor for specified isotope lambda.

Parameters

lambda	Negative lambda removes decay correction
frametime	Frame start time, or mid time if framedur<=0
framedur	If unknown, set <0 and give mid time for frametime

Returns: A negative value is returned in case of error.

Definition at line 98 of file halflife.c.

                                                                         {
  double cf, ff;
 
  if(frametime<0) return(-1.0);
  cf=exp(lambda*frametime);
  if(framedur>1.0E-5) {
    ff=fabs(lambda)*framedur/(1.0-exp(-fabs(lambda)*framedur));
    if(lambda<0.0) cf/=ff; else cf*=ff;
  }
  return(cf);
}

Referenced by atnMake(), dftDecayCorrection(), img2sif(), and noiseSD4Simulation().

◆ hlLambda2factor_float()

float hlLambda2factor_float	(	float	lambda,
		float	frametime,
		float	framedur )

extern

Calculate the decay correction factor for specified isotope lambda. Version for floats (mainly image data).

Parameters

lambda	Negative lambda removes decay correction
frametime	Frame start time, or mid time if framedur<=0
framedur	If unknown, set <0 and give mid time for frametime

Returns: A negative value is returned in case of error.

Definition at line 118 of file halflife.c.

                                                                           {
  float cf, ff;
 
  if(frametime<0) return(-1.0);
  cf=exp(lambda*frametime);
  if(framedur>1.0E-5) {
    ff=fabs(lambda)*framedur/(1.0-exp(-fabs(lambda)*framedur));
    if(lambda<0.0) cf/=ff; else cf*=ff;
  }
  return(cf);
}  

Referenced by imgDecayCorrection(), and imgSetDecayCorrFactors().

◆ id_crypt()

int id_crypt	(	const char *	string,
		const char *	key,
		char *	out,
		int	decrypt )

extern

Scramble characters in ASCII range 32-126 using the Vigenere Cipher. Other characters are discarded.

Returns: Returns 0 if successful

Parameters

string	Original string to be encrypted/decrypted
key	Keyword string
out	Encrypted/decrypted string
decrypt	Set to 1 when decrypting, or to 0 when encrypting

Definition at line 23 of file idcrypt.c.

  {
    char *keystr;
    unsigned int len, r;
 
    len=strlen(string);
    if(len==0) return 0;
 
    keystr = malloc(len);
    if(!keystr) return 1;
 
    if(len>strlen(key)) {
        for(r=0;r<len;r++) {
            keystr[r]=key[r%strlen(key)]-32;
            if(keystr[r]>94) keystr[r]=94;
        }
    } else {
        for(r=0;r<len;r++) {
            keystr[r]=key[r]-32;
            if(keystr[r]>94) keystr[r]=94;
        }
    }
 
    for(r=0;r<len;r++) {
        int c=(unsigned char)string[r]-32;
        if(c>94) c=94;
        if(decrypt) {
            c=c-keystr[r];
            if(c<0) c=(95-(-c)%95);
        } else {
            c=(c+keystr[r])%95;
        }
        out[r]=c+32;
    }
    out[r]=0;
    free(keystr);
    return 0;
}

◆ iftDeleteItem()

int iftDeleteItem	(	IFT *	ift,
		int	item,
		int	verbose )

extern

Remove the specified item from IFT.

Returns: 0 if ok.

See also: iftPut, iftEmpty, iftReplaceNthValue

Parameters

ift	Pointer to IFT structure.
item	Index [0..keyNr-1] of key and value to delete.
verbose	Verbose level; if zero, then only warnings are printed into stderr

Definition at line 169 of file ift.c.

  {
  int i;
 
  if(verbose>0) printf("%s(*ift, %d)\n", __func__, item);
  if(ift==NULL) return(1);
  iftSetStatus(ift, IFT_FAULT);
  if(ift==NULL) {return(1);}
  if(item<0 || item>=ift->keyNr) {return(2);}
 
  if(ift->item[item].key!=NULL) free(ift->item[item].key);
  if(ift->item[item].value!=NULL) free(ift->item[item].value);
  ift->item[item].key=ift->item[item].value=NULL;
  for(i=item+1; i<ift->keyNr; i++) {
    ift->item[i-1].type=ift->item[i].type;
    ift->item[i-1].sw=ift->item[i].sw;
    ift->item[i-1].key=ift->item[i].key;
    ift->item[i-1].value=ift->item[i].value;
    ift->item[i].key=ift->item[i].value=NULL;
  }
  ift->keyNr--;
  iftSetStatus(ift, IFT_OK);
  return(0);
}

Referenced by ecat7CopyHeadersNoQuant().

◆ iftdup()

int iftdup	(	IFT *	ift1,
		IFT *	ift2,
		int	verbose )

extern

Make a copy (duplicate) of IFT structure.

Returns: Returns 0 (IFT_OK) when successful, otherwise an appropriate ift status code.

See also: iftInit, iftRead, iftPut, iftEmpty

Parameters

ift1	Pointer to IFT structure to be copied.
ift2	Pointer to initiated IFT structure; any previous contents are deleted.
verbose	Verbose level; if zero, then only warnings are printed into stderr

Definition at line 235 of file ift.c.

  {
  int ret, li;
 
  if(verbose>0) printf("%s(*ift1, *ift2)\n", __func__);
  /* Check the input */
  if(ift1==NULL || ift2==NULL) return IFT_FAULT;
 
  /* Empty the new IFT */
  iftEmpty(ift2);
 
  /* Copy the contents */
  ift2->type=ift1->type;
  ift2->status=ift1->status;
  for(li=0; li<ift1->keyNr; li++) {
    ret=iftPut(ift2, ift1->item[li].key, ift1->item[li].value, NULL, verbose);
    if(ret!=IFT_OK) {iftEmpty(ift2); return(ret);} 
    ift2->item[li].type=ift1->item[li].type;
    ift2->item[li].sw=ift1->item[li].sw;
  }
  // keyNr was set by iftPut()
  return IFT_OK;
}

Referenced by imgCopyhdr().

◆ iftEmpty()

void iftEmpty ( IFT * ift )

extern

Free memory allocated for IFT. All contents are destroyed.

See also: iftInit

Parameters

ift	Pointer to target IFT structure.

Definition at line 60 of file ift.c.

  {
  if(ift==NULL) return;
  for(int i=0; i<ift->_memNr; i++) {
    if(ift->item[i].key!=NULL) free(ift->item[i].key);
    if(ift->item[i].value!=NULL) free(ift->item[i].value);
  }
  free(ift->item); ift->item=NULL; ift->_memNr=ift->keyNr=0;
  ift->data=NULL; ift->datasize=0;
}

Referenced by cptReadOne(), defRead(), dftRead(), ehdrEmpty(), ematEmpty(), hrrtMakeCalHdr(), iftdup(), iftRead(), iftReadValue(), imgEmpty(), imgReadMicropetFrame(), imgReadMicropetHeader(), interfileIsHeader(), irdRead(), niftiHeaderToIFT(), readEcat931Calibrationfile(), res2ift(), tsvRead(), vrdRead(), and xelRead().

◆ iftFindNthKey()

int iftFindNthKey	(	IFT *	ift,
		char *	str,
		int	n,
		int	verbose )

extern

Find the Nth item of IFT where the specified string is found in the key. Comparison is case sensitive.

Returns: -1 if key was not found, or other negative value in case of an error, and the index [0..keyNr-1] if matching key is found.

See also: iftFindNthValue, iftGetNth, iftGet, iftGetFrom, iftGetKeyNr, iftRead, iftReplaceNthValue

Parameters

ift	Pointer to IFT structure.
str	Pointer to the case-sensitive (partial) key string.
n	Nth (1..keyNr-1) incidence of value is searched.
verbose	Verbose level; if zero, then only warnings are printed into stderr

Definition at line 84 of file iftsrch.c.

  {
  int li, found_nr=0;
 
  if(verbose>0) printf("%s(*ift, \"%s\", %d)\n", __func__, str, n);
  if(ift==NULL) {return(-10);}
  if(str==NULL || strlen(str)<1) {iftSetStatus(ift, IFT_FAULT); return(-11);}
  if(n<1) {iftSetStatus(ift, IFT_FAULT); return(-11);}
  
  /* Search the list */
  for(li=0; li<ift->keyNr; li++) {
    if(strstr(ift->item[li].key, str)!=NULL) {
      found_nr++;
      if(n==found_nr) {iftSetStatus(ift, IFT_OK); return(li);}
    }
  }
  iftSetStatus(ift, IFT_KEYNOTFOUND);
  return(-1);
}

Referenced by cptReadOne().

◆ iftFindNthValue()

int iftFindNthValue	(	IFT *	ift,
		char *	str,
		int	n,
		int	verbose )

extern

Find the Nth item of IFT where the specified string is found in the value. Comparison is case sensitive.

Returns: -1 if key was not found, or other negative value in case of an error, and the index [0..keyNr-1] if matching value is found.

See also: iftFindNthKey, iftGetFrom, iftGetKeyNr, iftRead, iftReplaceNthValue

Parameters

ift	Pointer to IFT structure.
str	Pointer to the case-sensitive (partial) value string.
n	Nth (1..keyNr-1) incidence of value is searched..
verbose	Verbose level; if zero, then only warnings are printed into stderr.

Definition at line 120 of file iftsrch.c.

  {
  int li, found_nr=0;
 
  if(verbose>0) printf("%s(*ift, \"%s\", %d)\n", __func__, str, n);
  if(ift==NULL) {return(-10);}
  if(str==NULL || strlen(str)<1) {iftSetStatus(ift, IFT_FAULT); return(-11);}
  if(n<1) {iftSetStatus(ift, IFT_FAULT); return(-11);}
  
  /* Search the list */
  for(li=0; li<ift->keyNr; li++) {
    if(strstr(ift->item[li].value, str)!=NULL) {
      found_nr++;
      if(n==found_nr) {iftSetStatus(ift, IFT_OK); return(li);}
    }
  }
  iftSetStatus(ift, IFT_VALUENOTFOUND);
  return(-1);
}

Referenced by cptReadOne(), and xelRead().

◆ iftGet()

int iftGet	(	IFT *	ift,
		char *	key,
		int	verbose )

extern

Find the key in the IFT and return the index [0..keyNr-1]. Key is case insensitive.

See also: iftGetNth, iftFindNthKey, iftGetFrom, iftGetKeyNr, iftRead, iftReplaceNthValue

Returns: -1 if key was not found, or other negative value in case of an error

Parameters

ift	Pointer to IFT structure.
key	Pointer to the key string; contents are replaced by the correct key string.
verbose	Verbose level; if zero, then only warnings are printed into stderr

Definition at line 15 of file iftsrch.c.

  {
  int li;
 
  if(verbose>0) printf("%s(*ift, \"%s\")\n", __func__, key);
  if(ift==NULL) {return(-10);}
  if(key==NULL || strlen(key)<1) {iftSetStatus(ift, IFT_FAULT); return(-11);}
  
  /* Search the list */
  for(li=0; li<ift->keyNr; li++) {
    if(strcasecmp(ift->item[li].key, key)==0) {
      strcpy(key, ift->item[li].key);
      iftSetStatus(ift, IFT_OK); return(li);
    }
  }
  iftSetStatus(ift, IFT_KEYNOTFOUND);
  return(-1);
}

Referenced by cptReadOne(), dft_fill_hdr_from_IFT(), ecat7CopyHeadersNoQuant(), imgGetMicropetHeader(), irdRead(), niftiHeaderFromIFT(), tsvRead(), vrdRead(), and xelRead().

◆ iftGetDoubleValue()

int iftGetDoubleValue	(	IFT *	ift,
		int	si,
		const char *	key,
		double *	value,
		int	verbose )

extern

Finds the specified key string from IFT structure, and reads the corresponding value as double.

Returns: Returns the index of key/value, -1 if key or value was not found, and <-1 in case of an error.

See also: iftGetFloatValue, iftGetIntValue, iftRead

Parameters

ift	Pointer to existing IFT.
si	Index [0..keyNr-1] from which the search is started.
key	Pointer to the key string; search is case-insensitive.
value	Pointer to double variable where value is written; NaN is written in case of an error.
verbose	Verbose level; if zero, then only warnings are printed into stderr.

Definition at line 268 of file iftsrch.c.

  {
  int li;
 
  if(verbose>0) printf("%s(*ift, \"%s\", *value)\n", __func__, key);
  if(ift==NULL) {return(-10);}
  if(value==NULL) {iftSetStatus(ift, IFT_FAULT); return(-10);}
  *value=nan("");
  if(key==NULL || strlen(key)<1) {iftSetStatus(ift, IFT_FAULT); return(-11);}
  if(si<0) {iftSetStatus(ift, IFT_FAULT); return(-12);}
  iftSetStatus(ift, IFT_VALUENOTFOUND);
  /* Search the list */
  for(li=si; li<ift->keyNr; li++) {
    if(strcasecmp(ift->item[li].key, key)==0) {
      if(ift->item[li].value==NULL || strlen(ift->item[li].value)<1) return -1;
      (void)sscanf(ift->item[li].value, "%lf", value);
      if(isnan(*value)) return -1;
      iftSetStatus(ift, IFT_OK); return(li);
    }
  }
  iftSetStatus(ift, IFT_KEYNOTFOUND);
  return(-1);
}

Referenced by hrrtMakeCalHdr().

◆ iftGetFloatValue()

int iftGetFloatValue	(	IFT *	ift,
		int	si,
		const char *	key,
		float *	value,
		int	verbose )

extern

Finds the specified key string from IFT structure, and reads the corresponding value as float.

See also: iftGetDoubleValue, iftGetIntValue

Returns: Returns the index of key/value, -1 if key or value was not found, and <-1 in case of an error.

Parameters

ift	Pointer to existing IFT.
si	Index [0..keyNr-1] from which the search is started.
key	Pointer to the key string; search is case-insensitive.
value	Pointer to float variable where value is written; NaN is written in case of an error.
verbose	Verbose level; if zero, then only warnings are printed into stderr.

Definition at line 228 of file iftsrch.c.

  {
  int li;
 
  if(verbose>0) printf("%s(*ift, \"%s\", *value)\n", __func__, key);
  if(ift==NULL) {return(-10);}
  if(value==NULL) {iftSetStatus(ift, IFT_FAULT); return(-10);}
  *value=nanf("");
  if(key==NULL || strlen(key)<1) {iftSetStatus(ift, IFT_FAULT); return(-11);}
  if(si<0) {iftSetStatus(ift, IFT_FAULT); return(-12);}
  iftSetStatus(ift, IFT_VALUENOTFOUND);
  /* Search the list */
  for(li=si; li<ift->keyNr; li++) {
    if(strcasecmp(ift->item[li].key, key)==0) {
      if(ift->item[li].value==NULL || strlen(ift->item[li].value)<1) return -1;
      (void)sscanf(ift->item[li].value, "%f", value);
      if(isnan(*value)) return -1;
      iftSetStatus(ift, IFT_OK); return(li);
    }
  }
  iftSetStatus(ift, IFT_KEYNOTFOUND);
  return(-1);
}

Referenced by imgGetMicropetHeader(), imgReadMicropetFrame(), niftiHeaderFromIFT(), and upetReadImagedata().

◆ iftGetFrom()

int iftGetFrom	(	IFT *	ift,
		int	si,
		const char *	key,
		int	verbose )

extern

Finds the specified key in the IFT, starting from specified index. Key is case insensitive.

Returns: Returns the index of key/value, -1 if key or value was not found, and <-1 in case of an error.

See also: iftGetFullmatchFrom, iftFindNthValue

Parameters

ift	Pointer to existing IFT.
si	Index [0..keyNr-1] from which the search is started.
key	Pointer to the key string; search is case-insensitive.
verbose	Verbose level; if zero, then only warnings are printed into stderr.

Definition at line 156 of file iftsrch.c.

  {
  int li;
 
  if(verbose>0) printf("%s(*ift, %d, \"%s\")\n", __func__, si, key);
  if(ift==NULL) {return(-10);}
  if(key==NULL || strlen(key)<1) {iftSetStatus(ift, IFT_FAULT); return(-11);}
  if(si<0) {iftSetStatus(ift, IFT_FAULT); return(-12);}
  
  /* Search the list */
  for(li=si; li<ift->keyNr; li++) {
    if(strcasecmp(ift->item[li].key, key)==0) {
      iftSetStatus(ift, IFT_OK); return(li);
    }
  }
  iftSetStatus(ift, IFT_KEYNOTFOUND);
  return(-1);
}

Referenced by hrrtMakeCalHdr(), and interfileIsHeader().

◆ iftGetFullmatchFrom()

int iftGetFullmatchFrom	(	IFT *	ift,
		int	si,
		const char *	key,
		const char *	value,
		int	verbose )

extern

Finds the index with specified key and value in the IFT, starting from specified index. Key and value are case insensitive.

Returns: Returns the index of key/value, -1 if key or value was not found, and <-1 in case of an error.

See also: iftGetFrom, iftFindNthValue

Parameters

ift	Pointer to existing IFT.
si	Index [0..keyNr-1] from which the search is started.
key	Pointer to the key string; search is case-insensitive.
value	Pointer to the value string; search is case-insensitive.
verbose	Verbose level; if zero, then only warnings are printed into stderr.

Definition at line 191 of file iftsrch.c.

  {
  int li;
 
  if(verbose>0)  printf("%s(*ift, %d, \"%s\", \"%s\")\n", __func__, si, key, value);
  if(ift==NULL) {return(-10);}
  if(key==NULL) {iftSetStatus(ift, IFT_FAULT); return(-11);}
  if(value==NULL) {iftSetStatus(ift, IFT_FAULT); return(-12);}
  if(si<0) {iftSetStatus(ift, IFT_FAULT); return(-13);}
  
  /* Search the list */
  for(li=si; li<ift->keyNr; li++) {
    if(strcasecmp(ift->item[li].key, key)!=0) continue;
    if(strcasecmp(ift->item[li].value, value)!=0) continue;
    iftSetStatus(ift, IFT_OK); return(li);
  }
  iftSetStatus(ift, IFT_KEYNOTFOUND);
  return(-1);
}

Referenced by imgReadMicropetFrame(), and upetReadImagedata().

◆ iftGetIntValue()

int iftGetIntValue	(	IFT *	ift,
		int	si,
		const char *	key,
		int *	value,
		int	verbose )

extern

Finds the specified key string from IFT structure, and reads the corresponding value as int.

Returns: Returns the index of key/value, -1 if key or value was not found, and <-1 in case of an error.

See also: iftGetKeyNr, iftGetDoubleValue

Parameters

ift	Pointer to existing IFT.
si	Index [0..keyNr-1] from which the search is started.
key	Pointer to the key string; search is case-insensitive.
value	Pointer to int variable where value is written; -9999 is written in case of an error.
verbose	Verbose level; if zero, then only warnings are printed into stderr.

Definition at line 309 of file iftsrch.c.

  {
  int li;
 
  if(verbose>0) printf("%s(*ift, \"%s\", *value)\n", __func__, key);
  if(ift==NULL) {return(-10);}
  if(value==NULL) {iftSetStatus(ift, IFT_FAULT); return(-10);}
  *value=-9999;
  if(key==NULL || strlen(key)<1) {iftSetStatus(ift, IFT_FAULT); return(-11);}
  if(si<0) {iftSetStatus(ift, IFT_FAULT); return(-12);}
  /* Search the list */
  iftSetStatus(ift, IFT_VALUENOTFOUND);
  for(li=si; li<ift->keyNr; li++) {
    if(strcasecmp(ift->item[li].key, key)==0) {
      if(ift->item[li].value==NULL || strlen(ift->item[li].value)<1) return -1;
      (void)sscanf(ift->item[li].value, "%d", value);
      if(*value==-9999) return -1;
      iftSetStatus(ift, IFT_OK); return(li);
    }
  }
  iftSetStatus(ift, IFT_KEYNOTFOUND);
  return(-1);
}

Referenced by hrrtMakeCalHdr(), imgGetMicropetHeader(), niftiHeaderFromIFT(), and upetReadImagedata().

◆ iftGetKeyNr()

int iftGetKeyNr	(	IFT *	ift,
		const char *	key,
		int	verbose )

extern

Find the nr of occurrences of the specified key in the IFT. Key is case insensitive.

See also: iftFindNthValue

Returns: Returns the nr of occurrences of the key.

Parameters

ift	Pointer to existing IFT
key	Pointer to the key string
verbose	Verbose level; if zero, then only warnings are printed into stderr

Definition at line 350 of file iftsrch.c.

  {
  int li, found_nr=0;
 
  if(verbose>0) printf("%s(*ift, \"%s\")\n", __func__, key);
  if(ift==NULL) {return(0);}
  if(key==NULL || strlen(key)<1) {iftSetStatus(ift, IFT_FAULT); return(0);}
  iftSetStatus(ift, IFT_KEYNOTFOUND);
  
  /* Search the list */
  for(li=0; li<ift->keyNr; li++) if(strcasecmp(ift->item[li].key, key)==0) found_nr++;
  if(found_nr>0) iftSetStatus(ift, IFT_OK);
  else iftSetStatus(ift, IFT_KEYNOTFOUND);
  return(found_nr);
}

◆ iftGetNth()

int iftGetNth	(	IFT *	ift,
		char *	key,
		int	n,
		int	verbose )

extern

Find the Nth key with similar name in the IFT and return the index [0..keyNr-1]. Key is case insensitive.

See also: iftFindNthKey, iftGet, iftGetFrom, iftGetKeyNr, iftRead, iftReplaceNthValue

Returns: -1 if key was not found, or other negative value in case of an error

Parameters

ift	pointer to IFT structure.
key	Pointer to the key string; contents are replaced by the correct key string.
n	Nth (1..) incidence of key is searched.
verbose	Verbose level; if zero, then only warnings are printed into stderr.

Definition at line 48 of file iftsrch.c.

  {
  int li, found_nr=0;
 
  if(verbose>0) printf("%s(*ift, \"%s\", %d)\n", __func__, key, n);
  if(ift==NULL) {return(-10);}
  if(key==NULL || strlen(key)<1) {iftSetStatus(ift, IFT_FAULT); return(-11);}
  if(n<1) {iftSetStatus(ift, IFT_FAULT); return(-11);}
  
  /* Search the list */
  for(li=0; li<ift->keyNr; li++) {
    if(strcasecmp(ift->item[li].key, key)==0) {
      strcpy(key, ift->item[li].key); found_nr++;
      if(n==found_nr) {iftSetStatus(ift, IFT_OK); return(li);}
    }
  }
  iftSetStatus(ift, IFT_KEYNOTFOUND);
  return(-1);
}

Referenced by dft_fill_hdr_from_IFT(), and tsvRead().

◆ iftInit()

void iftInit ( IFT * ift )

extern

Initiate IFT structure. This should be called once before first use.

See also: iftEmpty, iftRead, iftPut

Parameters

ift	Pointer to target IFT structure.

Definition at line 45 of file ift.c.

  {
  if(ift==NULL) return;
  memset(ift, 0, sizeof(IFT));
  ift->_memNr=ift->keyNr=0; ift->item=NULL;
  ift->data=NULL; ift->datasize=0;
}

Referenced by cptReadOne(), dftRead(), ehdrInitiate(), ematInitiate(), hrrtMakeCalHdr(), iftReadValue(), imgInit(), interfileIsHeader(), irdRead(), readEcat931Calibrationfile(), tsvRead(), vrdRead(), and xelRead().

◆ iftPut()

int iftPut	(	IFT *	ift,
		char *	key,
		char *	value,
		char *	cmt_type,
		int	verbose )

extern

Add specified key and its value to the IFT.

Also comment type (first character pointed to) can be added. Either key or value can be empty, but not both of them.

Returns: Returns 0 if ok. Sets ift->status.

See also: iftInit, iftEmpty, iftPutDouble, iftDeleteItem

Parameters

ift	Pointer to initiated IFT; previous contents are not changed
key	Key string; can be empty ("" or NULL)
value	Value string; can be empty ("" or NULL)
cmt_type	Pointer to comment character, e.g. '#' or ';' or '!'; can be empty ("" or NULL)
verbose	Verbose level; if zero, then only warnings are printed into stderr

Definition at line 82 of file ift.c.

  {
  int i, add_nr=1000;
 
  if(ift==NULL) {return(1);}
  if((key==NULL || strlen(key)<1) && (value==NULL || strlen(value)<1)) {
    iftSetStatus(ift, IFT_FAULT); return(2);}
  if(verbose>0) {
    printf("%s(ift, ", __func__);
    if(key!=NULL) printf("\"%s\", ", key); else printf("NULL, ");
    if(value!=NULL) printf("\"%s\", ", value); else printf("NULL, ");
    if(cmt_type!=NULL) printf("\"%1.1s\")\n", cmt_type); else printf("NULL)\n");
  }
 
  /* If necessary, allocate more memory for items */
  if(ift->_memNr<=ift->keyNr) {
    ift->item=realloc(ift->item, (ift->_memNr+add_nr)*sizeof(IFT_KEY_AND_VALUE));
    if(ift->item==NULL) {iftSetStatus(ift, IFT_NOMEMORY); return(5);}
    ift->_memNr+=add_nr;
    for(i=ift->keyNr; i<ift->_memNr; i++) {
      ift->item[i].type=(char)0;
      ift->item[i].sw=(short int)0;
      ift->item[i].key=NULL;
      ift->item[i].value=NULL;
    }
  }
  
  /* Set the contents */
  /* type */
  if(cmt_type!=NULL) ift->item[ift->keyNr].type=cmt_type[0];
  /* key */
  if(key!=NULL) ift->item[ift->keyNr].key=strdup(key);
  else ift->item[ift->keyNr].key=strdup("");
  if(ift->item[ift->keyNr].key==NULL) {iftSetStatus(ift, IFT_NOMEMORY); return(7);}
  /* value */
  if(value!=NULL) ift->item[ift->keyNr].value=strdup(value);
  else ift->item[ift->keyNr].value=strdup("");
  if(ift->item[ift->keyNr].value==NULL) {iftSetStatus(ift, IFT_NOMEMORY); return(8);}
 
  ift->keyNr++;
  iftSetStatus(ift, IFT_OK);
  return(0);
}

Referenced by defRead(), ecat7ImageheaderToIFT(), ecat7MHeaderToIFT(), ecat7ScanheaderToIFT(), hrrtMakeCalHdr(), iftdup(), iftPutDouble(), iftRead(), niftiHeaderToIFT(), and res2ift().

◆ iftPutDouble()

int iftPutDouble	(	IFT *	ift,
		char *	key,
		double	value,
		char *	cmt_type,
		int	verbose )

extern

Add specified key and its floating point (double) value to the IFT.

Also comment type (first character pointed to) can be added. Key can be empty.

Returns: Returns 0 if ok. Sets ift->status.

See also: iftPut

Parameters

ift	Pointer to initiated IFT; previous contents are not changed
key	Key string; can be empty ("" or NULL)
value	Value as double
cmt_type	Pointer to comment character, e.g. '#' or ';' or '!'; can be empty ("" or NULL)
verbose	Verbose level; if zero, then only warnings are printed into stderr

Definition at line 145 of file ift.c.

  {
  char dstr[128];
  sprintf(dstr, "%g", value);
  return(iftPut(ift, key, dstr, cmt_type, verbose));
}

Referenced by hrrtMakeCalHdr().

◆ iftRead()

int iftRead	(	IFT *	ift,
		char *	filename,
		int	is_key_required,
		int	verbose )

extern

Read IFT file keys and values. Previous contents of IFT are preserved.

This function can read the initial ASCII part of files that contain also binary data.

Returns: Returns 0 if ok. Sets ift->status.

See also: iftdup, iftInit, iftEmpty, defRead

Parameters

ift	Pointer to initiated but empty IFT structure.
filename	Input file name.
is_key_required	0=key name is not required, 1=only lines with key and equals sign are read.
verbose	Verbose level; if zero, then only warnings are printed into stderr

Definition at line 24 of file iftfile.c.

  {
  int i, ret, nr=0, line=0, eq_type=0, initial_key_nr=0, nonprintable=0;
  char *allfile, *cptr, *key_ptr, *value_ptr, *eq_ptr, *eq_ptr2, *cmt_ptr;
  char empty_char=(char)0;
  FILE *fp;
 
 
  /* Check function input */
  if(verbose>0) printf("%s(*ift, %s)\n", __func__, filename);
  if(ift==NULL) return(1);
  if(filename==NULL || strlen(filename)<1) {
    iftSetStatus(ift, IFT_FAULT); return(1);
  }
  if(ift->keyNr>0) initial_key_nr=ift->keyNr;
 
  /* Open file */
  if(strcasecmp(filename, "stdin")==0) {
    fp=stdin;
  } else {
    fp=fopen(filename, "r");
    if(fp==NULL) {iftSetStatus(ift, IFT_CANNOTREAD); return(2);}
  }
 
  /* Get file size */
  nr=nonprintable=0; while((ret=fgetc(fp))!=EOF) {
    if(iscntrl(ret) && ret!=13 && ret!=10 && ret!=9) {
      nonprintable=1; break;}
    nr++;
  }
  if(nr<2) {
    if(strcasecmp(filename, "stdin")!=0) fclose(fp);
    if(nonprintable>0) {
      /* File contains non-printable characters; maybe binary file */
      iftSetStatus(ift, IFT_UNKNOWNFORMAT);
    } else {
      /* File just din't have any content */
      iftSetStatus(ift, IFT_NODATA);
    }
    return(3);
  }
  if(verbose>1) printf("  the size of file is %d bytes\n", nr);
  if(nr>5000000) {
    if(strcasecmp(filename, "stdin")!=0) fclose(fp);
    iftSetStatus(ift, IFT_UNKNOWNFORMAT); return(3);
  }
  rewind(fp);
 
  /* Allocate memory for file contents; calloc and extra space to end with 0. */
  allfile=(char*)calloc((nr+3), sizeof(char));
  if(allfile==NULL) {
    if(strcasecmp(filename, "stdin")!=0) fclose(fp);
    iftSetStatus(ift, IFT_NOMEMORY); return(4);
  }
 
  /* Read file contents and close the file */
  i=0; while((ret=fgetc(fp))!=EOF && i<nr) allfile[i++]=(char)ret;
  if(strcasecmp(filename, "stdin")!=0) fclose(fp);
  /* Make sure that file ends with CR or LF */
  if(allfile[i-1]!=10 && allfile[i-1]!=13) {allfile[i++]=(char)13; allfile[i++]=(char)10;}
 
  /* and then fill the list */
  /* separate the first line */
  cptr=strtok(allfile, "\n\r"); line=0;
  do {
    if(verbose>2) printf("line %d: '%s'\n", line, cptr);
    /* Remove initial spaces and tabs */
    i=strspn(cptr, " \t"); cptr+=i;
    if(strlen(cptr)<1) {cptr=strtok(NULL, "\n\r"); line++; continue;}
    /* Check if line starts with a comment character */
    if((cmt_ptr=strchr("#!;%", cptr[0]))!=NULL) {
      cmt_ptr=cptr; cptr++; i=strspn(cptr, " \t"); cptr+=i;
      if(strlen(cptr)<1) {cptr=strtok(NULL, "\n\r"); line++; continue;}
    }
    if(verbose>2) printf("  line %d: '%s'\n", line, cptr);
    /* Find the 'equals' sign */
    eq_ptr=strstr_noquotation(cptr, ":=");
    if(eq_ptr==NULL) eq_ptr=strstr_noquotation(cptr, "=");
    if(eq_ptr==NULL) {
      eq_ptr=strstr_noquotation(cptr, ":");
      /* do not accept time representation */
      if(eq_ptr!=NULL && strlen(cptr)>=strlen(eq_ptr)+2 && istime(eq_ptr-2)<=0) {
        /* ... but search for later equals sign */
        eq_ptr2=strstr_noquotation(eq_ptr+4, ":"); eq_ptr=NULL;
        if(eq_ptr2!=NULL && strlen(cptr)>=strlen(eq_ptr2)+2) {
          if(istime(eq_ptr2-2)<=0) eq_ptr2=NULL; else eq_ptr=eq_ptr2;
        }
      }
    }
    if(eq_ptr==NULL) {
      /* Equals sign not found; if required, then ignore this line */
      if(is_key_required) {cptr=strtok(NULL, "\n\r"); line++; continue;}
      /* If not required, then key="" */
      key_ptr=eq_ptr=&empty_char; value_ptr=cptr;
    } else {
      if(strncmp(eq_ptr, ":=", 2)==0) eq_type=1;
      else if(strncmp(eq_ptr, "=", 1)==0) eq_type=2;
      else if(strncmp(eq_ptr, ":", 1)==0) eq_type=3;
      else eq_type=0;
      *eq_ptr=(char)0; eq_ptr++; key_ptr=cptr;
      /* Find the end of the 'equals' sign; that is the start of value */
      i=strspn(eq_ptr, ":="); value_ptr=eq_ptr+i;
      /* Remove initial spaces and tabs */
      i=strspn(value_ptr, " \t"); value_ptr+=i;
    }
    /* Remove tail spaces and tabs */
    i=strlen(key_ptr); while(i>0 && isspace((int)key_ptr[i-1])) i--; key_ptr[i]=(char)0;
    if(i==0) { /* Length of key name is zero */
      if(is_key_required) {cptr=strtok(NULL, "\n\r"); line++; continue;}
    }
    i=strlen(value_ptr); while(i>0 && isspace((int)value_ptr[i-1])) i--; value_ptr[i]=(char)0;
    if(verbose>2) printf("  key='%s' value='%s'\n", key_ptr, value_ptr);
    /* Remove quotation marks */
    i=strlen(key_ptr)-1; if(i<0) i=0;
    if((key_ptr[0]=='\'' && key_ptr[i]=='\'') || (key_ptr[0]=='\"' && key_ptr[i]=='\"')) {
      key_ptr[i]=(char)0; if(i>0) key_ptr++;}
    if(strlen(key_ptr)<1) { /* Length of key name without comments is zero */
      if(is_key_required) {cptr=strtok(NULL, "\n\r"); line++; continue;}
    }
    i=strlen(value_ptr)-1; if(i<0) i=0;
    if((value_ptr[0]=='\'' && value_ptr[i]=='\'') || (value_ptr[0]=='\"' && value_ptr[i]=='\"')) {
      value_ptr[i]=(char)0; if(i>0) value_ptr++;}
    if(verbose>2) printf("    key='%s' value='%s'\n", key_ptr, value_ptr);
    /* Put key and value in the list */
    ret=iftPut(ift, key_ptr, value_ptr, cmt_ptr, verbose-1);
    if(ret) {
      free(allfile); iftEmpty(ift); iftSetStatus(ift, IFT_FAULT);
      return(10+ret);
    }
    /* separate the next line */
    cptr=strtok(NULL, "\n\r"); line++;
  } while(cptr!=NULL);
  free(allfile);
  if(verbose>2) printf("eq_type=%d\n", eq_type);
  ift->type=eq_type;
  /* Did we actually get any data? */
  if(ift->keyNr<=initial_key_nr) {iftSetStatus(ift, IFT_NODATA); return(7);}
 
  iftSetStatus(ift, IFT_OK);
  return(0);
}

Referenced by cptReadOne(), dftRead(), hrrtMakeCalHdr(), iftReadValue(), interfileIsHeader(), irdRead(), readEcat931Calibrationfile(), tsvRead(), vrdRead(), and xelRead().

◆ iftReadValue()

char * iftReadValue	(	char *	filename,
		char *	keystr,
		int	verbose )

extern

Read value string from IFT file.

Returns: Returns pointer to a locally allocated copy of the value string, or NULL if none found; free the returned pointer when no more needed.

Parameters

filename	File name.
keystr	String to search for in the key. If NULL, or key containing the string is not found, then if there is only one value in the file, pointer to that is returned.
verbose	Verbose level; if zero, then only warnings are printed into stderr

Definition at line 180 of file iftfile.c.

  {
  if(filename==NULL || !filename[0]) return(NULL);
  /* Read the file */
  IFT ift; iftInit(&ift); if(iftRead(&ift, filename, 0, verbose)!=0) return(NULL);
  /* If file contains just one value, then return pointer to a copy of that */
  if(ift.keyNr==1 && strlen(ift.item[0].value)>0) {
    char *s=strdup(ift.item[0].value);
    iftEmpty(&ift);
    return(s);
  }
  /* If key string to search for was not given, then we cannot do more */
  if(keystr==NULL || strlen(keystr)<1) return(NULL);
  /* Search the list for the key */
  int i=-1;
  for(int li=0; li<ift.keyNr; li++) {
    if(strcasestr(ift.item[li].key, keystr)!=NULL && strlen(ift.item[li].value)>0) {
      i=li; break;
    }
  }
  if(i<0) return(NULL);
  /* If found, then return pointer to a copy of that */
  char *s=strdup(ift.item[i].value);
  iftEmpty(&ift);
  return(s);
}

◆ iftReplaceNthValue()

int iftReplaceNthValue	(	IFT *	ift,
		int	item,
		char *	value,
		int	verbose )

extern

Replaces specified value in IFT with a new value.

Returns: 0 if ok.

See also: iftDeleteItem, iftPut

Parameters

ift	Pointer to IFT structure.
item	Index [0..keyNr-1] of key and value.
value	Value string; can be empty ("" or NULL).
verbose	Verbose level; if zero, then only warnings are printed into stderr

Definition at line 206 of file ift.c.

  {
  if(ift==NULL) {return(1);}
  if(item>=ift->keyNr) {iftSetStatus(ift, IFT_FAULT); return(2);}
  if(verbose>0) printf("%s(ift, %d, %s)\n", __func__, item, value);
  /* Delete old value */
  if(ift->item[item].value!=NULL) free(ift->item[item].value);
  /* Set new value */
  if(value!=NULL) ift->item[item].value=strdup(value);
  else ift->item[item].value=strdup("");
  if(ift->item[item].value==NULL) {iftSetStatus(ift, IFT_NOMEMORY); return(8);}
  iftSetStatus(ift, IFT_OK);
  return(0);
}

◆ iftSetStatus()

void iftSetStatus	(	IFT *	ift,
		int	status )

extern

Sets IFT status

Parameters

ift	Pointer to target IFT structure.
status	New status value.

Definition at line 29 of file ift.c.

  {
  if(ift==NULL) return;
  if(status<0 || status>8) ift->status=ift_status[IFT_FAULT];
  else ift->status=ift_status[status];
}

Referenced by defRead(), iftDeleteItem(), iftFindNthKey(), iftFindNthValue(), iftGet(), iftGetDoubleValue(), iftGetFloatValue(), iftGetFrom(), iftGetFullmatchFrom(), iftGetIntValue(), iftGetKeyNr(), iftGetNth(), iftPut(), iftRead(), iftReplaceNthValue(), iftWrite(), and iftWriteItem().

◆ iftWrite()

int iftWrite	(	IFT *	ift,
		char *	filename,
		int	verbose )

extern

Write all keys and values.

Returns: 0 if ok.

See also: iftInit, iftRead, iftWriteItem

Parameters

ift	Pointer to IFT structure.
filename	Output filename, or string "stdout" to print on console.
verbose	Verbose level; if zero, then only warnings are printed into stderr

Definition at line 282 of file iftfile.c.

  {
  int li, ret;
  FILE *fp;
 
  if(verbose>0) printf("%s(*ift, %s)\n", __func__, filename);
  if(ift==NULL) return(1);
  if(filename==NULL || strlen(filename)<1) {iftSetStatus(ift, IFT_FAULT); return(1);}
  if(ift->keyNr<1) return(0);
 
  /* Open file */
  if(strcasecmp(filename, "stdout")==0) {
    fp=stdout;
  } else {
    fp=fopen(filename, "w");
    if(fp==NULL) {iftSetStatus(ift, IFT_CANNOTWRITE); return(2);}
  }
 
  /* Write the contents */
  for(li=0, ret=0; li<ift->keyNr; li++) {
    ret=iftWriteItem(ift, li, fp, verbose);
    if(ret) break;
  }
  if(strcasecmp(filename, "stdout")!=0) fclose(fp);
  return(ret);
}

Referenced by cptReadOne(), ecat7MainheaderFromIFT(), ecat7ReadHeaders(), hrrtMakeCalHdr(), and niftiHeaderFromIFT().

◆ iftWriteItem()

int iftWriteItem	(	IFT *	ift,
		int	item,
		FILE *	fp,
		int	verbose )

extern

Write one item in IFT to the specified file pointer.

See also: iftInit, iftEmpty, iftRead, iftWrite

Returns: 0 if ok.

Parameters

ift	Pointer to initiated but empty IFT structure.
item	Index [0..keyNr-1] of key and value to write.
fp	Output file pointer.
verbose	Verbose level; if zero, then only warnings are printed into stderr

Definition at line 221 of file iftfile.c.

  {
  char eq_sign[3];
  int ret;
 
  if(verbose>0) printf("%s(*ift, %d, fp)\n", __func__, item);
  if(ift==NULL) {return(1);}
  iftSetStatus(ift, IFT_FAULT);
  if(fp==NULL) {return(2);}
  if(item<0 || item>=ift->keyNr) {return(3);}
 
  iftSetStatus(ift, IFT_OK);
  switch(ift->type) {
    case 1: strcpy(eq_sign, ":="); break;
    case 2: strcpy(eq_sign, "="); break;
    case 3: strcpy(eq_sign, ":"); break;
    case 4: strcpy(eq_sign, " "); break;
    case 5: strcpy(eq_sign, "\t"); break;
    case 6: strcpy(eq_sign, ","); break;
    case 7: strcpy(eq_sign, ";"); break;
    default: strcpy(eq_sign, ":="); break;
  }
  if(ift->item[item].type!=' ' && ift->item[item].type!=(char)0) {
    ret=fprintf(fp, "%c ", ift->item[item].type);
    if(ret<1) {iftSetStatus(ift, IFT_CANNOTWRITE); return(6);}
  }
  if(ift->item[item].key==NULL || strlen(ift->item[item].key)<1) {
    ret=fprintf(fp, "%s\n", ift->item[item].value);
  } else {
    if((IFT_SPACE_BEFORE_EQ_SIGN==0 && IFT_SPACE_AFTER_EQ_SIGN==0) ||
       ift->type==4 || ift->type==5 || ift->type==6 || ift->type==7)
      ret=fprintf(fp, "%s%s%s\n",
        ift->item[item].key, eq_sign, ift->item[item].value);
    else if(IFT_SPACE_BEFORE_EQ_SIGN==1 && IFT_SPACE_AFTER_EQ_SIGN==0)
      ret=fprintf(fp, "%s %s%s\n",
        ift->item[item].key, eq_sign, ift->item[item].value);
    else if(IFT_SPACE_BEFORE_EQ_SIGN==0 && IFT_SPACE_AFTER_EQ_SIGN==1)
      ret=fprintf(fp, "%s%s %s\n",
        ift->item[item].key, eq_sign, ift->item[item].value);
    else
      ret=fprintf(fp, "%s %s %s\n",
        ift->item[item].key, eq_sign, ift->item[item].value);
  }
  if(ret<1) {iftSetStatus(ift, IFT_CANNOTWRITE); return(6);}
  return(0);
}

Referenced by iftWrite().

◆ integerListAdd()

int integerListAdd	(	INTEGER_LIST *	l,
		int	v,
		int	ifnew )

extern

Add one integer to INTEGER_LIST.

See also: integerListInit, integerListSort, integerListAddFromString, integerListExpandFromString.

Returns: Returns the number of added integers (0 or 1), or <0 in case of an error.

Parameters

l	Pointer to initiated list
v	Integer value to add
ifnew	Add integer to the list only if it is new (0=no, 1=yes)

Definition at line 190 of file intex.c.

  {
  int i;
  if(l==NULL) return(-1);
 
  /* If only new value is to be added, then check that this is new */
  if(ifnew!=0) {for(i=0; i<l->nr; i++) if(l->list[i]==v) return(0);}
  /* Add value to list if there is space left, and quit */
  if(l->_allocNr>l->nr) {l->list[l->nr++]=v; return(1);}
  /* Allocate more space */
  if(l->_allocNr==0) l->list=(int*)malloc(10*sizeof(int));
  else l->list=(int*)realloc(l->list, (10+l->_allocNr)*sizeof(int));
  if(l->list==NULL) {l->nr=l->_allocNr=0; return(-2);}
  l->_allocNr+=10;
  /* Add value to list */
  l->list[l->nr++]=v;
  return(1);
}

Referenced by imgMaskRoiNr(), integerListAddFromString(), and integerListExpandFromString().

◆ integerListAddFromString()

int integerListAddFromString	(	const char *	s1,
		const char *	s2,
		INTEGER_LIST *	l,
		const int	ifnew )

extern

Read a list of integer values from given string with given delimiters.

See also: integerListInit, integerListSort, integerListAdd, integerListExpandFromString.

Returns: The number of added integer values, or <0 in case of an error.

Author: Vesa Oikonen

Parameters

s1	Pointer to string from which the integers are read, for example "2,3,6,8".
s2	String containing character delimiters, for example ", ".
l	Pointer to INTEGER_LIST struct; previous contents are preserved.
ifnew	Add integer to the list only if it is new (0=no, 1=yes)

Definition at line 242 of file intex.c.

  {
  if(l==NULL) return(-1);
  if(s1==NULL || s2==NULL) return(0);
  
  /* Get the nr of tokens */
  int i, j, m, n, v;
  n=strTokenNr((char*)s1, s2); if(n<1) return(0);
  /* Read the values */
  char tmp[128];
  for(i=j=0; i<n; i++) {
    if(strTokenNCpy(s1, s2, 1+i, tmp, 128)<1) return(-2);
    if(atoi_with_check(tmp, &v)) return(-3);
    m=integerListAdd(l, v, ifnew); if(m<0) return(-4);
    j+=m;
  }
  return(j);
}

◆ integerListEmpty()

int integerListEmpty ( INTEGER_LIST * l )

extern

Free the memory allocated in the INTEGER_LIST struct.

See also: integerListInit, integerListSort, integerListAdd.

Returns: Returns <>0 in case of an error.

Parameters

l	Pointer to INTEGER_LIST struct

Definition at line 175 of file intex.c.

  {
  if(l==NULL) return -1;
  if(l->_allocNr>0) free(l->list);
  l->nr=l->_allocNr=0;
  l->list=NULL;
  return 0;
}

Referenced by imgMaskRoiNr().

◆ integerListExpandFromString()

int integerListExpandFromString	(	const char *	s1,
		const char *	s2,
		INTEGER_LIST *	l,
		const int	ifnew )

extern

Read ranges and individual integer values from given string with given delimiters.

See also: integerListInit, integerListSort, integerListAdd, integerListAddFromString.

Returns: The number of added integer values, or <0 in case of an error.

Author: Vesa Oikonen

Parameters

s1	Pointer to string from which the integers are read, for example "0-8,12,32-28" or "0..8, 12, 28..34".
s2	String containing character delimiters, for example ", ".
l	Pointer to INTEGER_LIST struct; previous contents are preserved.
ifnew	Add integer to the list only if it is new (0=no, 1=yes)

Definition at line 278 of file intex.c.

  {
  if(l==NULL) return(-1);
  if(s1==NULL || s2==NULL) return(0);
  
  /* Get the nr of tokens */
  int n=strTokenNr((char*)s1, s2); if(n<1) return(0);
  /* Read the values */
  char tmp[128], tmp2[128], *t, *tail;
  int i, j, m, first, last, sw;
  for(i=j=0; i<n; i++) {
    if(strTokenNCpy(s1, s2, 1+i, tmp, 128)<1) return(-2);
    t=tmp; errno=0;
    first=strtol(t, &tail, 10); if(errno) return(-3);
    if(*tail) {
      strcpy(tmp2, tail); t=tmp2;
      if(*t=='-') t++;
      else if(*t=='.') {t++; if(*t=='.') t++; else return(-4);}
      else return(-5);
      if(!*t) return(-6);
      last=strtol(t, &tail, 10); if(errno) return(-7);
      if(*tail) return(-8);
    } else {
      last=first;
    }
 
    if(first>last) {sw=first; first=last; last=sw;}
    for(int v=first; v<=last; v++) {
      m=integerListAdd(l, v, ifnew); if(m<0) return(-10);
      j+=m;
    }
  }
  return(j);
}

◆ integerListInit()

int integerListInit ( INTEGER_LIST * l )

extern

Call this (once) before using INTEGER_LIST struct for the first time.

See also: integerListEmpty, integerListSort, integerListAddFromString, integerListExpandFromString.

Returns: Returns <>0 in case of an error.

Parameters

l	Pointer to INTEGER_LIST struct

Definition at line 161 of file intex.c.

  {
  if(l==NULL) return -1;
  l->nr=l->_allocNr=0;
  l->list=NULL;
  return 0;
}

◆ integerListSort()

int integerListSort ( INTEGER_LIST * l )

extern

Sort INTEGER_LIST

See also: integerListInit, integerListAdd, integerListAddFromString, integerListExpandFromString.

Returns: Returns <>0 in case of an error.

Parameters

l	Pointer to INTEGER_LIST struct

Definition at line 219 of file intex.c.

  {
  int i, j, v;
 
  if(l==NULL) return -1;
  if(l->nr<2) return 0;
  for(i=0; i<l->nr; i++) for(j=i+1; j<l->nr; j++) {
    if(l->list[i]>l->list[j]) {
      v=l->list[i]; l->list[i]=l->list[j]; l->list[j]=v;
    }
  }
  return 0;
}

Referenced by imgMaskRoiNr().

◆ intEmpty()

void intEmpty ( INT_list * l )

extern

Free the memory allocated in the INT_LIST struct.

Parameters

l	Pointer to INT_LIST struct

Definition at line 23 of file intex.c.

  {
  if(l==NULL) return;
  if(l->nr>0) free(l->i);
  l->nr=0;
  l->i=NULL;
  return;
}

Referenced by intExpand().

◆ intExpand()

int intExpand	(	char *	text,
		INT_list *	list )

extern

Existing list is freed and all data is cleared. Deprecated. Expanded integers are listed in list.i[] in increasing order.

Parameters

text	Integer expressions to be expanded, e.g. 0-8,12,34-28
list	Pointer for int list data

Returns: 0 if ok and at least one integer is listed.

See also: integerListAddFromString, integerListExpandFromString.

Definition at line 43 of file intex.c.

                                          {
  int j;
  char *p, *t;
  int first, last, swap, intMax=65536;
 
  /* Check the arguments */
  if(strlen(text)<1) return(1);
  intEmpty(list);
 
  /* Expand */
  p=strtok(text, " ,;.&\t\n\r\0");
  while(p!=NULL) {
    t=p; first=last=-1;
    while((*t!='-') && (!isdigit((int)*t)) && (*t)) t++;
    if(*t=='-') {
      while((!isdigit((int)*t)) && (*t)) t++;
      if(isdigit((int)*t)) {first=0; last=atoi(t);}
    } else if(isdigit((int)*t)) {
      first=atoi(t); /*if (first==0) first=1;*/
      while((isdigit((int)*t)) && (*t)) t++;
      if(*t == '-') {
        t++; while((!isdigit((int)*t)) && (*t)) t++;
        if(isdigit((int)*t)) last=atoi(t); else last=intMax;
      }
    }
    if((first>=0) && (last>=0)) {
      if(first>last) {swap=first; first=last; last=swap;}
      if(last>intMax) {if(first<=intMax) last=intMax; else last=0;}
      for(j=first; j<=last && list->nr<intMax; j++)
        if(_intexadd(list, j)<0) return(2);
    } else if(first>=0) {
      if(first<=intMax) if(_intexadd(list, first)<0) return(3);
    } else if(last>=0) {
      if(last>=intMax) if(_intexadd(list, last)<0) return(4);
    }
    p=strtok(NULL, " ,;.&\t\n\r\0");
  }
  if(list->nr<1) return(1);
  return(0);
}

◆ intInit()

void intInit ( INT_list * l )

extern

Call this (once) before using INT_LIST struct for the first time.

Parameters

l	Pointer to INT_LIST struct

Definition at line 12 of file intex.c.

  {
  if(l==NULL) return;
  l->nr=0;
  l->i=NULL;
  return;
}

◆ intMerge()

INT_list intMerge	(	INT_list *	list1,
		INT_list *	list2 )

extern

Merges two lists and returns the result. (the originals are not touched) Duplicate entries are removed. Deprecated.

Parameters

list1	The first list
list2	The second list

Returns: pointer to the new combined list.

Definition at line 124 of file intex.c.

                                                    {
  int r,count=0,l1=0,l2=0;
  INT_list newlist;
  int *tmplist;
  int found;
  tmplist=(int*)malloc(sizeof(int)*(list1->nr+list2->nr));
  while(l1<list1->nr || l2<list2->nr) {
    if(l1<list1->nr) {
      tmplist[count]=list1->i[l1];
      count++;
      l1++;
    }
    found=0;
    if(l2<list2->nr) {
      for(r=0;r<count;r++) {
        if(tmplist[r]==list2->i[l2]) found++;
      }
      if(found<1) {
        tmplist[count]=list2->i[l2]; count++;
      }
      l2++;
    }
  }
  newlist.i=(int*)malloc(sizeof(int)*count);
  memcpy(newlist.i,tmplist,count*sizeof(int));
  newlist.nr=count;
  free(tmplist);
  return newlist;
}

◆ inverfc()

double inverfc ( double x )

extern

Inverse complementary error function erfc^{-1}(x).

Based on the code by Acklam PJ, 2010.

Returns: Returns an approximate value of erfc^{-1}(x).

Parameters

x	Parameter for the inverse complementary error function.

Definition at line 205 of file doubleutil.c.

  {
  static const double a[] = {
    -3.969683028665376E+01, 2.209460984245205E+02, -2.759285104469687E+02,
    1.383577518672690E+02, -3.066479806614716E+01,  2.506628277459239
  };
  static const double b[] = {
    -5.447609879822406E+01, 1.615858368580409E+02, -1.556989798598866E+02, 
    6.680131188771972E+01, -1.328068155288572E+01
  };
  static const double c[] = {
    -7.784894002430293E-03, -3.223964580411365E-01, -2.400758277161838, 
    -2.549732539343734, 4.374664141464968, 2.938163982698783
  };
  static const double d[] = {
    7.784695709041462E-03, 3.224671290700398E-01, 2.445134137142996, 3.754408661907416
  };
  double y, e, u;
  x/=2.0;
  if(x<0.02425) {
    double q=sqrt(-2.0 * log(x));
    y = (((((c[0]*q + c[1])*q + c[2])*q + c[3])*q + c[4])*q + c[5])
         / ((((d[0]*q + d[1])*q + d[2])*q + d[3])*q + 1.0);
  } else if(x<=0.97575) {
    double q=x-0.5;
    double r=q*q;
    y = (((((a[0]*r + a[1])*r + a[2])*r + a[3])*r + a[4])*r + a[5])*q
         / (((((b[0]*r + b[1])*r + b[2])*r + b[3])*r + b[4])*r + 1.0);
  } else {
    double q=sqrt(-2.0 * log(1.0-x));
    y = -(((((c[0]*q + c[1])*q + c[2])*q + c[3])*q + c[4])*q + c[5])
          / ((((d[0]*q + d[1])*q + d[2])*q + d[3])*q + 1.0);
  }
  e=0.5*erfc(-y/M_SQRT2)-x;
  u=(e)*M_SQRT2*M_PI*exp(0.5*y*y);
  y-=u/(1.0 + 0.5*y*u);
  return(fabs(y/M_SQRT2));
}

Referenced by imgGaussianFIRFilter().

◆ isdate()

int isdate ( char * str )

extern

Check if specified string contains date in correct international format (YYYY-MM-DD). String must start with date, but any contents after it is ignored.

Returns: Returns 0 if date is in correct format, -1 if format is correct but date is invalid, and otherwise <>0.

Parameters

str	String to be checked; not changed in this routine

Definition at line 146 of file datetime.c.

  {
  int Y, M, D, n;
  if(strlen(str)<10) return 1;
  if(str[4]!='-' || str[7]!='-') return 2;
  if(strncasecmp(str, "YYYY-MM-DD", 10)==0) return -1;
  n=sscanf(str, "%4d-%2d-%2d", &Y, &M, &D); if(n!=3) return 3;
  if(M>12 || D>31) return -1;
  if(Y<0 || M<1 || D<1) return -1;
  return 0;
}

Referenced by dcmDA2intl(), get_date(), and isdatetime().

◆ isdate2()

int isdate2	(	char *	str,
		char *	intdate )

extern

Check if specified string contains date in correct format (DD.MM.YYYY or DD/MM/YYYY). String must start with date, but any contents after it is ignored.

Returns: Returns 0 if date is in correct format, -1 if format is correct but date is invalid, and otherwise <>0.

Parameters

str	String to be checked; not changed in this routine
intdate	Pointer to allocated string where date is written in international format; enter NULL, if not needed.

Definition at line 168 of file datetime.c.

  {
  int Y, M, D, n;
  if(strlen(str)<10) return 1;
  if(str[2]!='.' && str[2]!='/') return 2;
  if(str[5]!='.' && str[5]!='/') return 2;
  if(isspace(str[8])) return 2;
  if(strncasecmp(str, "DD.MM.YYYY", 10)==0 ||
     strncasecmp(str, "DD/MM/YYYY", 10)==0) return -1;
  n=sscanf(str, "%2d/%2d/%4d", &D, &M, &Y);
  if(n<3) n=sscanf(str, "%2d.%2d.%4d", &D, &M, &Y);
  if(n!=3) return 3;
  if(M>12 || D>31) return -1;
  if(Y<0 || M<1 || D<1) return -1;
  if(intdate!=NULL) sprintf(intdate, "%04d-%02d-%02d", Y, M, D);
  return 0;
}

Referenced by get_date(), and isdatetime().

◆ isdate3()

int isdate3	(	char *	str,
		char *	intdate )

extern

Check if specified string contains date in correct format (DD.MM.YY or DD/MM/YY). String must start with date, but any contents after it is ignored.

Returns: Returns 0 if date is in correct format, -1 if format is correct but date is invalid, and otherwise <>0.

Parameters

str	String to be checked; not changed in this routine
intdate	Pointer to allocated string where date is written in international format; enter NULL, if not needed.

Definition at line 199 of file datetime.c.

  {
  int Y, M, D, n;
  if(strlen(str)<8) return 1;
  if(str[2]!='.' && str[2]!='/') return 2;
  if(str[5]!='.' && str[5]!='/') return 2;
  if(strncasecmp(str, "DD.MM.YY", 8)==0 ||
     strncasecmp(str, "DD/MM/YY", 8)==0) return -1;
  n=sscanf(str, "%2d/%2d/%2d", &D, &M, &Y);
  if(n<3) n=sscanf(str, "%2d.%2d.%2d", &D, &M, &Y);
  if(n!=3) return 3;
  if(Y>99 || M>12 || D>31) return -1;
  if(Y<0 || M<1 || D<1) return -1;
  if(intdate!=NULL) {
    if(Y>=70) Y+=1900; else Y+=2000;
    sprintf(intdate, "%04d-%02d-%02d", Y, M, D);
  }
  return 0;
}

Referenced by get_date(), and isdatetime().

◆ isdate4()

int isdate4	(	int	dateint,
		int *	year,
		int *	month,
		int *	day )

extern

Check if specified integer contains date in format YYYYMMDD.

Returns: Returns 0 if date is in correct format, -1 if format is correct but date is invalid, and otherwise <>0.

Parameters

dateint	Integer to be checked; not changed in this routine
year	Year is written in this pointer; enter NULL if not needed
month	Month is written in this pointer; enter NULL if not needed
day	Day is written in this pointer; enter NULL if not needed

Definition at line 230 of file datetime.c.

  {
  int Y, M, D, n;
  if(dateint<1201 || dateint>99991231) return 1;
  n=dateint/100; D=dateint-100*n; Y=n/100; M=n-100*Y;
  if(M>12 || D>31) return -1;
  if(Y<1 || M<1 || D<1) return -1;
  if(year!=NULL) *year=Y;
  if(month!=NULL) *month=M;
  if(day!=NULL) *day=D;
  return 0;
}

Referenced by ecat7ReadMainheader().

◆ isdatetime()

int isdatetime	(	char *	str,
		char *	intdate )

extern

Check if specified string contains date and time in correct format (YYYY-MM-DD hh:mm:ss, DD.MM.YYYY hh:mm:ss, or DD.MM.YY hh:mm:ss). String must start with date, but any contents after time is ignored.

Returns: Returns 0 if date and time are in correct format, <0 if format is correct but date or time is invalid, and otherwise <>0.

Parameters

str	String to be checked; not changed in this routine
intdate	Pointer to allocated string where date and time is written in international format (YYYY-MM-DD hh:mm:ss); enter NULL, if not needed.

Definition at line 280 of file datetime.c.

  {
  int ret1, ret2, len, type=0;
  char correct_date[20], *time_ptr=NULL;
 
  len=strlen(str); if(len<17) return 1;
  if((ret1=isdate(str))<=0) {
    type=1; strncpy(correct_date, str, 10); correct_date[10]=(char)0; 
  } else if((ret1=isdate2(str, correct_date))<=0) type=2;
  else if((ret1=isdate3(str, correct_date))<=0) type=3;
  else {return 2;}
  if(type==1 || type==2) {
    if(str[10]!=' ' && str[10]!='\t') return 3;
    time_ptr=str+11;
  } else if(type==3) {
    if(str[8]!=' ' && str[8]!='\t') return 3;
    time_ptr=str+9;
  }
  ret2=istime(time_ptr);
  if(ret1>0) {
    if(ret2>=0) return(100*ret1+ret2); else return(100*ret2);
  }
  if(ret2>0) return(10*ret2);
  if(ret1<0) {if(ret2<0) return(-3); else return(-1);}
  if(ret2<0) return(-2);
  if(intdate!=NULL) {
    sprintf(intdate, "%10.10s %8.8s", correct_date, time_ptr);
  }
  return 0;
}

Referenced by dcmDT2intl(), and get_datetime().

◆ isleapyear()

int isleapyear ( long int year )

extern

Check for leap year.

Returns: Nonzero if 'year' is a leap year.

Parameters

year Year

Definition at line 423 of file datetime.c.

  {
  return (year) % 4 == 0 && ((year) % 100 != 0 || (year) % 400 == 0);
}

Referenced by time_to_tm().

◆ istime()

int istime ( char * str )

extern

Check if specified string contains time in correct format (hh:mm:ss). String must start with time, but any contents after it is ignored.

Returns: Returns 0 if time is in correct format, -1 if format is correct but time is invalid, and otherwise <>0.

Parameters

str	String to be checked; not changed in this routine

Definition at line 259 of file datetime.c.

  {
  int h, m, s, n;
  if(strlen(str)<8) return 1;
  if(str[2]!=':' || str[5]!=':') return 2;
  if(strncasecmp(str, "hh:mm:ss", 8)==0) return -1;
  n=sscanf(str, "%d:%d:%d", &h, &m, &s); if(n!=3) return 3;
  if(h<0 || h>23 || m<0 || m>59 || s<0 || s>59) return -1;
  return 0;
}

Referenced by dcmTM2intl(), iftRead(), and isdatetime().

◆ leaps_between()

long int leaps_between	(	long int	year1,
		long int	year2 )

extern

Calculates the number of leap years between year1 and year2.

Returns: Nr of leap years.

Parameters

year1	Year 1
year2	Year 2

Definition at line 435 of file datetime.c.

  {
  long int leaps1, leaps2;
 
  leaps1 = math_div(year1-1, 4) - math_div(year1-1, 100) + math_div(year1-1, 400);
  leaps2 = math_div(year2-1, 4) - math_div(year2-1, 100) + math_div(year2-1, 400);
  return(leaps2-leaps1);
}

Referenced by time_to_tm().

◆ libpet_idcrypt_version()

const char * libpet_idcrypt_version ( void )

extern

Return idcrypt module version info

Definition at line 15 of file idcrypt.c.

                                         {
    return "2004-12-14";
}

◆ little_endian()

int little_endian ( )

extern

Check whether current platform uses little endian byte order. See H&S Sec. 6.1.2 pp. 163-4.

Returns: Returns 1, if current platform is little endian, and 0 if not.

Definition at line 14 of file swap.c.

{
  int x=1;
  if(*(char *)&x==1) return(1); else return(0);
}

Referenced by anaReadHeader(), anaReadImagedata(), anaWriteHeader(), dcmFileWrite(), dcmitemGetInt(), dcmitemGetReal(), dcmReadFileTag(), dcmReadFileVL(), dcmReadFileVRVL(), dcmValueString(), dcmWriteFileTag(), dcmWriteFileVRVL(), ecat63Create(), ecat63Matenter(), ecat63ReadAttnheader(), ecat63ReadImageheader(), ecat63ReadMainheader(), ecat63ReadMatdata(), ecat63ReadMatlist(), ecat63ReadNormheader(), ecat63ReadScanheader(), ecat63WriteAttnheader(), ecat63WriteImageheader(), ecat63WriteMainheader(), ecat63WriteMatdata(), ecat63WriteNormheader(), ecat63WriteScanheader(), ecat7Create(), ecat7EnterMatrix(), ecat7Read2DNormheader(), ecat7Read2DScanheader(), ecat7ReadAttenheader(), ecat7ReadImageheader(), ecat7ReadMainheader(), ecat7ReadMatlist(), ecat7ReadMatrixdata(), ecat7ReadNormheader(), ecat7ReadPolmapheader(), ecat7ReadScanheader(), ecat7Write2DNormheader(), ecat7Write2DScanheader(), ecat7WriteAttenheader(), ecat7WriteImageheader(), ecat7WriteMainheader(), ecat7WriteMatrixdata(), ecat7WriteNormheader(), ecat7WritePolmapheader(), ecat7WriteScanheader(), hrrtMakeCalHdr(), imgSetNiftiHeader(), imgWriteAnalyze(), imgWriteAnalyzeFrame(), niftiReadHeader(), niftiReadImagedata(), niftiWriteHeader(), tiffWriteImg(), and upetReadImagedata().

◆ localtime_r()

struct tm * localtime_r	(	const time_t *	t,
		struct tm *	tm )

extern

Convert time_t to local time in struct tm.

This version of localtime_r function is here for systems (at least Windows) where it is not defined. Uses localtime, which is threadsafe in Windows.

See also: gmtime_r

Returns: Pointer to struct tm, or null in case of an error.

Parameters

t	Pointer to time_t; do not give pointer to int here, like &e7mhdr.scan_start_time
tm	Pointer to struct tm, to be filled here

Definition at line 46 of file datetime.c.

  {
  struct tm *ltm=localtime(t);
  if(ltm==NULL || tm==NULL) return(NULL);
  *tm=*ltm; tm->tm_isdst=-1;
  return tm;
}

◆ math_div()

long int math_div	(	long int	a,
		long int	b )

extern

Division for long integers.

Returns: a/b

Parameters

a	A/b
b	a/B

Definition at line 409 of file datetime.c.

  {
  return a/b - (a%b < 0);
}

Referenced by leaps_between(), and time_to_tm().

◆ petCunit()

char * petCunit ( int cunit )

extern

Return pointer to string describing the calibration data units.

See also: petCunitId, petTunitId, petTunit, cunitFromFilename

Parameters

cunit index of PET_data units_string[].

Definition at line 211 of file petunits.c.

  {
  int n=0;
  while(cunit_string[n]!=0) n++;
  if(cunit<0 || cunit>n-1) return(cunit_string[CUNIT_UNKNOWN]);
  else return(cunit_string[cunit]);
}

Referenced by csv2dft_linkset(), dftGetPmodTitle(), and imgUnit().

◆ petCunitId()

int petCunitId ( const char * unit )

extern

Identify the specified units string as PET data unit.

Returns: Returns the unit id number.

See also: petTunitId, petCunit, petTunit, cunitFromFilename

Definition at line 74 of file petunits.c.

{
  if(unit==NULL) return CUNIT_UNKNOWN;
  if(strlen(unit)==0)                             return CUNIT_UNKNOWN;
  else if(strcasecmp(unit, "unknown")==0)         return CUNIT_UNKNOWN;
  else if(strcasecmp(unit, "cnts/sec")==0)        return CUNIT_CPS;
  else if(strcasecmp(unit, "counts/sec")==0)      return CUNIT_CPS;
  else if(strcasecmp(unit, "ECAT counts/sec")==0) return CUNIT_CPS;
  else if(strcasecmp(unit, "cps")==0)             return CUNIT_CPS;
  else if(strcasecmp(unit, "counts")==0)          return CUNIT_COUNTS;
  else if(strcasecmp(unit, "cnts")==0)            return CUNIT_COUNTS;
  else if(strcasecmp(unit, "kBq/cc")==0)          return CUNIT_KBQ_PER_ML;
  else if(strcasecmp(unit, "kBqcc")==0)           return CUNIT_KBQ_PER_ML;
  else if(strcasecmp(unit, "kBq/mL")==0)          return CUNIT_KBQ_PER_ML;
  else if(strcasecmp(unit, "kBqmL")==0)           return CUNIT_KBQ_PER_ML;
  else if(strcasecmp(unit, "sec*kBq/cc")==0)      return CUNIT_SEC_KBQ_PER_ML;
  else if(strcasecmp(unit, "sec*kBq/mL")==0)      return CUNIT_SEC_KBQ_PER_ML;
  else if(strcasecmp(unit, "integral")==0)        return CUNIT_SEC_KBQ_PER_ML;
  else if(strcasecmp(unit, "1/sec")==0)           return CUNIT_PER_SEC;
  else if(strcasecmp(unit, "1/s")==0)             return CUNIT_PER_SEC;
  else if(strcasecmp(unit, "s-1")==0)             return CUNIT_PER_SEC;
  else if(strcasecmp(unit, "1/min")==0)           return CUNIT_PER_MIN;
  else if(strcasecmp(unit, "min-1")==0)           return CUNIT_PER_MIN;
  else if(strcasecmp(unit, "mL/mL")==0)           return CUNIT_ML_PER_ML;
  else if(strcasecmp(unit, "mL/cc")==0)           return CUNIT_ML_PER_ML;
  else if(strcasecmp(unit, "mL/dL")==0)           return CUNIT_ML_PER_DL;
  else if(strcasecmp(unit, "mL/100mL")==0)        return CUNIT_ML_PER_DL;
  else if(strcasecmp(unit, "mL/(mL*min)")==0)     return CUNIT_ML_PER_ML_PER_MIN;
  else if(strcasecmp(unit, "mL/(min*mL)")==0)     return CUNIT_ML_PER_ML_PER_MIN;
  else if(strcasecmp(unit, "mL/(cc*min)")==0)     return CUNIT_ML_PER_ML_PER_MIN;
  else if(strcasecmp(unit, "mL/(min*cc)")==0)     return CUNIT_ML_PER_ML_PER_MIN;
  else if(strcasecmp(unit, "mL/mL/min")==0)       return CUNIT_ML_PER_ML_PER_MIN;
  else if(strcasecmp(unit, "mL/min/mL")==0)       return CUNIT_ML_PER_ML_PER_MIN;
  else if(strcasecmp(unit, "mL/cc/min")==0)       return CUNIT_ML_PER_ML_PER_MIN;
  else if(strcasecmp(unit, "mL/min/cc")==0)       return CUNIT_ML_PER_ML_PER_MIN;
  else if(strcasecmp(unit, "mL/(dL*min)")==0)     return CUNIT_ML_PER_DL_PER_MIN;
  else if(strcasecmp(unit, "mL/(min*dL)")==0)     return CUNIT_ML_PER_DL_PER_MIN;
  else if(strcasecmp(unit, "mL/(100mL*min)")==0)  return CUNIT_ML_PER_DL_PER_MIN;
  else if(strcasecmp(unit, "mL/(min*100mL)")==0)  return CUNIT_ML_PER_DL_PER_MIN;
  else if(strcasecmp(unit, "mL/dL/min")==0)       return CUNIT_ML_PER_DL_PER_MIN;
  else if(strcasecmp(unit, "mL/min/dL")==0)       return CUNIT_ML_PER_DL_PER_MIN;
  else if(strcasecmp(unit, "mL/100mL/min")==0)    return CUNIT_ML_PER_DL_PER_MIN;
  else if(strcasecmp(unit, "mL/min/100mL")==0)    return CUNIT_ML_PER_DL_PER_MIN;
  else if(strcasecmp(unit, "unitless")==0)        return CUNIT_UNITLESS;
  else if(strcasecmp(unit, "1/1")==0)             return CUNIT_UNITLESS;
  else if(strcasecmp(unit, "Hounsfield Unit")==0) return CUNIT_HU;
  else if(strcasecmp(unit, "HU")==0)              return CUNIT_HU;
  else if(strcasecmp(unit, "nCi/cc")==0)          return CUNIT_NCI_PER_ML;
  else if(strcasecmp(unit, "nCicc")==0)           return CUNIT_NCI_PER_ML;
  else if(strcasecmp(unit, "nCi/mL")==0)          return CUNIT_NCI_PER_ML;
  else if(strcasecmp(unit, "nCimL")==0)           return CUNIT_NCI_PER_ML;
  else if(strcasecmp(unit, "MBq/cc")==0)          return CUNIT_MBQ_PER_ML;
  else if(strcasecmp(unit, "MBqcc")==0)           return CUNIT_MBQ_PER_ML;
  else if(strcasecmp(unit, "MBq/mL")==0)          return CUNIT_MBQ_PER_ML;
  else if(strcasecmp(unit, "MBqmL")==0)           return CUNIT_MBQ_PER_ML;
  else if(strcasecmp(unit, "Bq/cc")==0)           return CUNIT_BQ_PER_ML;
  else if(strcasecmp(unit, "Bqcc")==0)            return CUNIT_BQ_PER_ML;
  else if(strcasecmp(unit, "Bq/mL")==0)           return CUNIT_BQ_PER_ML;
  else if(strcasecmp(unit, "BqmL")==0)            return CUNIT_BQ_PER_ML;
  else if(strcasecmp(unit, "uCi/cc")==0)          return CUNIT_UCI_PER_ML;
  else if(strcasecmp(unit, "uCicc")==0)           return CUNIT_UCI_PER_ML;
  else if(strcasecmp(unit, "uCi/mL")==0)          return CUNIT_UCI_PER_ML;
  else if(strcasecmp(unit, "uCimL")==0)           return CUNIT_UCI_PER_ML;
  else if(strcasecmp(unit, "umol/(100g*min)")==0) return CUNIT_UMOL_PER_MIN_PER_100G;
  else if(strcasecmp(unit, "umol/(min*100g)")==0) return CUNIT_UMOL_PER_MIN_PER_100G;
  else if(strcasecmp(unit, "umol/100g/min")==0)   return CUNIT_UMOL_PER_MIN_PER_100G;
  else if(strcasecmp(unit, "umol/min/100g")==0)   return CUNIT_UMOL_PER_MIN_PER_100G;
  else if(strcasecmp(unit, "umol/(dL*min)")==0)   return CUNIT_UMOL_PER_MIN_PER_100G;
  else if(strcasecmp(unit, "umol/(min*dL)")==0)   return CUNIT_UMOL_PER_MIN_PER_100G;
  else if(strcasecmp(unit, "umol/dL/min")==0)     return CUNIT_UMOL_PER_MIN_PER_100G;
  else if(strcasecmp(unit, "umol/min/dL")==0)     return CUNIT_UMOL_PER_MIN_PER_100G;
  else if(strcasecmp(unit, "mg/(100g*min)")==0)   return CUNIT_MG_PER_MIN_PER_100G;
  else if(strcasecmp(unit, "mg/(min*100g)")==0)   return CUNIT_MG_PER_MIN_PER_100G;
  else if(strcasecmp(unit, "mg/100g/min")==0)     return CUNIT_MG_PER_MIN_PER_100G;
  else if(strcasecmp(unit, "mg/min/100g")==0)     return CUNIT_MG_PER_MIN_PER_100G;
  else if(strcasecmp(unit, "mg/(dL*min)")==0)     return CUNIT_MG_PER_MIN_PER_100G;
  else if(strcasecmp(unit, "mg/(min*dL)")==0)     return CUNIT_MG_PER_MIN_PER_100G;
  else if(strcasecmp(unit, "mg/dL/min")==0)       return CUNIT_MG_PER_MIN_PER_100G;
  else if(strcasecmp(unit, "mg/min/dL")==0)       return CUNIT_MG_PER_MIN_PER_100G;
  else if(strcasecmp(unit, "%")==0)               return CUNIT_PERCENTAGE;
  else if(strcasecmp(unit, "kcps")==0)            return CUNIT_KCPS;
  else if(strcasecmp(unit, "min*kBq/cc")==0)      return CUNIT_MIN_KBQ_PER_ML;
  else if(strcasecmp(unit, "min*kBq/mL")==0)      return CUNIT_MIN_KBQ_PER_ML;
  else if(strcasecmp(unit, "Bq")==0)              return CUNIT_BQ;
  else if(strcasecmp(unit, "kBq")==0)             return CUNIT_KBQ;
  else if(strcasecmp(unit, "MBq")==0)             return CUNIT_MBQ;
  else if(strcasecmp(unit, "GBq")==0)             return CUNIT_GBQ;
  else if(strcasecmp(unit, "nCi")==0)             return CUNIT_NCI;
  else if(strcasecmp(unit, "uCi")==0)             return CUNIT_UCI;
  else if(strcasecmp(unit, "mCi")==0)             return CUNIT_MCI;
  else if(strcasecmp(unit, "%ID")==0)             return CUNIT_PID;
  else if(strcasecmp(unit, "% ID")==0)            return CUNIT_PID;
  else if(strcasecmp(unit, "%ID/g")==0)           return CUNIT_PIDM;
  else if(strcasecmp(unit, "% ID/g")==0)          return CUNIT_PIDM;
  else if(strcasecmp(unit, "%ID/mL")==0)          return CUNIT_PIDV;
  else if(strcasecmp(unit, "% ID/mL")==0)         return CUNIT_PIDV;
  else if(strcasecmp(unit, "%ID/cc")==0)          return CUNIT_PIDV;
  else if(strcasecmp(unit, "% ID/cc")==0)         return CUNIT_PIDV;
  else if(strcasecmp(unit, "g/mL")==0)            return CUNIT_G_PER_ML;
  else if(strcasecmp(unit, "g/cc")==0)            return CUNIT_G_PER_ML;
  else if(strncasecmp(unit, "SUV", 3)==0)         return CUNIT_G_PER_ML;
  else if(strcasecmp(unit, "mL/g")==0)            return CUNIT_ML_PER_G;
  else if(strcasecmp(unit, "cc/g")==0)            return CUNIT_ML_PER_G;
 
  return CUNIT_UNKNOWN;
}

Referenced by csv2dft_a(), csv2dft_linkset(), cunit_check_dft_vs_img(), dftGetPmodTitle(), dftReadinput(), dftReadModelingData(), dftReadReference(), dftSetComments(), dftTimeIntegral(), dftUnitConversion(), and dftWrite().

◆ petTunit()

char * petTunit ( int tunit )

extern

Return pointer to string describing the time unit.

See also: petCunitId, petTunitId, petCunit, cunitFromFilename

Parameters

tunit index of PET_time unit_string[].

Definition at line 226 of file petunits.c.

  {
  int n=0;
  while(tunit_string[n]!=0) n++;
  if(tunit<0 || tunit>n-1) return(tunit_string[TUNIT_UNKNOWN]);
  else return(tunit_string[tunit]);
}

Referenced by check_times_dft_vs_dft(), csv2dft_linkset(), dftInterpolateForIMG(), dftMatchTimeunits(), dftSetComments(), dftWrite(), dftWriteHTML(), fitWrite(), imgReadModelingData(), res_allocate_with_dft(), and tsvRead().

◆ petTunitId()

int petTunitId ( const char * timeunit )

extern

Identifies the specified string as PET time (x axis) units.

Returns: Returns the timeunit id number.

See also: petCunitId, petCunit, petTunit, cunitFromFilename

Definition at line 187 of file petunits.c.

{
  if(timeunit==NULL) return TUNIT_UNKNOWN;
  if(strlen(timeunit)==0)                         return TUNIT_UNKNOWN;
  else if(strcasecmp(timeunit,  "unknown")==0)    return TUNIT_UNKNOWN;
  else if(strncasecmp(timeunit, "seconds", 3)==0) return TUNIT_SEC;
  else if(strcmp(timeunit,      "s")==0)          return TUNIT_SEC;
  else if(strncasecmp(timeunit, "minutes", 3)==0) return TUNIT_MIN;
  else if(strcasecmp(timeunit,  "um")==0)         return TUNIT_UM;
  else if(strcasecmp(timeunit,  "mm")==0)         return TUNIT_MM;
  else if(strcasecmp(timeunit,  "cm")==0)         return TUNIT_CM;
  else if(strcasecmp(timeunit,  "m")==0)          return TUNIT_M;
  else if(strcasecmp(timeunit,  "h")==0)          return TUNIT_HOUR;
  else if(strcasecmp(timeunit,  "months")==0)     return TUNIT_MONTH;
  else if(strcasecmp(timeunit,  "y")==0)          return TUNIT_YEAR;
  else if(strcasecmp(timeunit,  "msec")==0)       return TUNIT_MSEC;
  return TUNIT_UNKNOWN;
}

Referenced by csv2dft_linkset(), dft_fill_hdr_from_IFT(), dftGetPmodTitle(), dftRead(), dftTimeunitToDFT(), and fitRead().

◆ printf32bits()

void printf32bits ( void * buf )

extern

Printfs as bit string the 32-bit variable pointed to by buf. Far from being optimized, thus only for testing and development purposes.

Parameters

buf	Pointer to memory

Definition at line 133 of file swap.c.

                             {
  unsigned int u, i;
  int j;
 
  memcpy(&u, buf, 4);
  for(i=32; i>0; i--) {
    j=i-1;if(i<32 && (i%8)==0) printf(" ");
    if(u & (1L<<j)) printf("1"); else printf("0");
  }
  printf("\n");
}

◆ readStrtokens()

int readStrtokens	(	const char *	filename,
		char ***	toklist )

extern

Read list of string tokens from specified file Remember to free the memory of string list.

Returns: Returns the nr of tokens or <0 in case of an error.

Parameters

filename	Name of file to read
toklist	Pointer to list of strings read and allocated here

Definition at line 181 of file readfile.c.

  {
  int i, ret, nr=0;
  char *allfile, *cptr;
  static char **list;
  FILE *fp;
 
  /* Check arguments */
  if(strlen(filename)<1) return(-1);
  /* Open file */
  fp=fopen(filename, "r"); if(fp==NULL) return(-2);
  /* Get file size */
  nr=0; while((ret=fgetc(fp))!=EOF) nr++; rewind(fp);
  if(nr<1) {fclose(fp); return(0);}
  /* Allocate memory for file contents */
  allfile=(char*)malloc((nr+1)*sizeof(char));
  if(allfile==NULL) {fclose(fp); return(-4);}
  /* Read file contents */
  i=0; while((ret=fgetc(fp))!=EOF && i<nr) allfile[i++]=(char)ret;
  fclose(fp); allfile[i]=(char)0;
  /* Get the number of tokens */
  cptr=allfile; nr=0;
  while(1) {
    i=strspn(cptr, " ;,|\t\n\r"); cptr+=i;
    i=strcspn(cptr, " ;,|\t\n\r"); cptr+=i;
    if(i==0 || *cptr==0) break; else nr++;
  }
  /*printf("Nr of tokens in %s: %d\n", filename, nr);*/
  /* Allocate memory for array of strings */
  list=(char**)malloc(nr*sizeof(char*));
  if(list==NULL) {free(allfile); return(-5);}
  /* and then fill the list */
  cptr=strtok(allfile, " ;,|\t\n\r");
  for(i=0; i<nr; i++) {
    if(cptr==NULL) {
      for(--i; i>=0; i--) free(list[i]);
      free(list); free(allfile);
      return(-8);
    }
    list[i]=(char*)malloc( (strlen(cptr)+1)*sizeof(char) );
    if(list[i]==NULL) {
      for(--i; i>=0; i--) free(list[i]);
      free(list); free(allfile);
      return(-9);
    }
    strcpy(list[i], cptr);
    cptr=strtok(NULL, " ;,|\t\n\r");
  }
  free(allfile);
  *toklist=list;
  return(nr);
}

◆ rnameCatenate()

int rnameCatenate	(	char *	rname,
		int	max_rname_len,
		char *	name1,
		char *	name2,
		char *	name3,
		char	space )

extern

Construct full TAC name from up to three subnames.

Returns: Returns 0 when successful, <>0 if failed.

Parameters

rname	Pointer to string of length max_rname_len, in where the full name will be placed
max_rname_len	Length of full TAC name, not including null char
name1	Pointer to 1st subname (anatomical region); NULL if not available
name2	Pointer to 1st subname (usually hemisphere); NULL if not available
name3	Pointer to 1st subname (usually image plane); NULL if not available
space	Spacing character, for example ' ', '_', or '-'

Definition at line 189 of file rname.c.

  {
  unsigned int len, rlen;
 
  if(rname==NULL || max_rname_len<1) return 1;
  if(name1==NULL && name2==NULL && name3==NULL) return 2;
 
  strcpy(rname, "");
 
  len=strlen(name1);
  if(len>0 && len<(unsigned int)max_rname_len && strcmp(name1, ".")!=0)
    strcpy(rname, name1);
 
  len=strlen(name2); rlen=strlen(rname);
  if(len>0 && (len+strlen(rname)+1)<(unsigned int)max_rname_len 
           && strcmp(name2, ".")!=0) {
    if(rlen>0) {rname[rlen]=(char)space; rname[rlen+1]=(char)0;}
    strcat(rname, name2);
  }
 
  len=strlen(name3); rlen=strlen(rname);
  if(len>0 && (len+strlen(rname)+1)<(unsigned int)max_rname_len 
           && strcmp(name3, ".")!=0) {
    if(rlen>0) {rname[rlen]=(char)space; rname[rlen+1]=(char)0;}
    strcat(rname, name3);
  }
 
  if(strlen(rname)<1) return 3;
  return 0;
}

Referenced by dftRead(), dftRNameSimplify(), fitRead(), and resRead().

◆ rnameMatch()

int rnameMatch	(	char *	rname,
		int	rnr,
		char *	test_str )

extern

Test whether region name or number matches with a test string. Test string can contain wildcards. If test string contains only one subname, it is tested against whole rname. If it contains 2-3 subnames, those are tested against the corresponding tokens in rname. Subname '.' stands for empty name. Number is tested only if test string contains one token of all digits.

Returns: Returns 1, in case of match, or 0 if not matched.

Parameters

rname	Region name which is tested
rnr	Region number (1..)
test_str	Test string

Definition at line 144 of file rname.c.

  {
  char region[3][MAX_REGIONNAME_LEN+1];
  char test[3][MAX_REGIONNAME_LEN+1];
  unsigned int ni, m, rtoknr=0, ttoknr=0;
 
  /* Check the input */
  if(rname==NULL || strlen(rname)<1) return(0);
  if(test_str==NULL || strlen(test_str)<1) return(0);
  /* Split region names into 1-3 subnames */
  rtoknr=rnameSplit(rname, region[0], region[1], region[2], MAX_REGIONNAME_LEN);
  ttoknr=rnameSplit(test_str, test[0], test[1], test[2], MAX_REGIONNAME_LEN);
  if(rtoknr==0 || ttoknr==0) return(0);
  /* If more than one subname to test for, then test against corresponding
     subnames */
  if(ttoknr>1) {
    for(ni=0; ni<ttoknr; ni++) {
      if( strcmp(test[ni], ".") && fncasematch(test[ni], region[ni])==0 ) {
        return(0);}
    }
    return(1);
  }
  /* If just one subname to test for */
  /* If it contains only digits, then we assume that it is region number */
  for(ni=m=0; ni<strlen(test[0]); ni++)
    if(isdigit((int)test[0][ni])==0) {m++; break;}
  if(m==0 && (ni=atoi(test[0]))>0) { /* it indeed is an acceptable  number */
    if((unsigned int)rnr==ni) return(1); else return(0);
  }
  /* It is name; chack if it is found anywhere in the region name */
  for(ni=0; ni<rtoknr; ni++) if(fncasematch(test[0], region[ni])) return(1);
  return(0);
}

Referenced by dftSelectRegions(), and resSelectRegions().

◆ rnameRmDots()

int rnameRmDots	(	char *	rname1,
		char *	rname2 )

extern

Region name may contain dots marking non-existing identification of hemisphere or image plane etc. This function removes the dots and extra space characters.

Returns: Returns <0 in case of an error, otherwise the number of tokens that were left.

Parameters

rname1	String which contains the original region name; the modified string is written in next string, if pointer to it is given next.
rname2	Pointer to previously allocated string, into which the modified region name will be written. Enter NULL, if previous string is to be modified instead.

Definition at line 99 of file rname.c.

  {
  char *temp, *out, *cptr, *lptr;
  int len, c=0;
 
  if(rname1==NULL) return -1;
  len=strlen(rname1);
  if(len<1) {if(rname2!=NULL) strcpy(rname2, ""); return 0;}
  temp=malloc(len+1); if(temp==NULL) return -3;
  strcpy(temp, rname1);
  if(rname2==NULL) out=rname1; else out=rname2;
  strcpy(out, "");
 
  /* remove dots and extra spaces */
  lptr=temp;
  cptr=strtok(lptr, " \t\n\r");
  while(cptr!=NULL) {
    if(strcmp(cptr, ".")!=0) {
      if(strlen(out)>0) strcat(out, " ");
      strcat(out, cptr); c++;
    }
    cptr=strtok(NULL, " \t\n\r");
  }
  free(temp);
  return c;
}

Referenced by plot_fit_svg(), plot_fitrange_svg(), plot_svg(), res2ift(), and rnameSplit().

◆ rnameSplit()

int rnameSplit	(	char *	rname,
		char *	name1,
		char *	name2,
		char *	name3,
		int	max_name_len )

extern

Split region name into 1-3 subparts of given max length.

Returns: Returns the number of subparts.

Parameters

rname	Region name to split (string is not edited)
name1	Pointer to 1st subname (anatomical region)
name2	Pointer to 2nd subname (usually hemisphere)
name3	Pointer to 3rd subname (usually image plane)
max_name_len	Max length of subnames, excluding terminal null

Definition at line 14 of file rname.c.

  {
  char temp[MAX_REGIONNAME_LEN+1], temp2[MAX_REGIONNAME_LEN+1];
  char *cptr, *lptr, space[64];
  int nr;
 
  if(rname==NULL || name1==NULL || name2==NULL || name3==NULL) return(0);
#if(0)
      printf("  rname := '%s'\n", rname);  
#endif
  if(max_name_len<1) return(0);
  name1[0]=name2[0]=name3[0]=(char)0;
  strncpy(temp, rname, MAX_REGIONNAME_LEN); temp[MAX_REGIONNAME_LEN]=(char)0;
 
  /* Try to divide long name with only upper case letters as possible dividers */
  if((int)strlen(rname)>max_name_len) {
    int nr1, nr2;
    strcpy(space, " _-."); nr1=strChrCount(temp, space);
    nr2=strUppercaseCount(temp+1); // not counting first character if that is uppercase
    if(nr1==0 && nr2>0 && nr2<=3) {
      nr=0;
      int len=strlen(temp);
      int i=1, n=0;
      cptr=temp;
      while(i<len && !isupper(temp[i])) i++;
      n=i; if(n>max_name_len) n=max_name_len;
      strncpy(name1, cptr, n); name1[n]=(char)0; nr++;
      cptr=temp+i; n=0;
      while(i<len && (!isupper(temp[i]) || n==0)) {i++; n++;}
      if(n>max_name_len) n=max_name_len;
      if(n>0) {
        strncpy(name2, cptr, n); name2[n]=(char)0; nr++;
        cptr=temp+i; n=0;
        while(i<len && (!isupper(temp[i]) || n==0)) {i++; n++;}
        if(n>max_name_len) n=max_name_len;
        if(n>0) {strncpy(name3, cptr, n); name3[n]=(char)0; nr++;}
      }
#if(0)
      printf("  -> subnames := '%s' '%s' '%s'\n", name1, name2, name3);  
#endif
      return(nr);
    }
  }
 
  nr=rnameRmDots(temp, temp2);
#if(0)  
  printf("temp = '%s'\n", temp);  
  printf("temp2 = '%s'\n", temp2);  
#endif
  strcpy(space, " \t");
  if(nr<3) {
    if(strChrCount(temp2, space)<2) strcat(space, "_");
    if(strChrCount(temp2, space)<2) strcat(space, "-");
  }
  strcat(space, "\n\r");
  nr=0; lptr=temp;
  cptr=strtok(lptr, space); if(cptr==NULL) return(nr);
  strncpy(name1, cptr, max_name_len); name1[max_name_len]=(char)0; nr++;
  cptr=strtok(NULL, space); if(cptr==NULL) return(nr);
  strncpy(name2, cptr, max_name_len); name2[max_name_len]=(char)0; nr++;
  cptr=strtok(NULL, space); if(cptr==NULL) return(nr);
  strncpy(name3, cptr, max_name_len); name3[max_name_len]=(char)0; nr++;
#if(0)
  printf("  -> subnames := '%s' '%s' '%s'\n", name1, name2, name3);  
#endif
  return(nr);
}

Referenced by csv2dft_a(), csv2dft_b(), csv2dft_linkset(), csv2dft_mat(), dft_fill_hdr_from_IFT(), dftGetPmodTitle(), resRead(), and rnameMatch().

◆ roinameExists()

int roinameExists ( char * roiname )

extern

Verifies whether TAC name exists or not.

TAC name string may contain only delimiters like '.', '_', '-', or spaces. Those cases are interpreted as no name in this function.

Returns: 1 if valid name exists, 0 if not.

Parameters

roiname ROI name string; not edited

Definition at line 241 of file rname.c.

  {
  if(roiname==NULL) return(0);
  size_t len, n;
  len=strlen(roiname); if(len==0) return(0);
  n=strspn(roiname, " _-.\t"); if(n==len) return(0);
  return(1);
}

◆ statSortDouble()

void statSortDouble	(	double *	data,
		unsigned int	n,
		int	order )

extern

Sort the given double array into ascending or descending order.

Author: Vesa Oikonen

Parameters

data	Pointer to data array of size n
n	Length of data array
order	Ascending (0) or descending (<>0) order

Definition at line 267 of file doubleutil.c.

  {
  if(n<2 || data==NULL) return;
  if(order==0) qsort(data, n, sizeof(double), statDoubleCompAsc);
  else qsort(data, n, sizeof(float), statDoubleCompDesc);
}

◆ statSortFloat()

void statSortFloat	(	float *	data,
		unsigned int	n,
		int	order )

extern

Sort the given double list into ascending or descending order.

Author: Vesa Oikonen

Sort the given float array into ascending or descending order.

Author: Vesa Oikonen

Parameters

data	Pointer to data array of size n
n	Length of data array
order	Ascending (0) or descending (<>0) order

Definition at line 76 of file imgqntls.c.

  {
  if(n<2 || data==NULL) return;
  if(order==0) qsort(data, n, sizeof(float), statFloatCompAsc);
  else qsort(data, n, sizeof(float), statFloatCompDesc);
}

◆ str_token_list_add()

int str_token_list_add	(	STR_TOKEN_LIST *	lst,
		char *	new_item )

extern

Put a string in STR_TOKEN_LIST.

Returns: Returns 0, if successful, and <>0 in case of an error.

Parameters

lst	List that has to be initialized beforehand.
new_item	String that is added to list.

Definition at line 42 of file readfile.c.

  {
  int i;
  const int add_nr=10;
 
  if(lst==NULL || new_item==NULL || strlen(new_item)<1) return(1);
  if(lst->list_size<=lst->token_nr) {
    lst->tok=realloc(lst->tok, sizeof(char*)*(lst->list_size+add_nr));
    if(lst->tok==NULL) return(2);
    for(i=lst->list_size; i<lst->list_size+add_nr; i++) lst->tok[i]=NULL;
    lst->list_size+=add_nr;
  }
  lst->tok[lst->token_nr]=strdup(new_item);
  if(lst->tok[lst->token_nr]==NULL) return(3);
  lst->token_nr++;
 
  return(0);
}

Referenced by str_token_list_read(), and textfileReadLines().

◆ str_token_list_del()

int str_token_list_del	(	STR_TOKEN_LIST *	lst,
		int	item )

extern

Remove the specified string item from the STR_TOKEN_LIST.

Returns: Returns 0, if successful, and <>0 in case of an error.

Parameters

lst	List that has to be initialized beforehand.
item	Item number to remove (1..item_nr).

Definition at line 70 of file readfile.c.

  {
  //printf("str_token_list_del(list with %d items, %d)\n", lst->token_nr, item);
  if(lst==NULL || item<1 || item>lst->token_nr) return(1);
  if(lst->tok[item-1]!=NULL) free(lst->tok[item-1]);
  for(int i=item; i<lst->token_nr; i++) {
    //printf("  index %d -> %d\n", i, i-1);
    lst->tok[i-1]=lst->tok[i]; lst->tok[i]=NULL;
  }
  lst->token_nr--;
  return(0);
}

Referenced by sifRead().

◆ str_token_list_empty()

void str_token_list_empty ( STR_TOKEN_LIST * lst )

extern

Free memory allocated for STR_TOKEN_LIST. All contents are destroyed.

Parameters

lst	Pointer to list to be emptied.

Definition at line 26 of file readfile.c.

  {
  for(int i=0; i<lst->list_size; i++)
    if(lst->tok[i]!=NULL) free(lst->tok[i]);
  if(lst->tok!=NULL) free(lst->tok);
  lst->tok=NULL;
  lst->list_size=0; lst->token_nr=0;
}

Referenced by roi_read(), sifRead(), str_token_list_read(), and textfileReadLines().

◆ str_token_list_init()

void str_token_list_init ( STR_TOKEN_LIST * lst )

extern

Initiate STR_TOKEN_LIST structure. This should be called once before first use.

Parameters

lst	Pointer to list to be initiated.

Definition at line 13 of file readfile.c.

  {
  memset(lst, 0, sizeof(STR_TOKEN_LIST));
  lst->list_size=0; lst->token_nr=0; lst->tok=NULL;
}

Referenced by roi_read(), and sifRead().

◆ str_token_list_read()

int str_token_list_read	(	const char *	filename,
		STR_TOKEN_LIST *	lst )

extern

Read all string tokens from text file into STR_TOKEN_LIST. List needs to be initialized. Previous contents are deleted.

Returns: Returns 0, if successful, and <>0 in case of an error.

Parameters

filename	Name of text file to read
lst	Token list is allocated by this function.

Definition at line 93 of file readfile.c.

  {
  int i, ret, nr=0;
  char *allfile, *cptr;
  FILE *fp;
 
  if(lst==NULL || filename==NULL || strlen(filename)<1) return(1);
  str_token_list_empty(lst);
 
  /* Open file */
  fp=fopen(filename, "r"); if(fp==NULL) return(2);
  /* Get file size */
  nr=0; while((ret=fgetc(fp))!=EOF) nr++; rewind(fp);
  if(nr<1) {fclose(fp); return(0);}
  /* Allocate memory for file contents */
  allfile=(char*)malloc((nr+1)*sizeof(char));
  if(allfile==NULL) {fclose(fp); return(3);}
  /* Read file contents */
  i=0; while((ret=fgetc(fp))!=EOF && i<nr) allfile[i++]=(char)ret;
  fclose(fp); allfile[i]=(char)0;
  /* and then fill the list */
  cptr=strtok(allfile, " ;,|\t\n\r");
  if(cptr==NULL) {free(allfile); return(4);}
  do {
    if(str_token_list_add(lst, cptr)) {free(allfile); return(10);}
    cptr=strtok(NULL, " ;,|\t\n\r");
  } while(cptr!=NULL);
  free(allfile);
 
  return(0);
}

◆ strcasestr()

char * strcasestr	(	const char *	haystack,
		const char *	needle )

extern

Case-insensitive version of strstr().

Returns: a pointer to the beginning of the first occurrence, or NULL if not found.

Parameters

haystack	Pointer to string in which sub-string needle is searched.
needle	Pointer to sub-string which is searched for in source string haystack.

Definition at line 279 of file strext.c.

  {
  if(!haystack || !*haystack || !needle || !*needle) return 0;
 
  const char *s=haystack, *p=needle;
  do {
    if(!*p) return(char*)haystack;
    if((*p==*s) || (tolower(*p)==tolower(*s))) {
      p++; s++;
    } else {
      p=needle; if(!*s) return(NULL);
      s=++haystack;
    }
  } while(1);
  return *p ? NULL : (char*)haystack;
}

Referenced by cunitFromFilename(), dftFormat(), iftReadValue(), and resWriteHTML_table().

◆ strChrCount()

int strChrCount	(	const char *	str1,
		const char *	str2 )

extern

Count how many times specified characters are found in a string. Search is case-sensitive.

See also: strTokenNr, strReplaceChar, strUppercaseCount

Returns: Returns the nr of characters of str2 found in str1.

Parameters

str1	String to search for characters; not modified.
str2	String containing characters which are searched for; not modified.

Definition at line 126 of file strext.c.

  {
  unsigned int n=0, i, j;
  if(str1==NULL || str2==NULL || strlen(str1)==0 || strlen(str2)==0) return n;
  for(i=0; i<strlen(str1); i++)
    for(j=0; j<strlen(str2); j++)
      if(str1[i]==str2[j]) n++;
  return n;
}

Referenced by rnameSplit().

◆ strCleanForXML()

void strCleanForXML ( char * s )

extern

Remove from string those characters that would require encoding in XML. Replaced by character '-'.

See also: strEncodeForXML

Parameters

s	Pointer to the string to be cleaned.

Definition at line 402 of file strext.c.

  {
  if(s==NULL || *s=='\0') return;
  strReplaceChar(s, '&', '-');
  strReplaceChar(s, '\'', '-');
  strReplaceChar(s, '\"', '-');
  strReplaceChar(s, '<', '-');
  strReplaceChar(s, '>', '-');
  return;
}

◆ strCleanSpaces()

int strCleanSpaces ( char * s )

extern

Removes any initial and trailing space characters from specified string s.

Space characters in the middle of the string are not removed.

Returns: 0 when successful, otherwise >0.

Author: Vesa Oikonen

Parameters

s	Pointer to the string.

Definition at line 343 of file strext.c.

  {
  if(s==NULL) return 0;
  int len=strlen(s); if(len<0) return 0;
  char *s2; s2=strdup(s); if(s2==NULL) return(1);
  int n=strncpyCleanSpaces(s2, s, len+1);
  if(n<1) strcpy(s, ""); else strcpy(s, s2);
  free(s2);
  return 0;
}

Referenced by csvRead().

◆ strEncodeForXML()

char * strEncodeForXML ( const char * s )

extern

Encode special characters for XML, including SVG.

Postcondition: Free the memory of returned string pointer.

Returns: Returns pointer to encoded string. NULL is returned in case of an error, or if encoding is not necessary.

See also: svg_str_encode

Parameters

s	Pointer to the string to be encoded.

Definition at line 364 of file strext.c.

  {
  if(s==NULL) return(NULL);
  /* Count the characters needing encoding */
  int n=0;
  for(size_t i=0; i<strlen(s); i++) {
    if(s[i]=='&') {n++; continue;}
    if(s[i]=='\'') {n++; continue;}
    if(s[i]=='\"') {n++; continue;}
    if(s[i]=='<') {n++; continue;}
    if(s[i]=='>') {n++; continue;}
  }
  if(n==0) return(NULL);
  /* Allocate memory for new string (one char to max 6 chars) */
  n*=5; n+=strlen(s)+1;
  char *ns=(char*)malloc(n*sizeof(char));
  if(ns==NULL) return(NULL);
  /* Process the string */
  for(int i=0; i<n; i++) ns[i]=(char)0;
  for(size_t i=0; i<strlen(s); i++) {
    if(s[i]=='&') {strcat(ns, "&amp;"); continue;}
    if(s[i]=='\'') {strcat(ns, "&apos;"); continue;}
    if(s[i]=='\"') {strcat(ns, "&quot;"); continue;}
    if(s[i]=='<') {strcat(ns, "&lt;"); continue;}
    if(s[i]=='>') {strcat(ns, "&gt;"); continue;}
    ns[strlen(ns)]=s[i];
  }
  return(ns);
}

Referenced by dftWriteHTML(), and resWriteHTML_table().

◆ strlcat()

size_t strlcat	(	char *	dst,
		const char *	src,
		size_t	dstsize )

extern

Safer version of strncat. At most dstsize-1 characters are appended from the source string to destination string. Destination string will be NUL terminated, unless dstsize <= strlen(dst).

Remarks: Included in POSIX but not in GCC.

Returns: the size of the buffer that would have been needed for the destination string; if >=dstsize, then truncation occurred.

Parameters

dst	Destination string.
src	Source string.
dstsize	The actual length of buffer allocated for the destination string; for example, destination string has been allocated as char dst[dstsize];

Definition at line 206 of file strext.c.

  {
  char *d;
  const char *s=src;
  size_t dlen, n;
 
  /* Find the current length of dst */
  dlen=strnlen(dst, dstsize);
  if(s==NULL) return(dlen);
  n=dstsize-dlen;
  if(n==0) return(dlen+strlen(s));
  d=dst+dlen;
  while(*s!='\0') {
    if(n!=1) {*d=*s; d++; n--;}
    s++;
  }
  *d='\0';
  return(dlen+(s-src));
}

Referenced by anaExists(), backupExistingFile(), dftRead(), and hrrtMakeCalHdr().

◆ strlcpy()

size_t strlcpy	(	char *	dst,
		const char *	src,
		size_t	dstsize )

extern

Safer version of strncpy or strcpy.

At most dstsize-1 characters are copied from the source string to the destination string. Destination string will be NUL terminated.

Remarks: Included in POSIX but not in GCC.

Returns: the size of the buffer that would have been needed for the destination string; if >=dstsize, then truncation occurred.

Parameters

dst	Destination string.
src	Source string.
dstsize	The actual length of buffer allocated for the destination string; for example, destination string has been allocated as char dst[dstsize];

Definition at line 245 of file strext.c.

  {
  if(dstsize>0) dst[0]='\0';
  if(strlen(src)==0) return(0);
 
  char *d=dst;
  const char *s=src;
  size_t n;
 
  /* Copy as many chars as allowed */
  n=dstsize;
  if(n!=0) while(--n!=0) {*d=*s; if(*d=='\0') {d++; s++; break;} d++; s++;}
  if(n==0) { // not enough space, add NUL, and check how much space were needed
    if(dstsize!=0) *d='\0';
    while(*s++) {}
  }
  return(s-src-1);
}

Referenced by anaEditHeader(), anaExists(), backupExistingFile(), csv2dft_b(), csv2dft_mat(), dcmAddItem(), dcmFileRead(), dft_fill_hdr_from_IFT(), dftAllocateWithIMG(), dftRead(), dftWrite(), ecat63AddImg(), ecat63EditMHeader(), ecat63ReadAllToImg(), ecat63ReadPlaneToImg(), ecat63WriteAllImg(), ecat7EditMHeader(), ecat7EditVHeader(), ecatCopy63to7mainheader(), ecatCopy7to63mainheader(), fitRead(), hrrtMakeCalHdr(), img2sif(), imgGetAnalyzeHeader(), imgGetEcat63MHeader(), imgGetEcat7MHeader(), imgGetMicropetSIF(), imgReadAnalyze(), imgReadAnalyzeHeader(), imgReadNiftiHeader(), imgSetEcat63MHeader(), imgSetEcat7MHeader(), imgSetNiftiHeader(), imgWriteAnalyze(), niftiCreateFNames(), niftiExists(), niftiHeaderToIFT(), resRead(), selectEcat931Calibrationfile(), sif2img(), strncpyCleanSpaces(), strTokenDup(), strTokenNCpy(), studynr_from_fname2(), studynr_in_fname(), upetExists(), and xelRead().

◆ strnCopyClean()

int strnCopyClean	(	char *	str1,
		const char *	str2,
		int	maxlen )

extern

Copy str2 to str1, removing any quotation marks around the string, and making sure that string fits to str2.

Returns: Returns the length of the new string str2.

Parameters

str1	Pointer to pre-allocated result string with length of at least maxlen characters, including NULL character
str2	Pointer to the original string; not changed in this function
maxlen	Max length of str1, including the end NULL

Definition at line 52 of file quots.c.

  {
  char *cptr;
  int i;
 
  if(str1==NULL || maxlen<1) return(0);
  str1[0]=(char)0; if(str2==NULL || strlen(str2)<1) return(0);
  cptr=(char*)str2; cptr+=strspn(str2, "\"\'\t\n\r ");
  strncpy(str1, cptr, maxlen-1); str1[maxlen-1]=(char)0;
  i=strlen(str1); if(i<1) return(0);
  cptr=str1+(i-1);
  while(i>0) {
    if(*cptr!='\"' && *cptr!='\'' && *cptr!='\t' && *cptr!='\n' &&
       *cptr!='\r' && *cptr!=' ')
           break;
    i--; str1[i]=(char)0; cptr--;
  }
  return(i);
}

Referenced by upetHeaderReadParameter().

◆ strncpyCleanSpaces()

int strncpyCleanSpaces	(	char *	s1,
		const char *	s2,
		int	maxlen )

extern

Version of strncpy() which as usual copies s2 to s1, but without any space characters or line end characters that may be around the string s2.

Returns: the length of the new string s1.

Author: Vesa Oikonen

Parameters

s1	Pointer to pre-allocated result string with length of at least maxlen characters, including NULL character.
s2	Pointer to the original string.
maxlen	Max length of s1, including the trailing zero.

Definition at line 308 of file strext.c.

  {
  if(s1==NULL) return(0);
  s1[0]=(char)0; if(maxlen<1) return(0);
  if(maxlen<2) {strcpy(s1, ""); return(0);}
  if(s2==NULL || strlen(s2)<1) return(0);
 
  char *cptr;
  int i;
 
  cptr=(char*)s2; cptr+=strspn(s2, "\t\n\r ");
  strlcpy(s1, cptr, maxlen); i=strlen(s1); if(i<1) return(0);
  cptr=s1+(i-1);
  while(i>0) {
    if(*cptr!='\t' && *cptr!='\n' && *cptr!='\r' && *cptr!=' ') break;
    i--; s1[i]=(char)0; cptr--;
  }
  return(i);
}

Referenced by roi_read(), roiRead(), and strCleanSpaces().

◆ strnlen()

size_t strnlen	(	const char *	s,
		size_t	n )

extern

Safer version of strlen, in case the argument s is not NUL terminated string. Computes the length of string s, but never scans beyond the n first bytes of the string.

Remarks: Included in POSIX and GCC, so this implementation may not be needed.

Returns: same as strlen() or n, whichever is smaller.

Parameters

s	Pointer to string, or character array, that may not be NULL terminated.
n	The actual length of buffer allocated for the string; for example, string could have been allocated as char s[n];

Definition at line 181 of file strext.c.

  {
  if(s==NULL) return(0);
  char *ps=(char*)s;
  size_t i=0;
  while(i<n && *ps!='\0') {i++; ps++;}
  return(i);
}

Referenced by csv2dft_mat(), dcmAddItem(), dcmDA2intl(), dcmDT2intl(), dcmFileRead(), dcmFileWrite(), dcmSOPIdentify(), dcmSOPUIDName(), dcmTM2intl(), dcmTrUID(), dftWrite(), irdRead(), sifRead(), strlcat(), and tpcHtmlUsage().

◆ strPtrToNextValue()

char * strPtrToNextValue	(	char *	str,
		char **	nxtp )

extern

This function searches the given string for a string representation of numerical value, possibly with decimal and exponent part. Returns also pointer to the string right after where the numerical value ended, and from where the next number can be searched. Return Returns pointer to start of the next string, or NULL if not found.

Parameters

str	Pointer to string where numerical values are searched; not changed.
nxtp	Obligatory pointer to string pointer, which will be set to point to the first character after value string, or to NULL if string ends.

Definition at line 134 of file decpoint.c.

  {
  char *cptr, *strt;
  unsigned int i, j, n=0;
 
  // search the start
  if(str==NULL || strlen(str)<1) {*nxtp=NULL; return NULL;}
  // find the index where number might start
  cptr=str;
  do {
    i=strcspn(cptr, "+-0123456789"); if(i==strlen(cptr)) {*nxtp=NULL; return NULL;}
    strt=cptr+i;
    // any previous character must be a separator
    if(i>0) {cptr=strt-1; i=strcspn(cptr, " \t,;"); if(i>0) cptr+=i;}
  } while(i>0);
  // get past + and - ; if more than one, that is an error
  cptr=strt; if(*cptr=='-' || *cptr=='+') {n++; cptr++;}
  // get past first set of numbers
  i=strspn(cptr, "0123456789"); n+=i;
  if(i<1) {*nxtp=cptr; return NULL;}
  cptr+=i; n+=i;
  // get past decimal separator
  i=strspn(cptr, ",.");
  if(i==0) { // there was none, thus we're done
    *nxtp=cptr; return strt; 
  } else if(i>1) { // can't be more than one decimal separator
    *nxtp=cptr; return strt;
  }
  cptr+=1; n+=1;
  // after decimal separator we must have numbers again
  i=strspn(cptr, "0123456789");
  if(i<1) { // what was thought to be decimal separator was not that
    *nxtp=cptr-1; return strt;
  }
  cptr+=i; n+=i;
  // check if we have exponent part
  if(*cptr=='E' || *cptr=='e') i=1; else i=0;
  // if not, then we're done
  if(i==0) {*nxtp=cptr; return strt;}
  cptr+=1;
  // we may have signed exponent
  if(*cptr=='-' || *cptr=='+') {cptr+=1; i++;}
  // after exponent we must have numbers again
  j=strspn(cptr, "0123456789");
  if(j<1) { // what was thought to be exponent was not that
    *nxtp=cptr-i; return strt;
  }
  cptr+=j; n+=j;
  // now this must be the end
  *nxtp=cptr;
  return strt;
}

◆ strReplaceChar()

void strReplaceChar	(	char *	str,
		char	c1,
		char	c2 )

extern

Replace certain characters in string with another character.

Parameters

str	Pointer to string in which the character is replaced.
c1	Character to be replaced.
c2	Character to use instead. If NULL, then only the first character is replaced.

Definition at line 159 of file strext.c.

  {
  char *cptr;
  if(strlen(str)==0) return;
  while((cptr=strchr(str, c1))!=NULL) *cptr=c2;
  return;
}

Referenced by dftWrite(), fitRead(), roi_read(), roiRead(), and strCleanForXML().

◆ strstr_noquotation()

char * strstr_noquotation	(	const char *	str1,
		const char *	str2 )

extern

The strstr_noquotation() function returns a pointer to the first occurrence in the string pointed to by str1, excluding parts that are inside quotation marks "" or '', of the string pointed to by str2.

Returns: Returns NULL pointer if no match is found, or if found, then pointer to the first occurrence.

Parameters

str1	Pointer to string to be searched
str2	Pointer to string with quotation marks

Definition at line 17 of file quots.c.

  {
  unsigned int i, test_len;
  unsigned int single_quotation=0;
  unsigned int double_quotation=0;
  char *cptr;
  
  if(str1==NULL) return((char*)NULL);
  if(str2==NULL) return((char*)str1);
  test_len=strlen(str2); if(test_len<1) return((char*)str1);
  for(i=0, cptr=(char*)str1; i<strlen(str1); i++, cptr++) {
    if(*cptr=='\'') {
      if(single_quotation==0) single_quotation=1; else single_quotation=0;
      continue;
    }
    if(*cptr=='\"') {
      if(double_quotation==0) double_quotation=1; else double_quotation=0;
      continue;
    }
    if(single_quotation==1 || double_quotation==1) continue;
    if(strncmp(cptr, str2, test_len)==0) return(cptr);
  }
  return((char*)NULL);
}

Referenced by iftRead().

◆ strTokenDup()

char * strTokenDup	(	const char *	s1,
		const char *	s2,
		int *	next )

extern

Search the string s1 for the first token. The characters making up the string s2 are the delimiters that determine the tokens.

Returns: Returns pointer to a copy of the token string, or NULL in case of an error or if no token found.

Postcondition: Remember to free the memory from the returned pointer after last use.

Author: Vesa Oikonen

Parameters

s1	String from where tokens are searched; not modified in any way.
s2	String containing character delimiters.
next	Index of s1 where the token ended; set to NULL, if not needed.

Definition at line 89 of file strext.c.

  {
  if(next!=NULL) *next=0;
  if(s1==NULL) return NULL;
 
  char *s3=NULL, *cptr;
  size_t j;
  
  /* If no delimiters, then return copy of s1 */
  if(s2==NULL || strlen(s2)<1) {
    s3=strdup(s1); if(next!=NULL) *next=strlen(s1);
    return s3;
  } 
  /* Pass initial delimiter characters */
  cptr=(char*)s1; j=strspn(cptr, s2); cptr+=j; if(next!=NULL) *next=j;
  /* calculate characters between delimiters */
  j=strcspn(cptr, s2); if(j==0) {return NULL;}
  if(next!=NULL) *next+=j;
  /* Allocate space for token */
  s3=calloc(j+1, sizeof(char)); if(s3==NULL) return NULL;
  strlcpy(s3, cptr, j+1);
  return s3;
}

Referenced by fitRead(), and textfileReadLines().

◆ strTokenNCpy()

int strTokenNCpy	(	const char *	str1,
		const char *	str2,
		int	i,
		char *	str3,
		int	count )

extern

The strTokenNCpy() function copies the i'th token in the string pointed to by str1 into string pointed to by str3. The characters making up the string pointed to by str2 are the delimiters that determine the token.

See also: strTokenNr, strTokenDup

Returns: Returns the length of token, 0 if no token(s) found.

Parameters

str1	String from where tokens are searched; not modified in any way.
str2	String containing character delimiters.
i	Token number to copy (1..nr of tokens).
str3	String array into where the token is copied; string will be null terminated.
count	Length of str3, including terminal null

Definition at line 45 of file strext.c.

  {
  int j=0, n=0;
  char *cptr;
  if(str1==NULL || str2==NULL || strlen(str1)==0 || strlen(str2)==0) return(0);
  if(i<1 || str3==NULL || count<2) return(0);
 
  cptr=(char*)str1;
  do {
    // pass delimiter characters
    j=strspn(cptr, str2); cptr+=j;
    // pass characters between delimiters
    j=strcspn(cptr, str2); if(j>0) n++;
    // if this is the required token nr, then stop here
    if(n==i) {
      if(j>count-1) j=count-1;
      strlcpy(str3, cptr, j+1); //strncpy(str3, cptr, j); str3[j]=(char)0;
      break;
    }
    cptr+=j;
  } while(j>0);
  if(n>i) {str3[0]=(char)0; return(0);}
  return(j);
}

Referenced by fitRead(), integerListAddFromString(), integerListExpandFromString(), pxlRead(), resRead(), and roi_read().

◆ strTokenNr()

int strTokenNr	(	const char *	str1,
		const char *	str2 )

extern

The strTokenNr() function returns the number of tokens in the string pointed to by str1. The characters making up the string pointed to by str2 are the delimiters that determine the token.

See also: strTokenNCpy, strChrCount

Returns: Returns the nr of tokens.

Parameters

str1	String from where tokens are calculated; not modified in any way.
str2	String containing character delimiters.

Definition at line 17 of file strext.c.

  {
  int i=0, n=0;
  char *cptr;
  if(str1==NULL || str2==NULL || strlen(str1)==0 || strlen(str2)==0) return(0);
 
  cptr=(char*)str1;
  do {
    // pass delimiter characters
    i=strspn(cptr, str2); cptr+=i;
    // pass characters between delimiters
    i=strcspn(cptr, str2); cptr+=i; if(i>0) n++; 
  } while(i>0);
  return(n);
}

Referenced by fitRead(), integerListAddFromString(), integerListExpandFromString(), pxlRead(), resRead(), and roi_read().

◆ strUppercaseCount()

int strUppercaseCount ( const char * s )

extern

Count how many upper case characters are found in string.

Returns: Returns the nr of upper case characters found in s.

See also: strChrCount

Parameters

s	String to search for upper case characters; not modified.

Definition at line 146 of file strext.c.

  {
  unsigned int n=0, i;
  if(s==NULL || strlen(s)==0) return (int)n;
  for(i=0; i<strlen(s); i++) if(isupper(s[i])) n++;
  return (int)n;
}

Referenced by rnameSplit().

◆ studynr_from_fname()

int studynr_from_fname	(	char *	fname,
		char *	studynr )

extern

Extract study number (max MAX_STUDYNR_LEN chars) from filename. This function removes initial zeroes from the number part, and converts uppercase letters to lowercase, if necessary.

Returns: 0 if successful.

Parameters

fname	Filename may include path. Filename is not modified.
studynr	Pointer for the resulting study number. Memory (>=MAX_STUDYNR_LEN+1 chars) for it must be allocated before calling this.

Definition at line 119 of file studynr.c.

  {
  return studynr_from_fname2(fname, studynr, 0);
}

Referenced by cptReadOne(), dftRead(), imgMicropetCTToEcat7(), imgMicropetPETToEcat7(), resFName2study(), resRead(), tsvRead(), and xelRead().

◆ studynr_from_fname2()

int studynr_from_fname2	(	char *	fname,
		char *	studynr,
		int	force )

extern

Extract study number (max MAX_STUDYNR_LEN chars) from filename.

This function removes initial zeroes from the number part, and converts uppercase letters to lowercase, if necessary.

Returns: Returns 0 if successful.

Parameters

fname	Filename, which may include the path. Filename is not modified.
studynr	Pointer for the resulting study number. Memory (>=MAX_STUDYNR_LEN+1 chars) must be allocated before calling this.
force	Validity of studynr is verified (0) or not verified (1)

Definition at line 67 of file studynr.c.

  {
  unsigned int i;
  char *cptr;
 
  //printf("studynr_from_fname2()\n");
 
  if(fname==NULL || studynr==NULL) return(1);
  strcpy(studynr, "");
 
  /* At first try if studynr_in_fname() works */
  if(studynr_in_fname(fname, studynr)==0) return(0);
 
  /* Remove path */
  cptr=strrchr(fname, '/'); if(cptr==NULL) cptr=strrchr(fname, '\\');
  if(cptr==NULL) cptr=fname; else cptr++;
  //i=strlen(cptr); if(i>MAX_STUDYNR_LEN) i=MAX_STUDYNR_LEN;
  strlcpy(studynr, cptr, MAX_STUDYNR_LEN); //studynr[i]=(char)0;
  if(force!=0) return(0);
 
  /* Check that first character is a letter */
  if(!isalpha((int)studynr[0])) {strcpy(studynr, ""); return(4);}
  /* Remove everything after the letter+digit parts */
  for(i=1; i<strlen(studynr); i++) if(!isalpha((int)studynr[i])) break;
  for(; i<strlen(studynr); i++) if(!isdigit((int)studynr[i])) break;
  studynr[i]=(char)0;
  /* Check the length of the study number */
  if(strlen(studynr)<2) {strcpy(studynr, ""); return(5);}
  /* Remove initial zeroes from the number part */
  if(studynr_rm_zeroes(studynr)!=0) {strcpy(studynr, ""); return(6);}
  /* Convert characters to lower case */
  if(studynr_to_lowercase(studynr)!=0) {strcpy(studynr, ""); return(7);}
  /* Check the validity of the study number */
  if(!studynr_validity_check(studynr)) {strcpy(studynr, ""); return(8);}
  return(0);
}

Referenced by imgGetMicropetHeader(), and studynr_from_fname().

◆ studynr_in_fname()

int studynr_in_fname	(	char *	fname,
		char *	studynr )

extern

Find study number (max MAX_STUDYNR_LEN chars) inside filename; Study number must contain 1-5 letters followed by 1-5 digits, if such string is not found then error code is returned.

Initial zeroes are removed if necessary and uppercase letters are changed to lowercase.

Returns: Returns 0 if successful, otherwise <>0.

Parameters

fname	Filename may include path. Filename is not modified.
studynr	Pointer to string (>=MAX_STUDYNR_LEN+1 chars) where the resulting study number is written.

Definition at line 18 of file studynr.c.

  {
  unsigned int i;
  int ret;
  char *cptr, *lptr, temp[FILENAME_MAX], temp2[FILENAME_MAX];
 
  //printf("studynr_in_filename(%s, string)\n", fname);
  if(fname==NULL || studynr==NULL) return(1);
  for(i=0; i<=MAX_STUDYNR_LEN; i++) studynr[i]=(char)0;
  /* Remove path */
  cptr=strrchr(fname, '/'); if(cptr==NULL) cptr=strrchr(fname, '\\');
  if(cptr==NULL) cptr=fname; else cptr++;
  //i=strlen(cptr); if(i>FILENAME_MAX-1) i=FILENAME_MAX-1;
  strlcpy(temp, cptr, FILENAME_MAX); lptr=temp;
  /* Verify tokens in filename whether they are valid as study number */
  cptr=strtok(lptr, "_-+{}!~.()");
  while(cptr!=NULL && !studynr[0]) {
    strcpy(temp2, cptr);
    /* Remove everything after the letter+digit parts */
    for(i=1; i<strlen(temp2); i++) if(!isalpha((int)temp2[i])) break;
    for(; i<strlen(temp2); i++) if(!isdigit((int)temp2[i])) break;
    temp2[i]=(char)0;
    /* Only then check it */
    ret=studynr_validity_check2(temp2, 1);
    if(ret==1) {
      strcpy(studynr, temp2); 
      if(studynr_rm_zeroes(studynr)!=0) strcpy(studynr, "");
      if(studynr_to_lowercase(studynr)!=0) strcpy(studynr, "");
    }
    cptr=strtok(NULL, "_-+{}!~.()");
  }
  if(!studynr[0]) return(2);
  return 0;
}

Referenced by studynr_from_fname2().

◆ studynr_match()

int studynr_match	(	char *	studynr1,
		char *	studynr2 )

extern

Check whether two valid study numbers are the same. If either of study numbers is shorter than the other, the end parts are compared; thus study numbers that are changed by SPM can be matched. Argument strings are not modified.

Parameters

studynr1	compared number
studynr2	compared number

Returns: 1 if study numbers match exactly, 2 if the match is probable, and zero, if no match is found.

Definition at line 142 of file studynr.c.

                                                  {
  int len1, len2;
  char *cptr1, *cptr2;
 
  len1=strlen(studynr1); if(len1<2 || len1>MAX_STUDYNR_LEN) return(0);
  len2=strlen(studynr2); if(len2<2 || len2>MAX_STUDYNR_LEN) return(0);
  if(len1==len2 && strcmp(studynr1, studynr2)==0) return(1);
  if(len2>len1) cptr2=studynr2+len2-len1; else cptr2=studynr2;
  if(len1>len2) cptr1=studynr1+len1-len2; else cptr1=studynr1;
  if(strcasecmp(cptr1, cptr2)==0) return(2);
  return(0);
}

◆ studynr_rm_zeroes()

int studynr_rm_zeroes ( char * studynr )

extern

Remove zeroes from the number part of the PET study number.

Parameters

studynr modified study number.

Returns: nonzero in case of failure.

Definition at line 211 of file studynr.c.

                                     {
  int i, j, len;
 
  len=strlen(studynr); if(len<2) return(1);
  if(isdigit((int)studynr[0])) return(2);
  if(!isdigit((int)studynr[len-1])) return(3);
  for(i=1; i<len; i++) if(isdigit((int)studynr[i])) break;
  for(j=i; j<len; j++) if(studynr[j]!='0') break;
  if(i==j) return(0); /* no initial zeroes */
  for(;j<=len; i++, j++) studynr[i]=studynr[j];
  return(0);
}

Referenced by studynr_from_fname2(), and studynr_in_fname().

◆ studynr_to_lowercase()

int studynr_to_lowercase ( char * studynr )

extern

Convert the PET study number letters to lowercase. Conversion is not done to non-valid study number.

Parameters

studynr modified study number.

Returns: nonzero in case of failure.

Definition at line 233 of file studynr.c.

                                        {
  int i, len;
 
  len=strlen(studynr); if(len<2) return(1);
  if(isdigit((int)studynr[0])) return(2);
  if(!isdigit((int)studynr[len-1])) return(3);
  for(i=0; i<len; i++) {
    if(isdigit((int)studynr[i])) break;
    studynr[i]=(char)tolower((int)studynr[i]);
  }
  return(0);
}

Referenced by studynr_from_fname2(), and studynr_in_fname().

◆ studynr_validity_check()

int studynr_validity_check ( char * studynr )

extern

Check that the argument string is a valid TPC study number.

Valid study number here is defined as containing 1-5 letters (upper- or lowercase) followed by at least 1 digit, with total length of max 10 characters.

Returns: 1 if study number is valid, zero if not valid or in case of failure.

Parameters

studynr Evaluated study number; not modified.

Definition at line 196 of file studynr.c.

  {
  return studynr_validity_check2(studynr, 0);
}

Referenced by ecat63ReadAllToImg(), ecat63ReadPlaneToImg(), imgGetEcat63MHeader(), imgGetEcat7MHeader(), and studynr_from_fname2().

◆ studynr_validity_check2()

int studynr_validity_check2	(	char *	studynr,
		int	zero_ok )

extern

Check that the argument string is a valid TPC study number.

Valid study number here is defined as containing 1-5 letters (upper- or lowercase) followed by at least 1 digit, with total length of max 10 characters.

Returns: 1 if study number is valid, zero if not valid or in case of failure.

Parameters

studynr	String to be evaluated as study number. Not modified here.
zero_ok	Number part of study number may start with zero (1) or may not (0)

Definition at line 166 of file studynr.c.

  {
  int i, j, len;
 
  len=strlen(studynr); if(len<2 || len>MAX_STUDYNR_LEN) return(0);
  for(i=0; i<len; i++) if(!isalnum((int)studynr[i])) return(0);
  for(i=0; i<len; i++) if(!isalpha((int)studynr[i])) break;
  if(i<1 || i>5) return(0);
  if(zero_ok==0 && studynr[i]=='0') return(0); /* first digit must be >0 */
  for(j=0; (i+j)<len; j++) if(!isdigit((int)studynr[i+j])) return(0);
  //if(j<1 || j>5) return(0);
  if(j<1 || (i+j)>10) return(0);
  return(1);
}

Referenced by studynr_in_fname(), and studynr_validity_check().

◆ swabip()

void swabip	(	void *	buf,
		long long int	size )

extern

In-place swab, replaces the non-ANSI function swab(), which may not work if data is overlapping.

Parameters

buf	Pointer to memory
size	Size of buf in bytes

Definition at line 72 of file swap.c.

  {
  unsigned char c;
 
  for(long long i=1; i<size; i+=2) {
    c=*((unsigned char *)buf+i);
    *((unsigned char *)buf+i)=*((unsigned char *)buf+(i-1));
    *((unsigned char *)buf+(i-1))=c;
  }
}

Referenced by anaReadHeader(), anaReadImagedata(), anaWriteHeader(), dcmFileWrite(), dcmWriteFileTag(), ecat63ReadAttnheader(), ecat63ReadImageheader(), ecat63ReadMainheader(), ecat63ReadMatdata(), ecat63ReadNormheader(), ecat63ReadScanheader(), ecat63rFloat(), ecat63wFloat(), ecat63WriteAttnheader(), ecat63WriteImageheader(), ecat63WriteMainheader(), ecat63WriteMatdata(), ecat63WriteNormheader(), ecat63WriteScanheader(), ecat7Read2DNormheader(), ecat7Read2DScanheader(), ecat7ReadAttenheader(), ecat7ReadImageheader(), ecat7ReadMainheader(), ecat7ReadMatrixdata(), ecat7ReadNormheader(), ecat7ReadPolmapheader(), ecat7ReadScanheader(), ecat7rFloat(), ecat7Write2DNormheader(), ecat7Write2DScanheader(), ecat7WriteAttenheader(), ecat7WriteImageheader(), ecat7WriteMainheader(), ecat7WriteMatrixdata(), ecat7WriteNormheader(), ecat7WritePolmapheader(), ecat7WriteScanheader(), imgWriteAnalyze(), imgWriteAnalyzeFrame(), niftiReadHeader(), niftiReadImagedata(), niftiWriteHeader(), and upetReadImagedata().

◆ swap()

void swap	(	void *	from,
		void *	to,
		int	size )

extern

Swaps the specified short int, int, long int, float, or double from little endian to big endian or vice versa. Arguments are allowed to overlap.

Parameters

from	Pointer to a short int, int, long int, float, or double variable
to	Pointer to a short int, int, long int, float, or double variable
size	Size of from and to (byte nr) must be 1, 2, 4 or 8.

Definition at line 31 of file swap.c.

                                          {
  unsigned char c;
  unsigned short int s;
  unsigned long l;
  
  switch(size) {
    case 1:
      *(char *)to=*(char *)from;
      break;
    case 2:
      c=*(unsigned char *)from;
      *(unsigned char *)to = *((unsigned char *)from+1);
      *((unsigned char *)to+1) = c; 
      /*swab(from, to, size); // NOT ANSI */
      break;
    case 4:
      s=*(unsigned short *)from;
      *(unsigned short *)to = *((unsigned short *)from+1);
      *((unsigned short *)to+1) = s;
      swap((char*)to, (char*)to, 2);
      swap((char*)((unsigned short *)to+1), (char*)((unsigned short *)to+1), 2);
      break;
    case 8:
      l=*(unsigned long *)from;
      *(unsigned long *)to = *((unsigned long *)from+1);
      *((unsigned long *)to+1) = l;
      swap((char *)to, (char *)to, 4);
      swap((char*)((unsigned long *)to+1), (char*)((unsigned long *)to+1), 4);
      break;
  }
}

Referenced by dcmitemGetInt(), dcmitemGetReal(), dcmReadFileTag(), dcmReadFileVL(), dcmReadFileVRVL(), dcmValueString(), dcmWriteFileVRVL(), intExpand(), and swap().

◆ swawbip()

void swawbip	(	void *	buf,
		long long int	size )

extern

In-place swab and swaw, switches words and bytes from an array of 4-byte ints or floats.

Parameters

buf	Pointer to memory
size	Size of buf in bytes

Definition at line 93 of file swap.c.

  {
  unsigned char c, *cptr;
 
  cptr=(unsigned char*)buf;
  for(long long i=0; i<size; i+=4, cptr+=4) {
    c=cptr[0]; cptr[0]=cptr[3]; cptr[3]=c;
    c=cptr[1]; cptr[1]=cptr[2]; cptr[2]=c;
  }
}

Referenced by anaReadHeader(), anaReadImagedata(), anaWriteHeader(), ecat63Create(), ecat63Matenter(), ecat63ReadMatdata(), ecat63ReadMatlist(), ecat63rFloat(), ecat63rInt(), ecat63wFloat(), ecat63wInt(), ecat63WriteMatdata(), ecat7Create(), ecat7EnterMatrix(), ecat7Read2DNormheader(), ecat7Read2DScanheader(), ecat7ReadAttenheader(), ecat7ReadImageheader(), ecat7ReadMainheader(), ecat7ReadMatlist(), ecat7ReadMatrixdata(), ecat7ReadNormheader(), ecat7ReadPolmapheader(), ecat7ReadScanheader(), ecat7rFloat(), ecat7rInt(), ecat7Write2DNormheader(), ecat7Write2DScanheader(), ecat7WriteAttenheader(), ecat7WriteImageheader(), ecat7WriteMainheader(), ecat7WriteMatrixdata(), ecat7WriteNormheader(), ecat7WritePolmapheader(), ecat7WriteScanheader(), hrrtMakeCalHdr(), niftiReadHeader(), niftiReadImagedata(), niftiWriteHeader(), and upetReadImagedata().

◆ swawip()

void swawip	(	void *	buf,
		long long int	size )

extern

In-place swaw, switches words (but not bytes) from an array of 4-byte ints or floats.

Parameters

buf	Pointer to memory
size	Size of buf in bytes

Definition at line 114 of file swap.c.

  {
  unsigned short int s, *sptr;
 
  sptr=(unsigned short int*)buf;
  for(long long i=0; i<size; i+=4, sptr+=2) {
    s=sptr[0]; sptr[0]=sptr[1]; sptr[1]=s;
  }
}

Referenced by ecat63rFloat(), ecat63wFloat(), and ecat7rFloat().

◆ temp_roundf()

int temp_roundf ( float e )

extern

int roundf(float e) - Rounds up float e to nearest int

Parameters

e	float value

Returns: rounded integer

Definition at line 20 of file petc99.c.

{
#if defined(__STDC_VERSION__) && __STD_VERSION__>=199901L
  return(roundf(e));
#else
  if(e<0.0) {
      return (int)(e-0.5);
  } else {
      return (int)(e+0.5);
  }
#endif
}

Referenced by bootstrap(), ecat63AddImg(), ecat63WriteAllImg(), ecat63WriteImageMatrix(), ecat63WriteScanMatrix(), ecat7Write2DScanMatrix(), ecat7WriteImageMatrix(), ecat7WritePolarmapMatrix(), ecat7WriteScanMatrix(), ecatCopy63to7imageheader(), img2svol(), imgMaskRoiNr(), imgSetAnalyzeHeader(), imgVoiMaskTAC(), imgWrite2DEcat7(), imgWriteAnalyze(), imgWriteAnalyzeFrame(), imgWriteEcat63Frame(), imgWriteEcat7(), imgWriteEcat7Frame(), imgWritePolarmap(), roi_onoff(), roiFillGaps(), and tsvRead().

◆ textfileReadLines()

int textfileReadLines	(	const char *	filename,
		STR_TOKEN_LIST *	lst )

extern

Read all lines from text file into STR_TOKEN_LIST. List needs to be initialized. Previous contents are deleted.

Returns: Returns 0, if successful, and <>0 in case of an error.

See also: str_token_list_init, str_token_list_empty

Parameters

filename	Name of text file to read
lst	Token list is allocated by this function.

Definition at line 136 of file readfile.c.

  {
  int i, ret, nr=0;
  char *allfile, *cptr, *line;
  FILE *fp;
 
  if(lst==NULL || filename==NULL || strlen(filename)<1) return(1);
  str_token_list_empty(lst);
 
  /* Open file */
  fp=fopen(filename, "r"); if(fp==NULL) return(2);
  /* Get file size */
  nr=0; while((ret=fgetc(fp))!=EOF) nr++; rewind(fp);
  if(nr<1) {fclose(fp); return(0);} //printf("nr=%d\n", nr);
  /* Allocate memory for file contents */
  allfile=(char*)malloc((nr+1)*sizeof(char));
  if(allfile==NULL) {fclose(fp); return(3);}
  /* Read file contents */
  i=0; while((ret=fgetc(fp))!=EOF && i<nr) allfile[i++]=(char)ret;
  fclose(fp); allfile[i]=(char)0;
  /* and then fill the list */
  int npos=0; cptr=allfile;
  line=strTokenDup(cptr, "\n\r\0", &npos);
  if(line==NULL) {free(allfile); return(4);}
  do {
    //printf("line='%s'\n", line);
    if(str_token_list_add(lst, line)) {free(allfile); free(line); return(10);}
    free(line); if(npos==0) break; cptr=cptr+npos;
    line=strTokenDup(cptr, "\n\r\0", &npos);
  } while(line!=NULL);
  free(allfile);
 
  return(0);
}

Referenced by roi_read(), and sifRead().

◆ time_to_tm()

void time_to_tm	(	time_t	totalsecs,
		int	offset,
		struct tm *	result )

extern

Convert calendar time to local broken-down time. This function is copied from GNU C Library with tiny modifications.

Parameters

totalsecs	number of seconds elapsed since 00:00:00 on January 1, 1970, UTC; can be negative to represent times before 1970
offset	offset seconds adding to totalsecs (e.g. -timezone)
result	pointer to struct tm variable to receive broken-down time

Definition at line 452 of file datetime.c.

  {
  if(offset==0) {gmtime_r(&totalsecs, result); return;}
 
 
  long int days, rem, y, yg;
  const unsigned short *ip;
 
  days=totalsecs/86400; rem=totalsecs%86400; rem+=offset;
  while(rem<0) {rem+=86400; --days;}
  while(rem>=86400) {rem-=86400; ++days;}
  result->tm_hour=rem/3600;
  rem%=3600; result->tm_min=rem/60;
  result->tm_sec=rem%60;
 
  /* January 1, 1970 was a Thursday. */
  result->tm_wday=(4+days)%7;
  if(result->tm_wday<0) result->tm_wday+=7;
  y=1970;
  while(days<0 || days>=(isleapyear(y)?366:365)) {
    /* Guess a corrected year, assuming 365 days per year. */
    yg=y+math_div(days, 365);
    /* Adjust days and y to match the guessed year. */
    days-=(yg-y)*365 + leaps_between(y, yg);
    y=yg;
  }
  result->tm_year=y-1900; result->tm_yday=days;
  ip=__mon_yday[isleapyear(y)];
  for(y=11; days<ip[y]; y--) continue;
  days-=ip[y];
  result->tm_mon=y;
  result->tm_mday=days+1;
  result->tm_isdst=-1;
#if defined(HAVE_TM_GMTOFF)
  result->tm_gmtoff=0L;
#endif
}

◆ timegm()

time_t timegm ( struct tm * tm )

extern

Inverse of gmtime, converting struct tm to time_t.

Otherwise same as mktime, except that mktime uses local time. Uses gmtime_r or gmtime_s, if available, otherwise gmtime, which is threadsafe in Windows.

Returns: Returns the time_t, or -1 in case of an error.

Parameters

tm	Pointer to struct tm

Definition at line 69 of file datetime.c.

  {
#ifdef HAVE_GMTIME_R
  {
  time_t temp_lt;
  struct tm temp_gm;
  if(!tm) temp_lt=0; else temp_lt=mktime(tm);
  if(gmtime_r(&temp_lt, &temp_gm)==NULL) return((time_t)-1);
  return(time_t)(temp_lt + (temp_lt - mktime(&temp_gm)));
  }
#endif
 
#ifdef HAVE_GMTIME_S
  {
  time_t temp_lt;
  struct tm temp_gm;
  if(!tm) temp_lt=0; else temp_lt=mktime(tm);
  if(gmtime_s(&temp_gm, &temp_lt)!=0) return(time_t)-1);
  return(time_t)(temp_lt + (temp_lt - mktime(&temp_gm)));
  }
#endif
 
  time_t temp_lt;
  struct tm *temp_gm;
  if(!tm) temp_lt=0; else {tm->tm_isdst=-1; temp_lt=mktime(tm);}
  temp_gm=gmtime(&temp_lt); if(temp_gm) temp_gm->tm_isdst=-1;
  return(time_t)(temp_lt + (temp_lt - mktime(temp_gm)));
}

Referenced by ecat63Scanstarttime(), ecat63ScanstarttimeToTm(), ecat7EditMHeader(), fitRead(), get_datetime(), imgGetMicropetHeader(), resRead(), sifRead(), tmAdd(), and upetScanStart().

◆ tmAdd()

void tmAdd	(	int	s,
		struct tm *	d )

extern

Add given time in seconds to the date and time.

Parameters

s	Time to add in seconds; can be negative
d	Pointer to tm struct

Definition at line 514 of file datetime.c.

  {
  if(d==NULL) return;
  d->tm_sec+=s; 
  //mktime(d); // this automatically normalizes sec to hours etc if necessary
  timegm(d); // this automatically normalizes sec to hours etc if necessary
}

◆ tmDifference()

double tmDifference	(	struct tm *	tm1,
		struct tm *	tm0 )

extern

Calculate the difference in seconds between two given dates and times.

Returns: Returns tm1-tm0 in seconds.

Parameters

tm1	Pointer to tm struct
tm0	Pointer to tm struct

Definition at line 502 of file datetime.c.

  {
  return(difftime(mktime(tm1), mktime(tm0)));
}

Referenced by selectEcat931Calibrationfile().

◆ tpcHtmlUsage()

int tpcHtmlUsage	(	const char *	program,
		char *	text[],
		const char *	path )

extern

Write program usage given as argument, plus program name, tpcclib version, and default copyright text, into HTML file.

Any string @P, separated by space characters, is replaced by program name in the output. It can only be used once per line.
When line contains text 'stdoptions', in place of that the description of standard command-line options '-h, -v, -q, -s etc' are displayed.
Returns
Returns 0 when successful.

Author
Vesa Oikonen

See also
tpcPrintUsage

Parameters

program	Program name, may contain extension and path.
text	Program usage text.
path	Path name where to create file programname.html; path may contain trailing '/' or '\'.

Definition at line 213 of file proginfo.c.

  {
  unsigned int len, i, j;
  char *bprogram, *fname, *cptr, *line;
  FILE *fp;
 
  if(program==NULL || text==NULL || strnlen(program, 1)<1) return 1;
  
  /* Clean program name */
  bprogram=strdup(program);
  filenameRmPath(bprogram); filenameRmExtension(bprogram);
  
  /* Make file name */
  fname=calloc(strlen(path)+1+strlen(bprogram)+5, sizeof(char));
  if(fname==NULL) {free(bprogram); return 1;}
  strcpy(fname, path); len=strlen(fname); 
  if(len>0 && (fname[len-1]=='/' || fname[len-1]=='\\')) fname[len-1]=(char)0;
  len=strlen(fname); if(len>0) strcat(fname, "/"); 
  strcat(fname, bprogram); strcat(fname, ".html");
  
  //printf("fname := '%s'\n", fname);
 
  /* Open file for write */
  fp=stdout;
 
  /* Write HTML header */
  len=fprintf(fp, "<!DOCTYPE html>\n");
  if(len<10) {free(bprogram); free(fname); return 2;}
  fprintf(fp, "<html lang=\"en-GB\">\n");
  fprintf(fp, "<head>\n");
  fprintf(fp, "  <meta charset=\"utf-8\">\n");
  fprintf(fp, "  <meta name=\"viewport\" content=\"width=device-width, initial-scale=1.0\">\n");
  fprintf(fp, "  <title>%s</title>\n", bprogram);
  fprintf(fp, "  <style type=\"text/css\">\n");
  fprintf(fp, "    body {\n");
  fprintf(fp, "      margin-left: 2em;\n");
  fprintf(fp, "      font-family: monospace;\n");
  fprintf(fp, "      font-size: 1em;\n");
  fprintf(fp, "    }\n");
  fprintf(fp, "    h1 {\n");
  fprintf(fp, "      font-size: 1.3em;\n");
  fprintf(fp, "      margin-top: 1em;\n");
  fprintf(fp, "      margin-bottom: 1em;\n");
  fprintf(fp, "    }\n");
  fprintf(fp, "    footer {\n");
  fprintf(fp, "      border:1px solid gray;\n");
  fprintf(fp, "      font-size: 0.8em;\n");
  fprintf(fp, "    }\n");
  fprintf(fp, "    footer p {margin-left: 1em;}\n");
  fprintf(fp, "  </style>\n");
  fprintf(fp, "</head>\n\n");
  
  /* Write HTML body */
  fprintf(fp, "<body>\n");
  
  /* Write program name, version and copyright as title; */
  /* replace (c) with html code when necessary */
  fprintf(fp, "<h1>%s - tpcclib %d.%d.%d ", bprogram, tpcclib_VERSION_MAJOR,
          tpcclib_VERSION_MINOR, tpcclib_VERSION_PATCH);
  line=tpcclib_COPYRIGHT; len=strlen(line);
  for(j=0; j<len; j++) {
    if(strncasecmp(line+j, "(C)", 3)==0) {fputs("&copy;", fp); j+=2; continue;}
    fputc(line[j], fp);
  }
  fputs("</h1>\n\n", fp);
 
  /* Print usage */
  fprintf(fp, "<pre>\n");
  i=0; while(text[i]!=0) {
    line=text[i];
    /* If line contains string 'stdoptions' then print instead of it the
       description of standard command-line options */
    if(strstr(line, "stdoptions")) {
      int j=0;
      while(tpcstdoptions[j]!=0) fprintf(fp, "%s\n", tpcstdoptions[j++]);        
      i++; continue;
    }
 
    /* Process "See also" line: add links to other programs */
    if(strstr(line, "See also: ")!=NULL) {
      /* copy until the first ':' */
      j=0; while(line[j]!='\0') {
        fputc(line[j], fp);
        j++; if(line[j-1]==':') break;
      }
      /* the rest of line with token */
      char *tline; unsigned int n=0;
      tline=strdup(line+j); cptr=strtok(tline, ", :;\t\n\r");
      while(cptr!=NULL) {
        if(n>0) fputc(',', fp);
        fprintf(fp, " <a href=\"./%s.html\">%s</a>", cptr, cptr);
        cptr=strtok(NULL, ", :;\t\n\r"); 
        n++;
      }
      fputs("\n", fp);
      free(tline); 
      i++; continue;
    }
 
    
    /* Print the line one character at the time */
    len=strlen(line); j=0;
    while(j<len) {
      
      /* If WWW Address follows, then add a link */
      if(strncasecmp(line+j, "https://", 7)==0) {
        unsigned int li;
        cptr=line+j; len=strcspn(cptr, " ),;");
        fputs("<a href=\"", fp);
        for(li=0; li<len; li++) fputc(line[j+li], fp);
        fputs("\">", fp);
        for(li=0; li<len; li++) fputc(line[j+li], fp);
        fputs("</a>", fp);
        j+=len; continue;
      }
      
      /* If necessary, replace '@P' with program name */
      if(strncmp(line+j, " @P ", 4)==0) {
        fprintf(fp, " %s ", bprogram);
        j+=4; continue;
      }
      /* Replace (c) or (C) with html code */
      if(strncasecmp(line+j, "(C)", 3)==0) {
        fputs("&copy;", fp);
        j+=3; continue;
      }
 
      /* Replace <, >, and & characters with html codes */
      if(line[j]=='<') {fputs("&lt;", fp); j++; continue;}
      if(line[j]=='>') {fputs("&gt;", fp); j++; continue;}
      if(line[j]=='&') {fputs("&amp;", fp); j++; continue;}
 
      /* Just write the char normally */
      fputc(line[j], fp); j++;
    }
    fprintf(fp, "\n");
    i++; continue;    
 
  }
  fprintf(fp, "</pre>\n");
 
  /* Write footer */
  fprintf(fp, "\n<footer>\n");
  fprintf(fp, "<p>");
  i=0; while(tpclicense4html[i]!=0) fprintf(fp, "%s<br>\n", tpclicense4html[i++]);
  fprintf(fp, "</p>\n");
  fprintf(fp, "</footer>\n");
           
  /* Close HTML */
  fprintf(fp, "</body>\n");
  fprintf(fp, "</html>\n");
 
 
  free(bprogram); free(fname);
  return 0;
}

◆ tpcPrintBuild()

void tpcPrintBuild	(	const char *	program,
		FILE *	fp )

extern

Print tpctools build information.

Author: Vesa Oikonen

See also: tpcProgramName, tpcPrintUsage

Parameters

program	Program name; enter NULL, if not to be printed.
fp	File pointer where to print; usually stdout.

Definition at line 383 of file proginfo.c.

  {
  fprintf(fp, "\n");
  if(program!=NULL) {
    /* Print program name */
    char *s; s=strdup(program); 
    filenameRmPath(s); filenameRmExtension(s);
    fprintf(fp, " Program: %s\n", s);
    free(s);
  }
  /* Build time */
  fprintf(fp, " Build: %s %s\n",__DATE__,__TIME__);
  /* tpcclib (and program) version */
  fprintf(fp, " tpcclib version: %d.%d.%d\n", tpcclib_VERSION_MAJOR,
           tpcclib_VERSION_MINOR, tpcclib_VERSION_PATCH);
  /* Compiler information */
  fprintf(fp, " Build platform: %s\n", tpcclib_BUILD_HOST_SYSTEM_NAME);
  fprintf(fp, " Build processor: %s\n", tpcclib_BUILD_HOST_SYSTEM_PROCESSOR);
#if defined(__STDC_VERSION__)
  fprintf(fp, " Version of C: %ld\n", __STDC_VERSION__);
#endif
#if defined(__GNUC__) && defined(__VERSION__)
  fprintf(fp, " GNU C version: %s\n", __VERSION__);
#endif
#if defined(__clang__) && defined(__clang_version__)
  fprintf(fp, " Clang/LLVM version: %s\n", __clang_version__);
#endif
#ifdef _OPENMP
  fprintf(fp, " OpenMP version: %d\n", _OPENMP);
#endif
  /* Platform information */
#if defined(__x86_64__) || defined(__LP64__) || defined(__ppc64__) || \
    defined(__LLP64__) || defined(__ILP64__) 
  fprintf(fp, " Architecture: 64-bit\n");
#else
  fprintf(fp, " Architecture: 32-bit\n");
#endif
  /* Check if IEC 60559 floating-point standard is officially supported;
     that is 'always' supported, but PET data i/o functions rely on it.
     GCC does not define __STDC_IEC_559__ or implement the associated standard pragmas. */
#ifdef __STDC_IEC_559__
  fprintf(fp, " IEC 60559 floating-point standard is fully supported.\n");
#endif
  /* Large File Support extension (LFS) */
#if defined(_FILE_OFFSET_BITS)
  fprintf(fp, " _FILE_OFFSET_BITS: %d\n", _FILE_OFFSET_BITS);
#else
  fprintf(fp, " _FILE_OFFSET_BITS not defined\n");
#endif
#if defined(_LARGEFILE_SOURCE)
  fprintf(fp, " _LARGEFILE_SOURCE defined\n");
#endif
#if defined(_LARGEFILE64_SOURCE)
  fprintf(fp, " _LARGEFILE64_SOURCE defined\n");
#endif
}

◆ tpcPrintUsage()

void tpcPrintUsage	(	const char *	program,
		char *	text[],
		FILE *	fp )

extern

Print program usage given as argument, plus program name, tpcclib version, and default copyright text.

Any string @P, separated by space characters, is replaced by program name in the output. It can only be used once per line.
When line contains text 'stdoptions', in place of that the description of standard command-line options '-h, -v, -q, -s etc' are displayed.
Author
Vesa Oikonen

See also
tpcHtmlUsage, tpcPrintBuild, tpcProcessStdOptions

Parameters

program	Program name
text	Program usage text
fp	File pointer where to print; usually stdout

Definition at line 158 of file proginfo.c.

  {
  int i;
  char *cptr, *bprogram;
 
  /* Print program name, version, and copyright */
  if(strlen(program)>0) bprogram=strdup(program);
  else bprogram=strdup("unknown"); 
  filenameRmPath(bprogram); filenameRmExtension(bprogram);
  fprintf(fp, "\n  %s - tpcclib %d.%d.%d %s\n \n", bprogram, 
           tpcclib_VERSION_MAJOR, tpcclib_VERSION_MINOR, tpcclib_VERSION_PATCH,
           tpcclib_COPYRIGHT);
  /* Print usage */
  i=0; while(text[i]!=0) {
    /* If line contains string 'stdoptions' then print instead of it the
       description of standard command-line options */
    if(strstr(text[i], "stdoptions")) {
      int j=0;
      while(tpcstdoptions[j]!=0) fprintf(fp, "  %s\n", tpcstdoptions[j++]);        
      i++; continue;
    }
    /* If line does not contain program name, then just print it as it is */
    cptr=strstr(text[i], " @P ");
    if(cptr==NULL) {fprintf(fp, "  %s\n", text[i++]); continue;}
    /* Replace '@P' with program name */
    char *s; s=strdup(text[i]);
    s[strlen(text[i])-strlen(cptr)]=(char)0;
    fprintf(fp, "  %s %s %s\n", s, bprogram, cptr+4);
    free(s); i++;
  }
  fprintf(fp, " \n");
  /* Print licence info */
  i=0; while(tpclicense[i]!=0) fprintf(fp, "  %s\n", tpclicense[i++]);
  fprintf(fp, "\n");
  free(bprogram);
}

◆ tpcProcessStdOptions()

int tpcProcessStdOptions	(	const char *	s,
		int *	print_usage,
		int *	print_version,
		int *	verbose_level )

extern

Check if given command-line argument string is one of the standard command-line options of this project.

Returns: 0 if string was identified as standard option, otherwise 1.

Author: Vesa Oikonen

See also: tpcProgramName, tpcPrintUsage

Parameters

s	Pointer to command-line option string.
print_usage	If option string is either -h or –help, then this variable is set to 1.
print_version	If option string is either -v, -V, –version, or –build, then this variable is set to 1.
verbose_level	The level of debugging messages and listings: If option string is -d, –debug or –verbose, then +1 is added to this variable. If option string is -d[n], –debug[=n], or –verbose[=n], then +n is added to this variable. If options string is -q or –quiet, then this variable is set to 0. If options string is -s or –silent, then this variable is set to -1.

Definition at line 40 of file proginfo.c.

  {
  int n;
 
  char *cptr;
  /* Check that string is option, starting with '-' or '--' */
  if(s==NULL || strlen(s)<2 || s[0]!='-') return 1;
  /* Set pointer to the character after the first '-' */
  cptr=(char*)s+1;
  /* If also the next character is '-', then try the long forms of options */
  if(*cptr=='-') {
    cptr++; if(strlen(cptr)<1) return 1;
    if(strcasecmp(cptr, "help")==0) {*print_usage=1; return 0;}
    if(strcasecmp(cptr, "helphtml")==0) {*print_usage=2; return 0;}
    if(strcasecmp(cptr, "version")==0) {*print_version=1; return 0;}
    if(strcasecmp(cptr, "build")==0) {*print_version=1; return 0;}
    if(strcasecmp(cptr, "debug")==0) {*verbose_level+=1; return 0;}
    if(strcasecmp(cptr, "verbose")==0) {*verbose_level+=1; return 0;}
    if(strncasecmp(cptr, "debug=", 6)==0) {
      if(!isdigit(cptr[6])) return 1;
      n=atoi(cptr+6); *verbose_level+=n; return 0;
    }
    if(strncasecmp(cptr, "verbose=", 8)==0) {
      if(!isdigit(cptr[8])) return 1;
      n=atoi(cptr+8); *verbose_level+=n; return 0;
    }
    if(strcasecmp(cptr, "quiet")==0) {*verbose_level=0; return 0;}
    if(strcasecmp(cptr, "silent")==0) {*verbose_level=-1; return 0;}
    return 1;
  }
  /* So it is the short form, if anything */
  if(strcmp(cptr, "h")==0) {*print_usage=1; return 0;}
  if(strcasecmp(cptr, "v")==0) {*print_version=1; return 0;}
  if(strcmp(cptr, "d")==0) {*verbose_level+=1; return 0;}
  if(strncmp(cptr, "d", 1)==0 && strlen(cptr)>1) {
    if(!isdigit(cptr[1])) return 1;
    n=atoi(cptr+1); *verbose_level+=n; return 0;
  }
  if(strcmp(cptr, "q")==0) {*verbose_level=0; return 0;}
  if(strcmp(cptr, "s")==0) {*verbose_level=-1; return 0;}
  return 1;
}

◆ tpcProgramName()

void tpcProgramName	(	const char *	program,
		int	version,
		int	copyright,
		char *	prname,
		int	n )

extern

Process program name and optionally version into given string from argv[0].

See also: tpcPrintBuild, tpcPrintUsage, tpcProcessStdOptions

Parameters

program	Set to argv[0]
version	Add version (1) or do not add (0)
copyright	Add copyright (1) or do not add (0)
prname	Pointer to string where program name is written
n	Length of prname string, including trailing zero

Definition at line 101 of file proginfo.c.

  {
  char *tmp;
 
  /* Check the input */
  if(prname==NULL || n<1) return;
  prname[0]=(char)0;
  
  /* Remove path and extension */
  if(strlen(program)>0) tmp=strdup(program); else tmp=strdup("unknown"); 
  filenameRmPath(tmp); filenameRmExtension(tmp);
  /* Copy it if possible */
  n-=strlen(tmp); if(n>0) strcpy(prname, tmp); else {free(tmp); return;}
  free(tmp);
 
  /* Add version, if required */
  if(version!=0) {
    /* Create string with version number */
    char v[256];
    sprintf(v, "%d.%d.%d", tpcclib_VERSION_MAJOR, tpcclib_VERSION_MINOR, tpcclib_VERSION_PATCH);
    /* Copy it if possible */
    n--; // space
    n-=strlen(v); if(n>0) {strcat(prname, " "); strcat(prname, v);} 
  }
  
  /* Add copyright, if required */
  if(copyright!=0) {
    /* Copy it if possible */
    n--; // space
    n-=strlen(tpcclib_COPYRIGHT); 
    if(n>0) {strcat(prname, " "); strcat(prname, tpcclib_COPYRIGHT);} 
  }
  
  return;
}

Data Structures

Macros

Enumerations

Functions

Detailed Description

Macro Definition Documentation

◆ BRANCHING_C

◆ BRANCHING_Cu64

◆ BRANCHING_F

◆ BRANCHING_Ga

◆ BRANCHING_Ge

◆ BRANCHING_N

◆ BRANCHING_O

◆ BRANCHING_Rb

◆ HL_Br75

◆ HL_Br76

◆ HL_C11

◆ HL_Cu62

◆ HL_Cu64

◆ HL_F18

◆ HL_Fe52

◆ HL_Ga68

◆ HL_Ge68

◆ HL_I124

◆ HL_N13

◆ HL_Na22

◆ HL_O14

◆ HL_O15

◆ HL_Rb82

◆ HL_Zn62

◆ M_LN2

◆ MAX_REGIONNAME_LEN

◆ MAX_REGIONSUBNAME_LEN

◆ MAX_STUDYNR_LEN

◆ MAX_UNITS_LEN

Enumeration Type Documentation

◆ anonymous enum

◆ anonymous enum

◆ anonymous enum

◆ anonymous enum

Function Documentation

◆ _intexadd()

◆ asciiCommentLine()

◆ atof_dpi()

◆ atof_with_check()

◆ atoi_with_check()

◆ backupExistingFile()

◆ branchingFraction()

◆ ctime_r_int()

◆ cunitFromFilename()

◆ dec_comma_is()

◆ dec_nr()

◆ dec_separator()

◆ dec_separator_change()

◆ defRead()

◆ doubleCopy()

◆ doubleCSpanPositives()

◆ doubleMachEps()

◆ doubleMatch()

◆ doubleMatchRel()

◆ doubleMaxIndex()

◆ doubleMean()

◆ doubleSpanPositives()

◆ doubleSum()

◆ fileCopy()

◆ filenameGetExtension()

◆ filenameGetExtensions()

◆ filenameRmExtension()

◆ filenameRmExtensions()

◆ filenameRmPath()

◆ fncasematch()

◆ fnmatch()

◆ get_date()

◆ get_datetime()

◆ gmtime_r()

◆ hl2lambda()

◆ hlCorrectIsotopeCode()

◆ hlFromIsotope()

◆ hlIsotopeCode()

◆ hlIsotopeFromHalflife()