Header file for library libtpcextensions. More...

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

Data Structures
struct	TPCSTATUS

struct	INTLIST

Macros
#define	MAX_UNITS_LEN 32
	Define max units string length.

#define	MAX_TACNAME_LEN 31
	Max length of TAC ID name (not including trailing zero)

#define	MAX_PARNAME_LEN MAX_TACNAME_LEN
	Max string length for PAR name.

#define	MAX_STUDYNR_LEN 255
	Define max study number length.

Enumerations
enum	weights { WEIGHTING_UNKNOWN , WEIGHTING_OFF , WEIGHTING_ON_GENERAL , WEIGHTING_ON_COUNTS , WEIGHTING_ON_F , WEIGHTING_ON_FD , WEIGHTING_LAST }
	Is data weighted, or are weight factors available with data? More...

enum	tactype { TACTYPE_UNKNOWN , TACTYPE_PLASMA , TACTYPE_BLOOD , TACTYPE_PARENT , TACTYPE_METABOLITE , TACTYPE_VOI , TACTYPE_REF , TACTYPE_PIXEL , TACTYPE_LAST }
	TAC content types. More...

enum	unit { UNIT_UNKNOWN , UNIT_UNITLESS , UNIT_HU , UNIT_MSEC , UNIT_SEC , UNIT_MIN , UNIT_HOUR , UNIT_DAY , UNIT_MONTH , UNIT_YEAR , UNIT_UM , UNIT_MM , UNIT_CM , UNIT_M , UNIT_UL , UNIT_ML , UNIT_DL , UNIT_L , UNIT_UG , UNIT_MG , UNIT_G , UNIT_100G , UNIT_KG , UNIT_PMOL , UNIT_NMOL , UNIT_UMOL , UNIT_MMOL , UNIT_MOL , UNIT_COUNTS , UNIT_KCOUNTS , UNIT_BQ , UNIT_KBQ , UNIT_MBQ , UNIT_GBQ , UNIT_NCI , UNIT_UCI , UNIT_MCI , UNIT_CI , UNIT_CPS , UNIT_KCPS , UNIT_CPM , UNIT_KCPM , UNIT_BQ_PER_ML , UNIT_KBQ_PER_ML , UNIT_MBQ_PER_ML , UNIT_NCI_PER_ML , UNIT_UCI_PER_ML , UNIT_BQ_PER_G , UNIT_KBQ_PER_G , UNIT_MBQ_PER_G , UNIT_NCI_PER_G , UNIT_UCI_PER_G , UNIT_PMOL_PER_ML , UNIT_NMOL_PER_ML , UNIT_UMOL_PER_ML , UNIT_MMOL_PER_ML , UNIT_PMOL_PER_L , UNIT_NMOL_PER_L , UNIT_UMOL_PER_L , UNIT_MMOL_PER_L , UNIT_MOL_PER_L , UNIT_SEC_KBQ_PER_ML , UNIT_MIN_KBQ_PER_ML , UNIT_SEC_BQ_PER_ML , UNIT_MIN_BQ_PER_ML , UNIT_PERCENTAGE , UNIT_ML_PER_ML , UNIT_ML_PER_DL , UNIT_G_PER_ML , UNIT_ML_PER_G , UNIT_PER_SEC , UNIT_PER_MIN , UNIT_PER_HOUR , UNIT_ML_PER_ML_SEC , UNIT_ML_PER_ML_MIN , UNIT_ML_PER_DL_MIN , UNIT_ML_PER_G_MIN , UNIT_ML_PER_100G_MIN , UNIT_UMOL_PER_ML_MIN , UNIT_UMOL_PER_DL_MIN , UNIT_UMOL_PER_G_MIN , UNIT_UMOL_PER_100G_MIN , UNIT_MMOL_PER_ML_MIN , UNIT_MMOL_PER_DL_MIN , UNIT_MMOL_PER_G_MIN , UNIT_MMOL_PER_100G_MIN , UNIT_MG_PER_100G_MIN , UNIT_MG_PER_DL_MIN , UNIT_PID , UNIT_PID_PER_G , UNIT_PID_PER_ML , UNIT_PID_PER_KG , UNIT_PID_PER_L , UNIT_MBQ_PER_NMOL , UNIT_GBQ_PER_NMOL , UNIT_LAST }

enum	tpcerror { TPCERROR_OK , TPCERROR_FAIL , TPCERROR_OUT_OF_MEMORY , TPCERROR_NO_FILE , TPCERROR_FILE_EXISTS , TPCERROR_INVALID_FILENAME , TPCERROR_CANNOT_OPEN , TPCERROR_CANNOT_READ , TPCERROR_CANNOT_WRITE , TPCERROR_CANNOT_DELETE , TPCERROR_INVALID_FORMAT , TPCERROR_UNSUPPORTED , TPCERROR_NO_DATA , TPCERROR_DUPLICATE_DATA , TPCERROR_NO_WEIGHTS , TPCERROR_NO_DATETIME , TPCERROR_MISSING_DATA , TPCERROR_TOO_FEW , TPCERROR_TOO_BIG , TPCERROR_OVERLAPPING_DATA , TPCERROR_LARGE_GAP , TPCERROR_NO_X , TPCERROR_INVALID_X , TPCERROR_INVALID_XRANGE , TPCERROR_INCOMPATIBLE_UNIT , TPCERROR_INCOMPATIBLE_DATA , TPCERROR_NO_REFERENCE , TPCERROR_NO_KEY , TPCERROR_NO_VALUE , TPCERROR_INVALID_VALUE , TPCERROR_INVALID_SEPARATOR , TPCERROR_INVALID_PARNR , TPCERROR_INVALID_HEADER , TPCERROR_MISSING_HEADER , TPCERROR_UNKNOWN_ISOTOPE , TPCERROR_UNKNOWN_UNIT , TPCERROR_NO_SOLUTION , TPCERROR_BAD_FIT , TPCERROR_LAST }
	Error codes. Error code tpcerror, same as error message position in the table (table can not be accessed directly outside the c file). More...

enum	{ ENDIAN_UNKNOWN , ENDIAN_BIG , ENDIAN_LITTLE , ENDIAN_BIG_WORD , ENDIAN_LITTLE_WORD }

Functions
int	strTokenNr (const char s1, const char s2)

int	strTokenNCpy (const char s1, const char s2, int i, char *s3, int count)

int	strChrCount (const char s1, const char s2)

int	strUppercaseCount (const char *s)

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

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

char *	strdup (const char *s)

char *	strndup (const char *s, size_t n)

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 *	strstrNoQuotation (const char haystack, const char needle)

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

int	strCleanSpaces (char *s)

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

int	strClean (char *s)

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

int	strInPars (char *s)

void	strCleanPars (char *s)

char *	strncatInt (char *str1, const int n, size_t count)

char *	strncatDouble (char *str1, const double d, size_t count)

char *	strncatIntZ (char *str1, const int n, const int maxn, size_t count)

int	strIsSpaceOnly (char *s)

char *	strdelstr (char s1, const char s2)

char *	strTrimLeft (char *s, const size_t t)

char *	strEncodeForXML (const char *s)

void	strCleanForXML (char *s)

int	strIsValidNumber (const char *s)

double	atofVerified (const char *s)

int	atofCheck (const char s, double v)

int	strHaveDecimalComma (const char *s)

int	strHaveDecimalSeparator (const char *s)

int	atofList (const char s1, const char s2, double *x, int maxn)

void	filenameRmPath (char *s)

void	filenameRmFile (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 (const char *s)
	Get the last extension of a file name.

char *	filenameGetExtensions (const char *s)
	Get all extensions of a file name.

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

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)

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

int	tpcYesNo (const char *s)

char *	unitName (int unit_code)

int	unitIdentify (const char *s)

int	unitIdentifyFilename (const char *s)

int	unitIsDistance (int u)

int	unitIsTime (int u)

int	unitIsVolume (int u)

int	unitIsMass (int u)

int	unitIsMole (int u)

int	unitIsRadioactivity (int u)

int	unitIsRAConc (int u)

int	unitIsCombinatorial (int u)

double	unitConversionFactor (const int u1, const int u2)

int	unitInverse (int u)

int	unitDividerHasVolume (int u)

int	unitDividerHasMass (int u)

int	unitDividendHasRadioactivity (int u)

int	unitCombination (const int u1, const int u2, const int v1, const int v2)

int	unitDividerMassVolumeConversion (int u)

int	unitMultiply (int ua, int ub)

void	statusPrint (FILE fp, TPCSTATUS s)

void	statusInit (TPCSTATUS *s)

void	statusFree (TPCSTATUS *s)

void	statusSet (TPCSTATUS s, const char func, const char *srcfile, int srcline, tpcerror error)

char *	errorMsg (tpcerror e)

size_t	asciiFileSize (FILE fp, int nonprintable)

char *	asciiFileRead (FILE fp, char data, size_t maxlen)

int	asciiCommentLine (const char line, int cont)

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 calendar 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	strDateValid (const char *str)

int	strDateValid2 (const char str, char intdate)

int	strDateValid3 (const char str, char intdate)

int	strDateValid4 (int dateint, char intdate, int year, int month, int day)

int	strTimeValid (const char *str)

int	strDateTimeValid (const char str, char intdate)

int	strDateTimeRead (const char str, struct tm date)

int	strDateRead (const char str, struct tm date)

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)

double	strDateTimeDifference (const char dt1, const char dt0)

int	strDateTimeAdd (int s, char *dt)

int	studynrVerify (const char *s, int zero_ok)

int	studynrStandardize (char *s)

int	studynrFromFilename (const char fname, char studynr, int force)

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

int	doubleArrayMatch (const double a1, const double a2, const unsigned int n, 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	doubleCopyFinite (double t, double s, const unsigned int n)

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

unsigned int	doubleRange (double a, const unsigned int n, double amin, double *amax)

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

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

double	doubleWMean (double a, double w, const unsigned int n)

int	doubleGetWithUnit (const char s, double v, int *u)

int	doubleSpanPositives (double *a, const int n)

int	doubleCSpanPositives (double *a, const int n)

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

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

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

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

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

unsigned int	doubleGEIndex (double *a, const unsigned int n, double lim)

unsigned int	doubleGTIndex (double *a, const unsigned int n, double lim)

double	inverfc (double x)

int	floatMatch (const float v1, const float v2, const float lim)

int	floatMatchRel (const float v1, const float v2, const float lim)

float	floatMachEps ()

void	floatCopy (float t, float s, const unsigned int n)

unsigned int	floatMaxIndex (float *a, const unsigned int n)

float	floatSum (float *a, const unsigned int n)

float	floatMean (float *a, const unsigned int n)

int	floatGetWithUnit (const char s, float v, int *u)

int	floatSpanPositives (float *a, const int n)

int	floatCSpanPositives (float *a, const int n)

unsigned int	floatNonzeroes (float *a, const unsigned int n)

int	atoiCheck (const char s, int v)

unsigned int	factorial (unsigned int n)

unsigned long long int	lfactorial (unsigned long long int n)

void	intlistInit (INTLIST *l)

void	intlistFree (INTLIST *l)

int	intlistAdd (INTLIST *l, const int v, const int ifnew)

void	intlistSort (INTLIST *l)

int	intlistAddFromString (const char s1, const char s2, INTLIST *l, const int ifnew)

int	intlistExpandFromString (const char s1, const char s2, INTLIST *l, const int ifnew)

char *	roinameSubpart (const char roiname, const char dlm, const unsigned int si, char *subpart, const unsigned int slen)

char *	roinameEditByTemplate (const char template, const char currname, char *newname, const unsigned int count)

char *	roinameAddField (char roiname, const char field, const unsigned int in, const unsigned int count)

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

int	roinameMatch (const char roiname, const char test_str, TPCSTATUS *status)

int	endianness ()

int	endianLittle ()

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

void	swabip (void *buf, int size)

void	swawbip (void *buf, int size)

void	swawip (void *buf, int size)

void	swap64ip (void *buf, unsigned long long size)

void	swap32ip (void *buf, unsigned long long size)

void	swap16ip (void *buf, unsigned long long size)

Detailed Description

Header file for library libtpcextensions.

Author: Vesa Oikonen

Copyright: (c) Turku PET Centre

Definition in file tpcextensions.h.

Macro Definition Documentation

◆ MAX_PARNAME_LEN

#define MAX_PARNAME_LEN MAX_TACNAME_LEN

Max string length for PAR name.

Max length of parameter names and units

Definition at line 39 of file tpcextensions.h.

Referenced by parEnsureNames(), parFromIFT(), parReadCSV(), parReadLimits(), parReadRES(), parWriteCSV(), and parWriteLimits().

◆ MAX_STUDYNR_LEN

#define MAX_STUDYNR_LEN 255

Define max study number length.

Max string length for PET study number

Definition at line 46 of file tpcextensions.h.

Referenced by imgCopyHeader(), imgInit(), imgReadDICOM(), tacGetHeaderStudynr(), tacRead(), tacReadSIF(), and tacWriteSIF().

◆ MAX_TACNAME_LEN

#define MAX_TACNAME_LEN 31

Max length of TAC ID name (not including trailing zero)

Max string length for TAC name

Definition at line 35 of file tpcextensions.h.

Referenced by parEnsureNames(), parFromIFT(), parReadCSV(), parReadFIT(), parReadRES(), tacEnsureNames(), tacRead4DM(), tacReadCarimasTxt(), tacReadCSV(), tacReadDFT(), tacReadInveonCSV(), tacReadMat(), tacReadPMOD(), tacReadQView(), and tacReadSimple().

◆ MAX_UNITS_LEN

#define MAX_UNITS_LEN 32

Define max units string length.

Max string length for PET units

Definition at line 28 of file tpcextensions.h.

Referenced by parWriteCSV().

Enumeration Type Documentation

◆ anonymous enum

anonymous enum

Definition at line 488 of file tpcextensions.h.

     { // used by endianness()
  ENDIAN_UNKNOWN,
  ENDIAN_BIG,
  ENDIAN_LITTLE,
  ENDIAN_BIG_WORD,   /* Middle-endian, Honeywell 316 style */
  ENDIAN_LITTLE_WORD /* Middle-endian, PDP-11 style */
};

◆ tactype

enum tactype

TAC content types.

Enumerator
TACTYPE_UNKNOWN	Content type not known.
TACTYPE_PLASMA	Plasma TAC.
TACTYPE_BLOOD	Blood TAC.
TACTYPE_PARENT	Parent tracer in plasma.
TACTYPE_METABOLITE	Metabolite in plasma.
TACTYPE_VOI	Image volume-of-interest TAC.
TACTYPE_REF	Reference region TAC.
TACTYPE_PIXEL	Image pixel TAC.

Definition at line 66 of file tpcextensions.h.

             {
  TACTYPE_UNKNOWN,      
  TACTYPE_PLASMA,       
  TACTYPE_BLOOD,        
  TACTYPE_PARENT,       
  TACTYPE_METABOLITE,   
  TACTYPE_VOI,          
  TACTYPE_REF,          
  TACTYPE_PIXEL,        
  TACTYPE_LAST
} tactype;

◆ tpcerror

enum tpcerror

Error codes. Error code tpcerror, same as error message position in the table (table can not be accessed directly outside the c file).

See also: tpcerrormsg, statusmsg.c, TPCSTATUS

Enumerator
TPCERROR_OK	No error.
TPCERROR_FAIL	General error.
TPCERROR_OUT_OF_MEMORY	Cannot allocate memory.
TPCERROR_NO_FILE	File does not exist.
TPCERROR_FILE_EXISTS	File exists.
TPCERROR_INVALID_FILENAME	Invalid file name.
TPCERROR_CANNOT_OPEN	Cannot open file.
TPCERROR_CANNOT_READ	Cannot read file.
TPCERROR_CANNOT_WRITE	Cannot write file.
TPCERROR_CANNOT_DELETE	Cannot delete.
TPCERROR_INVALID_FORMAT	Invalid file format.
TPCERROR_UNSUPPORTED	Unsupported file type.
TPCERROR_NO_DATA	File contains no data.
TPCERROR_DUPLICATE_DATA	File contains duplicate data.
TPCERROR_NO_WEIGHTS	File contains no weights.
TPCERROR_NO_DATETIME	File contains no date or time.
TPCERROR_MISSING_DATA	File contains missing values.
TPCERROR_TOO_FEW	File contains too few samples.
TPCERROR_TOO_BIG	File is too big.
TPCERROR_OVERLAPPING_DATA	Overlapping data.
TPCERROR_LARGE_GAP	Large gap in data.
TPCERROR_NO_X	No sample times.
TPCERROR_INVALID_X	Invalid sample time.
TPCERROR_INVALID_XRANGE	Invalid sample time range.
TPCERROR_INCOMPATIBLE_UNIT	Incompatible units.
TPCERROR_INCOMPATIBLE_DATA	Incompatible data.
TPCERROR_NO_REFERENCE	Reference not found.
TPCERROR_NO_KEY	Key not found.
TPCERROR_NO_VALUE	Value not found.
TPCERROR_INVALID_VALUE	Invalid value.
TPCERROR_INVALID_SEPARATOR	Invalid field delimiter.
TPCERROR_INVALID_PARNR	Invalid number of parameters.
TPCERROR_INVALID_HEADER	Invalid header contents.
TPCERROR_MISSING_HEADER	Missing header.
TPCERROR_UNKNOWN_ISOTOPE	Unknown isotope.
TPCERROR_UNKNOWN_UNIT	Unknown data unit.
TPCERROR_NO_SOLUTION	No solution.
TPCERROR_BAD_FIT	Fitting not successful.
TPCERROR_LAST	marking the end of list

Definition at line 189 of file tpcextensions.h.

             {
// tpcerrormsg[] in statusmsg.c must be edited accordingly
  TPCERROR_OK,                 
  TPCERROR_FAIL,               
  TPCERROR_OUT_OF_MEMORY,      
  TPCERROR_NO_FILE,            
  TPCERROR_FILE_EXISTS,        
  TPCERROR_INVALID_FILENAME,   
  TPCERROR_CANNOT_OPEN,        
  TPCERROR_CANNOT_READ,        
  TPCERROR_CANNOT_WRITE,       
  TPCERROR_CANNOT_DELETE,      
  TPCERROR_INVALID_FORMAT,     
  TPCERROR_UNSUPPORTED,        
  TPCERROR_NO_DATA,            
  TPCERROR_DUPLICATE_DATA,     
  TPCERROR_NO_WEIGHTS,         
  TPCERROR_NO_DATETIME,        
  TPCERROR_MISSING_DATA,       
  TPCERROR_TOO_FEW,            
  TPCERROR_TOO_BIG,            
  TPCERROR_OVERLAPPING_DATA,   
  TPCERROR_LARGE_GAP,          
  TPCERROR_NO_X,               
  TPCERROR_INVALID_X,          
  TPCERROR_INVALID_XRANGE,     
  TPCERROR_INCOMPATIBLE_UNIT,  
  TPCERROR_INCOMPATIBLE_DATA,  
  TPCERROR_NO_REFERENCE,       
  TPCERROR_NO_KEY,             
  TPCERROR_NO_VALUE,           
  TPCERROR_INVALID_VALUE,      
  TPCERROR_INVALID_SEPARATOR,  
  TPCERROR_INVALID_PARNR,      
  TPCERROR_INVALID_HEADER,     
  TPCERROR_MISSING_HEADER,     
  TPCERROR_UNKNOWN_ISOTOPE,    
  TPCERROR_UNKNOWN_UNIT,       
  TPCERROR_NO_SOLUTION,        
  TPCERROR_BAD_FIT,            
  TPCERROR_LAST                
} tpcerror;

◆ unit

enum unit

unit code, same as unit position in the table (table can not be accessed directly outside the c file).

See also: units.c unitName, unitIdentify, unitConversionFactor

Enumerator
UNIT_UNKNOWN	Unknown unit.
UNIT_UNITLESS	Unitless.
UNIT_HU	Hounsfield Unit.
UNIT_MSEC	milliseconds
UNIT_SEC	seconds
UNIT_MIN	minutes
UNIT_HOUR	hours
UNIT_DAY	days
UNIT_MONTH	months
UNIT_YEAR	years
UNIT_UM	micrometer
UNIT_MM	millimeter
UNIT_CM	centimeter
UNIT_M	meter
UNIT_UL	microliter
UNIT_ML	milliliter
UNIT_DL	desiliter
UNIT_L	liter
UNIT_UG	microgram
UNIT_MG	milligram
UNIT_G	gram
UNIT_100G	100 grams
UNIT_KG	kilogram
UNIT_PMOL	picomole
UNIT_NMOL	nanomole
UNIT_UMOL	micromole
UNIT_MMOL	millimole
UNIT_MOL	mole
UNIT_COUNTS	counts
UNIT_KCOUNTS	kilocounts
UNIT_BQ	Becquerel.
UNIT_KBQ	kiloBecquerel
UNIT_MBQ	megaBecquerel
UNIT_GBQ	gigaBecquerel
UNIT_NCI	nanoCurie
UNIT_UCI	microCurie
UNIT_MCI	milliCurie
UNIT_CI	Curie.
UNIT_CPS	counts/s
UNIT_KCPS	kilocounts/s
UNIT_CPM	counts/min
UNIT_KCPM	kilocounts/min
UNIT_BQ_PER_ML	Bq/mL.
UNIT_KBQ_PER_ML	kBq/mL
UNIT_MBQ_PER_ML	MBq/mL.
UNIT_NCI_PER_ML	nCi/mL
UNIT_UCI_PER_ML	uCi/mL
UNIT_BQ_PER_G	Bq/g.
UNIT_KBQ_PER_G	kBq/g
UNIT_MBQ_PER_G	MBq/g.
UNIT_NCI_PER_G	nCi/g
UNIT_UCI_PER_G	uCi/g
UNIT_PMOL_PER_ML	pmol/mL
UNIT_NMOL_PER_ML	nmol/mL
UNIT_UMOL_PER_ML	umol/mL
UNIT_MMOL_PER_ML	mmol/mL
UNIT_PMOL_PER_L	pmol/L
UNIT_NMOL_PER_L	nmol/L
UNIT_UMOL_PER_L	umol/L
UNIT_MMOL_PER_L	mmol/L
UNIT_MOL_PER_L	mol/L
UNIT_SEC_KBQ_PER_ML	s*kBq/mL
UNIT_MIN_KBQ_PER_ML	min*kBq/mL
UNIT_SEC_BQ_PER_ML	s*Bq/mL
UNIT_MIN_BQ_PER_ML	min*Bq/mL
UNIT_PERCENTAGE	Percentage (%)
UNIT_ML_PER_ML	mL/mL
UNIT_ML_PER_DL	mL/dL
UNIT_G_PER_ML	g/mL
UNIT_ML_PER_G	mL/g
UNIT_PER_SEC	1/s
UNIT_PER_MIN	1/min
UNIT_PER_HOUR	1/h
UNIT_ML_PER_ML_SEC	mL/(mL*sec)
UNIT_ML_PER_ML_MIN	mL/(mL*min)
UNIT_ML_PER_DL_MIN	mL/(dL*min)
UNIT_ML_PER_G_MIN	mL/(g*min)
UNIT_ML_PER_100G_MIN	mL/(100g*min)
UNIT_UMOL_PER_ML_MIN	umol/(mL*min)
UNIT_UMOL_PER_DL_MIN	umol/(dL*min)
UNIT_UMOL_PER_G_MIN	umol/(g*min)
UNIT_UMOL_PER_100G_MIN	umol/(100g*min)
UNIT_MMOL_PER_ML_MIN	mmol/(mL*min)
UNIT_MMOL_PER_DL_MIN	mmol/(dL*min)
UNIT_MMOL_PER_G_MIN	mmol/(g*min)
UNIT_MMOL_PER_100G_MIN	mmol/(100g*min)
UNIT_MG_PER_100G_MIN	mg/(100g*min)
UNIT_MG_PER_DL_MIN	mg/(dL*min)
UNIT_PID	Percent of injected dose.
UNIT_PID_PER_G	Percent of injected dose / g.
UNIT_PID_PER_ML	Percent of injected dose / mL.
UNIT_PID_PER_KG	Percent of injected dose / kg.
UNIT_PID_PER_L	Percent of injected dose / L.
UNIT_MBQ_PER_NMOL	MBq/nmol (specific activity)
UNIT_GBQ_PER_NMOL	GBq/nmol (specific activity)
UNIT_LAST	end of list

Definition at line 84 of file tpcextensions.h.

             {
  UNIT_UNKNOWN,              
  UNIT_UNITLESS,             
  UNIT_HU,                   
  UNIT_MSEC,                 
  UNIT_SEC,                  
  UNIT_MIN,                  
  UNIT_HOUR,                 
  UNIT_DAY,                  
  UNIT_MONTH,                
  UNIT_YEAR,                 
  UNIT_UM,                   
  UNIT_MM,                   
  UNIT_CM,                   
  UNIT_M,                    
  UNIT_UL,                   
  UNIT_ML,                   
  UNIT_DL,                   
  UNIT_L,                    
  UNIT_UG,                   
  UNIT_MG,                   
  UNIT_G,                    
  UNIT_100G,                 
  UNIT_KG,                   
  UNIT_PMOL,                 
  UNIT_NMOL,                 
  UNIT_UMOL,                 
  UNIT_MMOL,                 
  UNIT_MOL,                  
  UNIT_COUNTS,               
  UNIT_KCOUNTS,              
  UNIT_BQ,                   
  UNIT_KBQ,                  
  UNIT_MBQ,                  
  UNIT_GBQ,                  
  UNIT_NCI,                  
  UNIT_UCI,                  
  UNIT_MCI,                  
  UNIT_CI,                   
  UNIT_CPS,                  
  UNIT_KCPS,                 
  UNIT_CPM,                  
  UNIT_KCPM,                 
  UNIT_BQ_PER_ML,            
  UNIT_KBQ_PER_ML,           
  UNIT_MBQ_PER_ML,           
  UNIT_NCI_PER_ML,           
  UNIT_UCI_PER_ML,           
  UNIT_BQ_PER_G,             
  UNIT_KBQ_PER_G,            
  UNIT_MBQ_PER_G,            
  UNIT_NCI_PER_G,            
  UNIT_UCI_PER_G,            
  UNIT_PMOL_PER_ML,          
  UNIT_NMOL_PER_ML,          
  UNIT_UMOL_PER_ML,          
  UNIT_MMOL_PER_ML,          
  UNIT_PMOL_PER_L,           
  UNIT_NMOL_PER_L,           
  UNIT_UMOL_PER_L,           
  UNIT_MMOL_PER_L,           
  UNIT_MOL_PER_L,            
  UNIT_SEC_KBQ_PER_ML,       
  UNIT_MIN_KBQ_PER_ML,       
  UNIT_SEC_BQ_PER_ML,        
  UNIT_MIN_BQ_PER_ML,        
  UNIT_PERCENTAGE,           
  UNIT_ML_PER_ML,            
  UNIT_ML_PER_DL,            
  UNIT_G_PER_ML,             
  UNIT_ML_PER_G,             
  UNIT_PER_SEC,              
  UNIT_PER_MIN,              
  UNIT_PER_HOUR,             
  UNIT_ML_PER_ML_SEC,        
  UNIT_ML_PER_ML_MIN,        
  UNIT_ML_PER_DL_MIN,        
  UNIT_ML_PER_G_MIN,         
  UNIT_ML_PER_100G_MIN,      
  UNIT_UMOL_PER_ML_MIN,      
  UNIT_UMOL_PER_DL_MIN,      
  UNIT_UMOL_PER_G_MIN,       
  UNIT_UMOL_PER_100G_MIN,    
  UNIT_MMOL_PER_ML_MIN,      
  UNIT_MMOL_PER_DL_MIN,      
  UNIT_MMOL_PER_G_MIN,       
  UNIT_MMOL_PER_100G_MIN,    
  UNIT_MG_PER_100G_MIN,      
  UNIT_MG_PER_DL_MIN,        
  UNIT_PID,                  
  UNIT_PID_PER_G,            
  UNIT_PID_PER_ML,           
  UNIT_PID_PER_KG,           
  UNIT_PID_PER_L,            
  UNIT_MBQ_PER_NMOL,         
  UNIT_GBQ_PER_NMOL,         
  UNIT_LAST                  
} unit;

◆ weights

enum weights

Is data weighted, or are weight factors available with data?

Enumerator
WEIGHTING_UNKNOWN	Not known; usually assumed that not weighted.
WEIGHTING_OFF	Not weighted or weights not available (weights for all included samples are 1.0).
WEIGHTING_ON_GENERAL	Weighted or weights are available, but not specified.
WEIGHTING_ON_COUNTS	Weights based on counts (Mazoyer et al, 1986).
WEIGHTING_ON_F	Weights based on sample frequency or frame length.
WEIGHTING_ON_FD	Weights based on decay and sample frequency or frame length (Thiele et al, 2008).

Definition at line 52 of file tpcextensions.h.

             {
  WEIGHTING_UNKNOWN,    
  WEIGHTING_OFF,        
  WEIGHTING_ON_GENERAL, 
  WEIGHTING_ON_COUNTS,  
  WEIGHTING_ON_F,       
  WEIGHTING_ON_FD,      
  // possibility to add more weighting schemes
  WEIGHTING_LAST
} weights;

Function Documentation

◆ 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

See also: asciiFileRead, asciiFileSize, strTrimLeft, strstrNoQuotation

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 149 of file readasciifile.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 iftPutFromString(), iftPutFromStringWithSpaceSeparator(), and tacFormatDetermine().

◆ asciiFileRead()

char * asciiFileRead	(	FILE *	fp,
		char *	data,
		size_t	maxlen )

extern

Read at most maxlen-1 characters from given file pointer to ASCII text file.

Author: Vesa Oikonen

Returns: Pointer to the string, or NULL in case of an error.

See also: asciiFileSize, asciiCommentLine, fileExist, filenameGetExtension

Parameters

fp	File pointer; rewind() is NOT called here to allow reading data in chunks.
data	Pointer to character array in which the data is written; NULL, if memory is to allocated here (in that case remember to free it after use). Terminal zero is added.
maxlen	Nr of characters to read, and the size of data array (including terminal zero).

Definition at line 117 of file readasciifile.c.

  {
  /* Check the input */
  if(fp==NULL || maxlen<2) return NULL;
  /* Allocate memory if needed */
  char *cptr;
  if(data!=NULL) cptr=data; else {
    cptr=malloc(maxlen*sizeof(char)); if(cptr==NULL) return NULL;
  }
  size_t i=0; int c;
  while((c=fgetc(fp))!=EOF && i<maxlen-1) cptr[i++]=(char)c;
  cptr[i]=(char)0; if(i==0) {if(data==NULL) free(cptr); cptr=NULL;}
  return cptr;
}

Referenced by csvRead(), iftRead(), and tacFormatDetermine().

◆ asciiFileSize()

size_t asciiFileSize	(	FILE *	fp,
		int *	nonprintable )

extern

Get the size of ASCII text file, or size of initial ASCII part of file containing binary data after the text part (in certain image formats).

One NUL character is not counted as binary data, but two already is; this, because equipment failures sometimes lead to data ending with NUL.

Author: Vesa Oikonen

Returns: The number of characters in the initial ASCII part of file. This does not include the terminal zero that would be added to a string.

See also: asciiFileRead, asciiCommentLine, fileExist

Parameters

fp	File pointer; rewind() is called before returning.
nonprintable	Value will be set to 1 if reading was stopped at a non-printable character meaning that file may contain a binary part, or to 0 if file contains only printable characters and size of all contents was counted. Enter NULL, if not needed.

Definition at line 26 of file readasciifile.c.

  {
  if(nonprintable!=NULL) *nonprintable=0;
  if(fp==NULL) return(0);
  size_t n, bn;
  int c;
  n=bn=0;
  while((c=fgetc(fp))!=EOF && n<SIZE_MAX-1 && bn<2) {
    if(iscntrl(c) && c!=13 && c!=10 && c!=9 && c!=169) {
      bn++;
      if(c!=0) {
        if(nonprintable!=NULL) *nonprintable=1; 
        break;
      }
    }
    n++;
  }
  rewind(fp);
  if(nonprintable!=NULL && bn>1) *nonprintable=1; 
  return n;
}

Referenced by csvRead(), iftRead(), tacFormatDetermine(), and tacRead().

◆ atofCheck()

int atofCheck	(	const char *	s,
		double *	v )

extern

Verifies that argument string at least seems like a valid number.

Both decimal point and comma are accepted. Optional result double value is set to NaN if string was not valid value.

See also: atofVerified, atoiCheck, doubleGetWithUnit, atofList, strCleanSpaces, strIsValidNumber

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

Parameters

s	String which is converted to a double; string must not contain any space characters.
v	Pointer to the double float; enter NULL, if not needed.

Definition at line 94 of file decpoint.c.

  {
  if(v!=NULL) *v=nan(""); 
  if(s==NULL) return 1;
  if(!strIsValidNumber(s)) return 1;
  char *p; p=strchr(s, ','); if(p==NULL) {*v=atof(s); return 0;}
  char *s2=strdup(s); p=strchr(s2, ','); *p='.';
  *v=atof(s2); free(s2); return 0;
}

Referenced by doubleGetWithUnit(), floatGetWithUnit(), iftGetDoubleValue(), tacReadCarimasTxt(), and tacReadInveonCSV().

◆ atofList()

int atofList	(	const char *	s1,
		const char *	s2,
		double *	x,
		int	maxn )

extern

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

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

Author: Vesa Oikonen

See also: strTokenNr, strTokenNCpy, intlistAddFromString, atofCheck, tpcYesNo, doubleCopy, doubleSpanPositives

Parameters

s1	Pointer to string from which double values are read.
s2	String containing character delimiters
x	Pointer to a pre-allocated array of doubles.
maxn	Size of double list; obligatory. If string contains more values, only maxn values will be read.

Definition at line 150 of file decpoint.c.

  {
  if(x==NULL || maxn<=0) return(-1);
  for(int i=0; i<maxn; i++) x[i]=nan("");
  if(s1==NULL || s2==NULL) return(0);
  
  /* Get the nr of tokens */
  int i, n;
  n=strTokenNr((char*)s1, s2); if(n<1) return(0);
  /* Read the values */
  char tmp[256];
  double v; 
  for(i=0; i<n && i<maxn; i++) {
    if(strTokenNCpy(s1, s2, 1+i, tmp, 256)<1) return(-2);
    v=atofVerified(tmp); if(isnan(v)) return(-3);
    x[i]=v;
  }
  return(i);
}

◆ atofVerified()

double atofVerified ( const char * s )

extern

Version of atof() which verifies that argument string at least seems like a valid number.

Both decimal point and comma are accepted. Result value is set to NaN if string was not valid value.

See also: atofCheck, doubleGetWithUnit, atoiCheck, strIsValidNumber, tpcYesNo

Returns: converted double value, or NaN in case of an error.

Parameters

s	String which is converted to a double; string must not contain any space characters.

Definition at line 75 of file decpoint.c.

  {
  if(s==NULL || !strIsValidNumber(s)) return nan("");
  char *p; p=strchr(s, ','); if(p==NULL) return atof(s);
  char *s2=strdup(s); p=strchr(s2, ','); *p='.';
  double t=atof(s2); free(s2); return t;
}

Referenced by atofList(), dcmImgIsotope(), iftGetDoubleWithUnit(), parReadCSV(), parReadFIT(), parReadRES(), tacRead4DM(), tacReadCSV(), tacReadDFT(), tacReadHRPLUSHC(), tacReadHRRTHC(), tacReadMat(), tacReadPMOD(), tacReadQView(), tacReadSIF(), and tacReadSimple().

◆ atoiCheck()

int atoiCheck	(	const char *	s,
		int *	v )

extern

Verifies that argument string at least seems like a valid integer. '\0' must follow the integer part in the string. Optional result integer value is set to 0 if string was not valid value.

See also: atofCheck, intlistAddFromString, strIsValidNumber

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

Parameters

s	String which is converted to an integer; string must not contain any space characters.
v	Pointer to the integer; enter NULL, if not needed.

Definition at line 25 of file intutil.c.

  {
  if(v!=NULL) *v=0; 
  if(s==NULL || strnlen(s, 2)<1) return(1);
  errno=0; char *tail;
  *v=strtol(s, &tail, 10);
  if(errno) {*v=0; return(2);}
  if(*tail) {*v=0; return(3);}
  return 0;
}

Referenced by intlistAddFromString(), and tacReadQView().

◆ ctime_r_int()

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

extern

Convert calendar 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.

See also: time_to_tm, strDateTimeValid, strDateTimeRead

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 119 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 parWriteRES().

◆ doubleAbsMaxIndex()

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

extern

Find the maximum absolute value in given double array.

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

See also: doubleMaxIndex, doubleGTIndex, doubleGEIndex, doubleRange

Parameters

a	Pointer to double array.
n	Length of array.

Definition at line 377 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) || fabs(a[i])>mv) {mv=fabs(a[i]); mi=i;}
  return(mi);
}

◆ doubleAbsMinIndex()

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

extern

Find the minimum absolute value in given double array.

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

See also: doubleAbsMaxIndex, doubleMinIndex, doubleGTIndex, doubleGEIndex

Parameters

a	Pointer to double array.
n	Length of array.

Definition at line 416 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) || fabs(a[i])<mv) {mv=fabs(a[i]); mi=i;}
  return(mi);
}

Referenced by spectralBFExtract().

◆ doubleArrayMatch()

int doubleArrayMatch	(	const double *	a1,
		const double *	a2,
		const unsigned int	n,
		const double	lim )

extern

Verifies that given two double arrays have the same values 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: doubleMatch, doubleMatchRel, doubleMachEps, statMeanSD

Parameters

a1	First array.
a2	Second array.
n	Length of arrays.
lim	Limit for absolute difference (if <0 then test will fail every time).

Definition at line 51 of file doubleutil.c.

  {
  if(a1==NULL || a2==NULL || n<1) return 0;
  for(unsigned int i=0; i<n; i++) if(!doubleMatch(a1[i], a2[i], lim)) return 0;
  return 1;
}

◆ doubleCopy()

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

extern

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

See also: floatCopy, doubleCopyFinite

Parameters

t	Target array.
s	Source array.
n	Length of arrays.

Definition at line 117 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];
}

Referenced by nloptIATGO(), nloptMPSO(), and nloptSimplexARRS().

◆ doubleCopyFinite()

unsigned int doubleCopyFinite	(	double *	t,
		double *	s,
		const unsigned int	n )

extern

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

See also: doubleCopy

Returns: The number of finite values that were copied.

Parameters

t	Target array.
s	Source array.
n	Length of arrays.

Definition at line 136 of file doubleutil.c.

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

◆ 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: floatCSpanPositives, doubleSpanPositives, atofList

Parameters

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

Definition at line 320 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);
}

Referenced by spectralBFExtract(), and spectralBFNr().

◆ doubleGEIndex()

unsigned int doubleGEIndex	(	double *	a,
		const unsigned int	n,
		double	lim )

extern

Get the index of first array value greater or equal than given limit.

Returns: The index, or n, if not found.

See also: doubleGTIndex, doubleMaxIndex, doubleMinIndex, fitLine, statMeanSD

Parameters

a	Pointer to an array of values.
n	The number of samples (length of a[]).
lim	The limit value.

Definition at line 435 of file doubleutil.c.

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

◆ doubleGetWithUnit()

int doubleGetWithUnit	(	const char *	s,
		double *	v,
		int *	u )

extern

Reads numerical value and its unit from argument string.

Both decimal point and comma are accepted. Optional result double value is set to NaN if string was not valid value. Optional unit is set to UNIT_UNKNOWN, if not valid or not found.

Returns: 0 if numerical value successfully read, and 1 in case of an error.

See also: floatGetWithUnit, atofVerified, atofCheck

Parameters

s	String from where double and unit is read; string must not contain any extra space characters.
v	Pointer to the double float; enter NULL, if not needed.
u	Pointer to int for unit code; enter NULL, if not needed.

Definition at line 272 of file doubleutil.c.

  {
  if(v!=NULL) *v=nan(""); 
  if(u!=NULL) *u=UNIT_UNKNOWN; 
  if(s==NULL) return 1;
  int n=strTokenNr(s, " \t"); if(n==0 || n>2) return 1;
  if(n==1) return(atofCheck(s, v)); // no unit
  /* separate number and unit */
  char buf[64];
  n=strTokenNCpy(s, " \t", 1, buf, 64);
  if(atofCheck(buf, v)) return 1;
  if(u==NULL) return 0;
  n=strTokenNCpy(s, " \t", 2, buf, 64);
  *u=unitIdentify(buf);
  return 0;
}

Referenced by iftCheckKeyValue(), and iftCheckKeyValues().

◆ doubleGTIndex()

unsigned int doubleGTIndex	(	double *	a,
		const unsigned int	n,
		double	lim )

extern

Get the index of first array value greater than given limit.

Returns: The index, or n, if not found.

See also: doubleGEIndex, doubleMaxIndex, doubleMinIndex, fitLine, statMeanSD, doubleMatch

Parameters

a	Pointer to an array of values.
n	The number of samples (length of a[]).
lim	The limit value.

Definition at line 452 of file doubleutil.c.

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

◆ 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: floatMachEps, doubleMatchRel, doubleMatch, factorial

Definition at line 105 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(), nlopt1D(), nloptFixedNr(), nloptLimitFixedNr(), and nloptPowellBrent().

◆ 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: floatMatch, doubleMatchRel, doubleMachEps, atofCheck, statMeanSD

Parameters

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

Definition at line 27 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;
}

Referenced by doubleArrayMatch(), parCompareParameters(), parCompareWSS(), parIsFitRange(), tacCheckX1X2X(), tacCompareConc(), tacCompareTimes(), tacCompareWeights(), and tacReadQView().

◆ 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: floatMatchRel, 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 77 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;
}

Referenced by parCompareParameters(), parCompareWSS(), tacCompareConc(), tacCompareTimes(), and tacCompareWeights().

◆ 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.

See also: floatMaxIndex, doubleSpanPositives, doubleCSpanPositives, doubleGTIndex, doubleRange

Parameters

a	Pointer to double array.
n	Length of array.

Definition at line 357 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: floatMean, doubleWMean, doubleSum, doubleMaxIndex, statMeanSD, statMedian

Parameters

a	Pointer to double array.
n	Length of array.

Definition at line 224 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);
}

◆ doubleMinIndex()

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

extern

Find the minimum value in given double array.

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

See also: doubleMaxIndex, doubleAbsMinIndex, doubleRange

Parameters

a	Pointer to double array.
n	Length of array.

Definition at line 396 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);
}

◆ doubleNaNs()

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

extern

Check the number of NaNs in a double array.

See also: doubleNonzeroes, doubleCopyFinite, doubleRange, doubleMean

Returns: The number of NaNs.

Parameters

a	Pointer to double array.
n	Length of array.

Definition at line 156 of file doubleutil.c.

  {
  if(a==NULL || n<1) return(0);
  unsigned int i, nn=0;
  for(i=0; i<n; i++) if(isnan(a[i])) nn++;
  return(nn);
}

◆ doubleNonzeroes()

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

extern

Find the number of non-zero values in given double array. Negative values are counted, but NaNs are not.

Returns: The number of non-zero values in array; 0 is returned also in case of errors.

See also: floatNonzeroes, doubleSpanPositives, doubleNaNs

Parameters

a	Pointer to double array.
n	Length of array.

Definition at line 339 of file doubleutil.c.

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

◆ doubleRange()

unsigned int doubleRange	(	double *	a,
		const unsigned int	n,
		double *	amin,
		double *	amax )

extern

Get the range of finite values in a double array.

See also: doubleMaxIndex, doubleMinIndex, doubleNonzeroes, doubleNaNs, doubleMean, statMedian

Returns: The number of finite values in the array, or 0 in case of an error.

Parameters

a	Pointer to double array.
n	Length of array.
amin	Pointer to the double float where minimum value is stored; enter NULL, if not needed.
amax	Pointer to the double float where maximum value is stored; enter NULL, if not needed.

Definition at line 174 of file doubleutil.c.

  {
  if(amin!=NULL) *amin=nan("");
  if(amax!=NULL) *amax=nan("");
  if(a==NULL || n<1) return(0);
 
  double mi=nan(""), ma=nan("");
  unsigned int i, nn=0;
  for(i=0; i<n; i++) if(isfinite(a[i])) {
    nn++;
    if(!isfinite(mi) || mi>a[i]) mi=a[i];
    if(!isfinite(ma) || ma<a[i]) ma=a[i];
  }
  if(amin!=NULL) *amin=mi;
  if(amax!=NULL) *amax=ma;
  return(nn);
}

◆ 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 non-positive value in the given array.

See also: floatSpanPositives, doubleCSpanPositives, atofList, doubleNonzeroes

Parameters

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

Definition at line 302 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);
}

Referenced by spectralBFExtract(), and spectralBFNr().

◆ 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: floatSum, doubleMean, doubleMaxIndex

Parameters

a	Pointer to double array.
n	Length of array.

Definition at line 206 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);
}

Referenced by spectralBFExtract().

◆ doubleWMean()

double doubleWMean	(	double *	a,
		double *	w,
		const unsigned int	n )

extern

Calculate the weighted mean of values in given double array.

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

See also: doubleMean, floatMean, doubleSum, doubleMaxIndex, statMeanSD

Parameters

a	Pointer to double array to calculated weighted mean from.
w	Pointer to double array containing the weights. Negative or NaN weight is interpreted as zero weight.
n	Length of arrays.

Definition at line 244 of file doubleutil.c.

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

Referenced by spectralBFExtract().

◆ endianLittle()

int endianLittle ( )

extern

Check whether current platform uses little endian byte order.

See also: swap, swabip, swawbip, swawip

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

Definition at line 53 of file endian.c.

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

Referenced by anaReadHeader(), dcmFileWrite(), dcmitemGetInt(), dcmitemGetReal(), dcmReadFileTag(), dcmReadFileVL(), dcmReadFileVRVL(), dcmValueString(), dcmWriteFileTag(), dcmWriteFileVRVL(), ecatMListRead(), ecatReadMainheader(), ecatWriteMainheader(), imgReadNifti(), imgSetNiftiHeader(), niftiReadHeader(), and niftiWriteHeader().

◆ endianness()

int endianness ( )

extern

Detect endianness at run-time. From https://sourceforge.net/p/predef/wiki/Endianness/

Returns: Returns ENDIAN_UNKNOWN, ENDIAN_BIG, ENDIAN_LITTLE, ENDIAN_BIG_WORD, or ENDIAN_LITTLE_WORD.

Definition at line 26 of file endian.c.

{
  union {
    uint32_t value;
    uint8_t data[sizeof(uint32_t)];
  } number;
 
  number.data[0] = 0x00;
  number.data[1] = 0x01;
  number.data[2] = 0x02;
  number.data[3] = 0x03;
 
  switch (number.value) {
    case UINT32_C(0x00010203): return ENDIAN_BIG;
    case UINT32_C(0x03020100): return ENDIAN_LITTLE;
    case UINT32_C(0x02030001): return ENDIAN_BIG_WORD;
    case UINT32_C(0x01000302): return ENDIAN_LITTLE_WORD;
    default:                   return ENDIAN_UNKNOWN;
  }
}

Referenced by endianLittle().

◆ errorMsg()

char * errorMsg ( tpcerror e )

extern

Return pointer to error message specified with the error code.

Returns: pointer to the error message string.

See also: tpcerror, statusPrint

Author: Vesa Oikonen

Parameters

e	Error code

See also: tpcerror.

Definition at line 68 of file statusmsg.c.

  {
  if(e<TPCERROR_OK || e>=TPCERROR_LAST) return NULL;
  return tpcerrormsg[e];
}

Referenced by imgRead(), and nloptIATGO().

◆ factorial()

unsigned int factorial ( unsigned int n )

extern

Calculate factorial of given number.

See also: lfactorial, igam, igamc

Returns: Returns factorial n!, or zero if factorial would cause wrap-around.

Parameters

n	Integer n, from which the factorial is calculated.

Note: Wrap-around will occur if n>12.

Definition at line 46 of file intutil.c.

  {
  if(n<1) return(1);
  if(n>12) return(0);
  return(n*factorial(n-1));
}

Referenced by factorial().

◆ filenameGetExtension()

char * filenameGetExtension ( const char * s )

extern

Get the last extension of a file name.

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

See also: filenameGetExtensions, filenameRmExtension, strTrimLeft

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

Parameters

s	Pointer to string; string is not edited here.

Definition at line 178 of file filename.c.

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

Referenced by abssIdFromFName(), imgWriteNifti(), micropetExists(), niftiBasename(), and tacRead().

◆ filenameGetExtensions()

char * filenameGetExtensions ( const char * s )

extern

Get all extensions of a file name.

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

See also: filenameGetExtension, filenameRmExtensions, strTrimLeft

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

Parameters

s	Pointer to string; string is not edited here.

Definition at line 203 of file filename.c.

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

Referenced by imgFormatFromFName(), and parFormatFromExtension().

◆ filenameRmExtension()

int filenameRmExtension ( char * s )

extern

Remove the last extension from file name.

See also: filenameGetExtension, filenameRmExtensions, strTrimLeft

Author: Vesa Oikonen

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

Parameters

s	Pointer to string.

Definition at line 71 of file filename.c.

  {
  if(s==NULL || strnlen(s, 1)<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(), micropetExists(), tpcHtmlUsage(), tpcPrintBuild(), tpcPrintUsage(), and tpcProgramName().

◆ filenameRmExtensions()

void filenameRmExtensions ( char * s )

extern

Remove all extensions from file name.

See also: filenameGetExtensions, filenameRmExtension, strdelstr, strTrimLeft

Author: Vesa Oikonen

Parameters

s	Pointer to string.

Definition at line 89 of file filename.c.

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

Referenced by studynrFromFilename().

◆ filenameRmFile()

void filenameRmFile ( char * s )

extern

Remove file name from file name string including path, leaving only path without trailing '/'. If there is no path, then result is ".". For root path the result string is empty. See also POSIX function dirname() in libgen.h.

See also: filenameRmPath, filenameRmExtensions, pathExist

Author: Vesa Oikonen

Parameters

s	Pointer to string.

Definition at line 44 of file filename.c.

  {
  if(s==NULL) return;
  unsigned int len=strnlen(s, FILENAME_MAX);
  if(len==0) return;
  //if(len==1) {s[0]='.'; return;}
  char *cptr;
  cptr=strrchr(s, '/'); if(cptr==NULL) cptr=strrchr(s, '\\');
  if(cptr!=NULL) {
    *cptr=(char)0;
    if(strcmp(s, "/")==0 || strcmp(s, "\\")==0) s[0]=(char)0;
    return;
  }
  if(strcmp(s, "..")==0) return;
  s[0]='.'; s[1]=(char)0;
  return;
}

Referenced by dcmFileList(), and imgWriteNifti().

◆ filenameRmPath()

void filenameRmPath ( char * s )

extern

Remove path from file name. See also POSIX function basename() in libgen.h.

See also: filenameRmFile, filenameRmExtensions, strTrimleft

Author: Vesa Oikonen

Parameters

s	Pointer to string.

Definition at line 20 of file filename.c.

  {
  if(s==NULL || strnlen(s, 1)<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 imgWriteNifti(), niftiCreateFNames(), studynrFromFilename(), tpcHtmlUsage(), tpcPrintBuild(), tpcPrintUsage(), and tpcProgramName().

◆ floatCopy()

void floatCopy	(	float *	t,
		float *	s,
		const unsigned int	n )

extern

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

See also: doubleCopy

Parameters

t	Target array.
s	Source array.
n	Length of arrays.

Definition at line 94 of file floatutil.c.

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

◆ floatCSpanPositives()

int floatCSpanPositives	(	float *	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: doubleCSpanPositives, floatSpanPositives, atofList

Parameters

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

Definition at line 228 of file floatutil.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);
}

◆ floatGetWithUnit()

int floatGetWithUnit	(	const char *	s,
		float *	v,
		int *	u )

extern

Reads numerical value and its unit from argument string.

Both decimal point and comma are accepted. Optional result float value is set to NaN if string was not valid value. Optional unit is set to UNIT_UNKNOWN, if not valid or not found.

Returns: 0 if numerical value successfully read, and 1 in case of an error.

See also: doubleGetWithUnit, atofVerified, atofCheck

Parameters

s	String from where float and unit is read; string must not contain any extra space characters.
v	Pointer to the float; enter NULL, if not needed.
u	Pointer to int for unit code; enter NULL, if not needed.

Definition at line 174 of file floatutil.c.

  {
  if(v!=NULL) *v=nanf(""); 
  if(u!=NULL) *u=UNIT_UNKNOWN; 
  if(s==NULL) return 1;
  int n=strTokenNr(s, " \t"); if(n==0 || n>2) return 1;
  if(n==1) { // no unit
    double dv;
    if(atofCheck(s, &dv)) return 1;
    *v=(float)dv; return 0;
  }
  /* separate number and unit */
  char buf[64];
  n=strTokenNCpy(s, " \t", 1, buf, 64);
  double dv;
  if(atofCheck(buf, &dv)) return 1;
  *v=(float)dv;
  if(u==NULL) return 0;
  n=strTokenNCpy(s, " \t", 2, buf, 64);
  *u=unitIdentify(buf);
  return 0;
}

◆ floatMachEps()

float floatMachEps ( )

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 FLT_EPSILON and DBL_EPSILON in float.h.

Returns: Estimate of machine epsilon.

Author: Vesa Oikonen

See also: doubleMachEps, floatMatchRel, floatMatch

Definition at line 82 of file floatutil.c.

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

◆ floatMatch()

int floatMatch	(	const float	v1,
		const float	v2,
		const float	lim )

extern

Verifies that given two floats 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: doubleMatch, floatMatchRel, floatMachEps, atofCheck

Parameters

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

Definition at line 27 of file floatutil.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(fabsf(v1-v2)<=lim) return 1;
  return 0;
}

Referenced by imgCompareConc(), imgCompareTimes(), and imgReadDICOM().

◆ floatMatchRel()

int floatMatchRel	(	const float	v1,
		const float	v2,
		const float	lim )

extern

Verifies that given two floats 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: doubleMatchRel, floatMatch, floatMachEps

Parameters

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

Definition at line 54 of file floatutil.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;
  float mean;
  mean=0.5*(v1+v2); if(!isnormal(mean)) return 0;
  if(fabsf((v1-v2)/mean)<=lim) return 1;
  return 0;
}

Referenced by imgCompareConc(), and imgCompareTimes().

◆ floatMaxIndex()

unsigned int floatMaxIndex	(	float *	a,
		const unsigned int	n )

extern

Find the maximum value in given float array.

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

See also: doubleMaxIndex, floatSpanPositives, floatCSpanPositives

Parameters

a	Pointer to float array.
n	Length of array.

Definition at line 113 of file floatutil.c.

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

◆ floatMean()

float floatMean	(	float *	a,
		const unsigned int	n )

extern

Calculate the mean of values in given float array.

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

See also: doubleMean, floatSum, floatMaxIndex, fstatMeanSD

Parameters

a	Pointer to float array.
n	Length of array.

Definition at line 150 of file floatutil.c.

  {
  if(a==NULL || n<1) return(nanf(""));
  float 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(nanf(""));
  return(s/(float)ci);
}

◆ floatNonzeroes()

unsigned int floatNonzeroes	(	float *	a,
		const unsigned int	n )

extern

Find the number of non-zero values in given float array. Negative values are counted, but NaNs are not.

Returns: The number of non-zero values in array; 0 is returned also in case of errors.

See also: doubleNonzeroes, floatSpanPositives

Parameters

a	Pointer to float array.
n	Length of array.

Definition at line 247 of file floatutil.c.

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

Referenced by imgContents().

◆ floatSpanPositives()

int floatSpanPositives	(	float *	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: doubleSpanPositives, floatCSpanPositives, atofList, floatNonzeroes

Parameters

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

Definition at line 210 of file floatutil.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);
}

◆ floatSum()

float floatSum	(	float *	a,
		const unsigned int	n )

extern

Calculate the sum of values in given float array.

Returns: The sum of array values.

See also: doubleSum, floatMean, floatMaxIndex

Parameters

a	Pointer to float array.
n	Length of array.

Definition at line 132 of file floatutil.c.

  {
  float 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);
}

◆ fncasematch()

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

extern

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

See also: fnmatch, fileExist, pathExist

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

Parameters

fname	file name that is evaluated.
key	key string which may contain wild-cards '?' and '*'.

Definition at line 139 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 roinameMatch().

◆ fnmatch()

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

extern

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

See also: fncasematch, fileExist, pathExist, roinameMatch, strncpyClean

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

Parameters

fname	file name that is evaluated.
key	key string which may contain wild cards '?' and '*'.

Definition at line 104 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().

◆ 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.

See also: localtime_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 31 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(), tacReadGEMS(), and timegm().

◆ intlistAdd()

int intlistAdd	(	INTLIST *	l,
		const int	v,
		const int	ifnew )

extern

Add one integer to the integer array in INTLIST structure.

See also: intlistInit, intlistFree, intlistSort, intlistAddFromString, strncatInt

Author: Vesa Oikonen

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

Parameters

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

Definition at line 108 of file intutil.c.

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

Referenced by intlistAddFromString(), and intlistExpandFromString().

◆ intlistAddFromString()

int intlistAddFromString	(	const char *	s1,
		const char *	s2,
		INTLIST *	l,
		const int	ifnew )

extern

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

See also: strTokenNr, strTokenNCpy, intlistAdd, intlistInit, intlistFree, intlistSort, intlistExpandFromString.

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 INTLIST structure; previous contents are preserved.
ifnew	Add integer to the list only if it is new (0=no, 1=yes).

Definition at line 156 of file intutil.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(atoiCheck(tmp, &v)) return(-3);
    m=intlistAdd(l, v, ifnew); if(m<0) return(-4);
    j+=m;
  }
  return(j);
}

◆ intlistExpandFromString()

int intlistExpandFromString	(	const char *	s1,
		const char *	s2,
		INTLIST *	l,
		const int	ifnew )

extern

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

See also: strTokenNCpy, intlistAddFromString, intlistInit, intlistFree, intlistSort, strIsValidNumber

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 INTLIST struct; previous contents are preserved.
ifnew	Add integer to the list only if it is new (0=no, 1=yes).

Definition at line 190 of file intutil.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=intlistAdd(l, v, ifnew); if(m<0) return(-10);
      j+=m;
    }
  }
  return(j);
}

Referenced by parSelectParameters(), and parSelectTACs().

◆ intlistFree()

void intlistFree ( INTLIST * l )

extern

Free memory allocated for INTLIST. All data is cleared.

See also: intlistInit

Author: Vesa Oikonen

Parameters

l	Pointer to INTLIST structure

Definition at line 92 of file intutil.c.

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

Referenced by parSelectParameters(), and parSelectTACs().

◆ intlistInit()

void intlistInit ( INTLIST * l )

extern

Initiate the INTLIST structure before any use.

See also: intlistFree, intlistAdd, intlistSort

Author: Vesa Oikonen

Parameters

l	Pointer to INTLIST.

Definition at line 80 of file intutil.c.

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

Referenced by intlistFree(), parSelectParameters(), and parSelectTACs().

◆ intlistSort()

void intlistSort ( INTLIST * l )

extern

Sort the integer array in INTLIST structure.

See also: intlistInit, intlistFree, intlistAdd

Author: Vesa Oikonen

Parameters

l	Pointer to INTLIST structure.

Definition at line 137 of file intutil.c.

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

Referenced by parSelectParameters(), and parSelectTACs().

◆ 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).

See also: doubleMachEps, igam, igamc

Parameters

x	Parameter for the inverse complementary error function.

Definition at line 473 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));
}

◆ lfactorial()

unsigned long long int lfactorial ( unsigned long long int n )

extern

Calculate factorial of given number.

See also: factorial, igam, igamc

Returns: Returns factorial n!, or zero if factorial would cause wrap-around.

Parameters

n	Integer n, from which the factorial is calculated.

Note: Wrap-around will occur if n>20.

Definition at line 63 of file intutil.c.

  {
  if(n<1) return(1);
  if(n>20) return(0);
  return(n*lfactorial(n-1));
}

Referenced by lfactorial(), mfEvalInt(), and mfEvalY().

◆ 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 54 of file datetime.c.

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

Referenced by tacReadScanditronics().

◆ roinameAddField()

char * roinameAddField	(	char *	roiname,
		const char *	field,
		const unsigned int	in,
		const unsigned int	count )

extern

Add a new field into given position inside an existing ROI name. Space character and '_' are accepted as field separators.

Returns: Returns pointer to the edited ROI name, or NULL in case of an error.

See also: roinameSubpart, roinameEditByTemplate, roinameExists, roinameMatch, strdelstr

Parameters

roiname	ROI name string which will be edited.
field	String which will be added as a new field in roiname; can be empty in order to add empty field, but NULL is not accepted.
in	Field index [0..n] for which position new field will be added; enter a large number to add field to the end.
count	Max length of the roiname, including space for terminal null. Function returns NULL (error) in case there is not enough space.

Definition at line 106 of file roiname.c.

  {
  if(roiname==NULL) return(roiname);
  if(field==NULL) return((char*)NULL);
  unsigned int slen=strlen(roiname);
  unsigned int flen=strlen(field);
  if(count < (slen+flen+2)) return((char*)NULL);
  if(slen==0) {
    if(flen==0) {strlcpy(roiname, "_", count); return(roiname);}
    else {strlcpy(roiname, field, count); return(roiname);}
  }
 
  char *c, subpart[slen+1], tmp[slen+1];
  unsigned int i=0;
  strcpy(tmp, roiname); strcpy(roiname, "");
  c=roinameSubpart(tmp, " _", i, subpart, slen+1);
  while(c) {
    if(i==in) {
      if(i>0) strlcat(roiname, "_", count);
      strlcat(roiname, field, count);
    }  
    if(i>0 || i==in) strlcat(roiname, "_", count);
    strlcat(roiname, subpart, count);
    c=roinameSubpart(tmp, " _", ++i, subpart, slen+1);
  }
  if(i<=in) {
    if(i>0) strlcat(roiname, "_", count);
    strlcat(roiname, field, count);
  }
  return(roiname);
}  

◆ roinameEditByTemplate()

char * roinameEditByTemplate	(	const char *	template,
		const char *	currname,
		char *	newname,
		const unsigned int	count )

extern

Make a new ROI name based on template and old ROI name.

Returns: Returns pointer to the new ROI name, or NULL in case of an error.

See also: roinameSubpart, roinameAddField, roinameExists, roinameMatch

Parameters

template	Template ROI name, where '_' is used to separate sub-parts, and '@' is replaced by the corresponding sub-part from the 'currname'.
currname	Current ROI name, from which sub-parts will be added as replacement of '@' in the template. Space character and '_' are assumed to separate the sub-parts.
newname	Pointer to pre-allocated string of length count, into which the new ROI name will be written.
count	Max length of new name, including terminal null.

Definition at line 65 of file roiname.c.

  {
  if(newname==NULL || count<1) return((char*)NULL);
  newname[0]=(char)0;
  if(template==NULL || currname==NULL) return((char*)NULL);
  
  unsigned int i, len;
  len=strlen(template); i=strlen(currname); if(i>len) len=i;
  if(len<1) return(newname);
  char subpart[++len], *c;
  
  i=0;
  c=roinameSubpart(template, "_", i, subpart, len);
  while(c) {
    /* If not '@' or if '@' can be replaced with sub-part from original tacname, then copy */
    if(strcmp(subpart, "@")!=0 || roinameSubpart(currname, " _", i, subpart, 256)) {
      if(i>0) strlcat(newname, "_", count);
      strlcat(newname, subpart, count);
    }
    c=roinameSubpart(template, "_", ++i, subpart, 256);
  }
  return(newname);
}

◆ 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.

See also: roinameSubpart, roinameEditByTemplate, roinameAddField, roinameMatch

Parameters

roiname ROI name string; not edited.

Definition at line 158 of file roiname.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);
}

Referenced by tacEnsureNames(), and tacIndividualNames().

◆ roinameMatch()

int roinameMatch	(	const char *	roiname,
		const char *	test_str,
		TPCSTATUS *	status )

extern

Test whether ROI (TAC) name matches with a test string.

Test string can contain wild cards. If both names contain ROI name separators (' ', '_', or '-'), the corresponding sub-names are tested against each other. If test_str does not contain separators but ROI name does, then names are considered matching if test_str matches with any of sub-names.

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

Author: Vesa Oikonen

See also: roinameSubpart, roinameEditByTemplate, roinameAddField, roinameExists, fnmatch, strstrNoQuotation

Parameters

roiname	Pointer to ROI name to be tested.
test_str	Test string.
status	Pointer to status data; enter NULL if not needed.

Definition at line 183 of file roiname.c.

  {
  int verbose=0; if(status!=NULL) verbose=status->verbose;
  if(verbose>0) printf("%s()\n", __func__);
  statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_OK);
  /* If both string are empty, then that is a match */
  if((roiname==NULL || strlen(roiname)<1) && (test_str==NULL || strlen(test_str)<1)) {
    return(1);
  } 
  /* If just either of those is empty, then that is a no match */
  if(roiname==NULL || strlen(roiname)<1) {
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_NO_DATA);
    return(0);
  }
  if(test_str==NULL || strlen(test_str)<1) {
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_NO_DATA);
    return(0);
  }
 
  /* Check the number of name sub-parts */
  int i, len;
  len=strlen(roiname); i=strlen(test_str); if(i>len) len=i;
  char subn1[len+1], subn2[len+1];
  int roisubnr=0, refsubnr=0;
  for(i=0; i<6; i++) {
    if(roinameSubpart(roiname, " _-", i, subn1, len+1)==NULL) break; 
    roisubnr++;
  }
  for(i=0; i<6; i++) {
    if(roinameSubpart(test_str, " _-", i, subn2, len+1)==NULL) break; 
    refsubnr++;
  }
  if(verbose>3) printf("roisubnr := %d\nrefsubnr := %d\n", roisubnr, refsubnr);
  if(roisubnr==0 && refsubnr==0) return(1);
  if(roisubnr==0 || refsubnr==0) return(0);
  if(refsubnr>1 && roisubnr!=refsubnr) return(0);
  
  /* If more than one sub-name to test for, then test against corresponding sub-names with 
     wild cards */
  if(refsubnr>1) {
    for(i=0; i<refsubnr; i++) {
      roinameSubpart(roiname, " _-", i, subn1, len+1); 
      if(!strcmp(subn1, ".")) strcpy(subn1, "");
      roinameSubpart(test_str, " _-", i, subn2, len+1);
      if(!strcmp(subn2, ".")) strcpy(subn2, "");
      if(verbose>3) printf(" '%s' vs '%s'\n", subn1, subn2);
      if(strnlen(subn1, 1)==0 && strnlen(subn2, 1)==0) continue;
      if(strnlen(subn2, 1)<1 || fncasematch(subn1, subn2)==0 ) return(0);
    }
    return(1);
  }
 
  /* So, we have just one sub-name in the test string */
  /* If it is found in any of sub names, accepting wild cards, then we consider it a match */
  if(roinameSubpart(test_str, " _-", 0, subn2, len+1)==NULL) return(0);
  if(!strcmp(subn2, ".")) strcpy(subn2, "");
  for(i=0; i<roisubnr; i++) {
    if(roinameSubpart(roiname, " _-", i, subn1, len+1)==NULL) continue;
    if(!strcmp(subn1, ".")) strcpy(subn1, "");
    if(fncasematch(subn1, subn2)) return(1);
  }
  return(0);
}

Referenced by parSelectTACs(), and tacSelectTACs().

◆ roinameSubpart()

char * roinameSubpart	(	const char *	roiname,
		const char *	dlm,
		const unsigned int	si,
		char *	subpart,
		const unsigned int	slen )

extern

Divide ROI name into its sub-parts, based on user-defined delimiter characters (usually underscore and/or space).

Author: Vesa Oikonen

Returns: Returns pointer to the required sub-part, or NULL if not found.

See also: roinameEditByTemplate, roinameAddField, roinameExists, roinameMatch

Parameters

roiname	Pointer to ROI name. Not modified.
dlm	String containing the sub-part delimiters. For example, " _".
si	Index [0..n-1] of sub-part to search for.
subpart	Pointer to allocated string for sub-part.
slen	Length of the sub-part string, including terminating null. Resulting sub-part length will be at maximum slen-1.

Definition at line 20 of file roiname.c.

  {
  if(subpart==NULL || slen<1) return(subpart);
  subpart[0]=(char)0;
  if(roiname==NULL || strnlen(roiname, 1)==0) return((char*)NULL);
  if(dlm==NULL || strnlen(dlm, 1)==0) return((char*)NULL);
 
  unsigned int i=0;
  char *s=strdup(roiname);
  char *sp=s; //printf(" si=%d\n", si);
  do {
    //printf("  sp='%s'\n", sp);
    char *p=strpbrk(sp, dlm);
    //printf("  i=%d  p='%s'\n", i, p);
    if(!p) {
      if(si==i) {strlcat(subpart, sp, slen); free(s); return(subpart);}
      free(s); return((char*)NULL);
    }
    *p=(char)0;
    if(si==i) {strlcat(subpart, sp, slen); free(s); return(subpart);}
    sp=p+1; i++;
  } while(*sp);  
 
  if(si==i) {strlcat(subpart, sp, slen); free(s); return(subpart);}
  free(s);
  return((char*)NULL);
}

Referenced by parReadFIT(), parReadRES(), parWriteFIT(), parWriteRES(), roinameAddField(), roinameEditByTemplate(), and roinameMatch().

◆ statusFree()

void statusFree ( TPCSTATUS * s )

extern

After last use, free any memory allocated in the structure.

See also: errorMsg, statusSet, statusInit

Parameters

s	Pointer to status.

Definition at line 126 of file statusmsg.c.

  {
  if(s==NULL) return;
  free(s->last_file); free(s->last_func);
  free(s->current_file); free(s->current_func);
  if(s->fp!=NULL && s->fp!=stdout && s->fp!=stderr) fclose(s->fp);
  statusInit(s);
}

◆ statusInit()

void statusInit ( TPCSTATUS * s )

extern

Before first use, initialize the structure contents using this function. After last use of status structure, call statusFree(), which also calls this function to re-initiate the structure contents.

See also: errorMsg, statusSet, statusFree

Parameters

s	Pointer to status.

Definition at line 104 of file statusmsg.c.

  {
  if(s==NULL) return;
  s->last_line=0;
  s->last_file=NULL;
  s->last_func=NULL;
  s->current_line=0;
  s->current_file=NULL;
  s->current_func=NULL;
  s->verbose=0;
  s->fp=stdout;
  s->forgiving=0;
  s->error=0;
}

Referenced by statusFree().

◆ statusPrint()

void statusPrint	(	FILE *	fp,
		TPCSTATUS *	s )

extern

Print the contents of status structure.

Author: Vesa Oikonen

See also: statusSet, errorMsg, statusInit, statusFree

Parameters

fp	File pointer for the output.
s	Pointer to status.

Definition at line 82 of file statusmsg.c.

  {
  if(fp==NULL || s==NULL) return;
  if(s->current_file!=NULL && s->current_func!=0 && s->current_line>0)
    fprintf(fp, "In %s() on line %d of file %s\n", s->current_func, s->current_line, s->current_file);
  if(s->verbose<=0) return;
  if(s->last_file!=NULL && s->last_func!=0 && s->last_line>0)
    fprintf(fp, "called in %s() on line %d of file %s\n", s->last_func, s->last_line, s->last_file);
  return;
}

◆ statusSet()

void statusSet	(	TPCSTATUS *	s,
		const char *	func,
		const char *	srcfile,
		int	srcline,
		tpcerror	error )

extern

Set the status contents.

See also: statusInit, statusFree, errorMsg, tpcProgramName, tpcProcessStdOptions

Parameters

s	Pointer to status.
func	Pointer to function name, usually func.
srcfile	Pointer to source file name, usually FILE.
srcline	Pointer to source file line number, usually LINE.
error	Error code, 0 if status is ok.

Definition at line 142 of file statusmsg.c.

  {
  if(s==NULL) return;
  s->error=error;
  /* Delete any previous content */
  s->last_line=0;
  free(s->last_file); free(s->last_func);
  s->last_file=NULL; s->last_func=NULL;
  /* Move current content to previous */
  s->last_line=s->current_line;
  s->last_file=s->current_file;
  s->last_func=s->current_func;
  /* Set current contents */
  s->current_line=srcline;
  if(func!=NULL) s->current_func=strdup(func);
  if(srcfile!=NULL) s->current_file=strdup(srcfile);
  s->error=error;
}

Referenced by abssWrite(), anaExists(), bfm1TCM(), bfmSRTM(), csvPutLine(), csvPutLineWithSpaces(), csvRead(), csvReorg(), csvWrite(), dcmFileList(), dcmFileRead(), dcmFileRemove(), dcmFileWrite(), dcmMListRead(), dcmmlSortByPlane(), ecatMListRead(), ecatReadMainheader(), ecatWriteMainheader(), iftCopyItems(), iftDeleteDuplicateKeys(), iftPut(), iftPutFromString(), iftPutFromStringWithSpaceSeparator(), iftRead(), iftReplaceKey(), iftReplaceValue(), iftWrite(), iftWriteItem(), imgAllocate(), imgCompareConc(), imgCompareTimes(), imgCompareUnit(), imgFillOHeader(), imgRead(), imgReadDICOM(), imgReadNifti(), imgWrite(), imgWriteDICOM(), imgWriteNifti(), micropetExists(), micropetHeaderRead(), mtgaPlotSVG(), niftiExists(), nlopt1D(), nloptIATGO(), nloptITGO1(), nloptITGO2(), nloptMPSO(), nloptPowellBrent(), nloptSimplex(), nloptSimplexARRS(), parAllocateWithTAC(), parCombineTACs(), parCompareParameterNames(), parCompareParameters(), parCompareParameterUnits(), parCompareTacNames(), parCompareWSS(), parExamplePerfectBolus(), parExampleTTACs(), parFromIFT(), parRead(), parReadCSV(), parReadFIT(), parReadRES(), parSelectByAnother(), parSelectParameters(), parSelectTACs(), parSortByName(), parToIFT(), parWrite(), parWriteCSV(), parWriteFIT(), parWriteIFT(), parWriteRES(), parWriteXML(), roinameMatch(), sifWeight(), spectralDExp(), spectralKRange(), tacAddZeroSample(), tacAllocateWithPAR(), tacAUC(), tacCompareConc(), tacCompareNames(), tacCompareTimes(), tacCompareUnit(), tacCompareWeights(), tacCorrectFrameOverlap(), tacDecayCorrection(), tacDelay(), tacFittime(), tacFormatDetermine(), tacFramesToSteps(), tacGetHeaderInjectiontime(), tacGetHeaderIsotope(), tacGetHeaderScanstarttime(), tacGetHeaderStudynr(), tacGetHeaderTimeunit(), tacGetHeaderUnit(), tacInput2sim(), tacIntegrate(), tacInterpolate(), tacInterpolateInto(), tacInterpolateToEqualLengthFrames(), tacPlotFitSVG(), tacPlotHistogramSVG(), tacPlotLineSVG(), tacRead(), tacRead4DM(), tacReadAllogg(), tacReadCarimasTxt(), tacReadCSV(), tacReadDFT(), tacReadGEMS(), tacReadHRPLUSHC(), tacReadHRRTHC(), tacReadInveonCSV(), tacReadMat(), tacReadModelingData(), tacReadModelingInput(), tacReadOldAllogg(), tacReadPMOD(), tacReadQView(), tacReadReference(), tacReadScanditronics(), tacReadSIF(), tacReadSimple(), tacSelectTACs(), tacSetWeights(), tacSetX(), tacSortByAUC(), tacSortByConc(), tacSortByName(), tacSortByTime(), tacToPAR(), tacVb(), tacVerifyTimeOrder(), tacWByFreq(), tacWeightModerate(), tacWeightNorm(), tacWMove(), tacWrite(), tacWriteCSV(), tacWriteDFT(), tacWritePMOD(), tacWriteSheetIntoXML(), tacWriteSIF(), tacWriteSimple(), tacWriteXML(), tacXUnitConvert(), tacYUnitConvert(), tacYUnitMass2Volume(), and tacYUnitVolume2Mass().

◆ 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.

See also: strdelstr, strstrNoQuotation

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 155 of file stringext.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 csvSearchField(), dcmImgIsotope(), dcmMListRead(), elementIdentify(), iftSearchKey(), iftSearchValue(), imgFormatFromFName(), parFormatFromExtension(), parReadRES(), parWriteRES(), strdelstr(), tacReadInveonCSV(), tacSelectBestReference(), tacWritePMOD(), tacWriteXML(), and unitIdentifyFilename().

◆ strChrCount()

int strChrCount	(	const char *	s1,
		const char *	s2 )

extern

Count how many times characters specified in string s2 are found in string s1.

Search is case-sensitive.

Author: Vesa Oikonen

Returns: Returns the nr of characters of s2 found in s1.

See also: strReplaceChar, strstrNoQuotation, strUppercaseCount

Parameters

s1	String to search for characters; not modified.
s2	String containing characters which are searched for; not modified.

Definition at line 98 of file stringext.c.

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

Referenced by strncpyClean().

◆ strClean()

int strClean ( char * s )

extern

Removes any initial and trailing space characters from specified string s. Removes also quotation marks ("" or '') that may be around the string s. Quotes or space characters in the middle of the string are not removed.

Returns: 0 when successful, otherwise >0.

Author: Vesa Oikonen

See also: strncpyCleanSpaces, strCleanSpaces, strIsSpaceOnly, strCleanPars

Parameters

s	Pointer to the string.

Definition at line 389 of file stringext.c.

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

Referenced by iftPutFromString(), iftPutFromStringWithSpaceSeparator(), imgReadDICOM(), tacNameSplit(), and tacRead4DM().

◆ 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 769 of file stringext.c.

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

◆ strCleanPars()

void strCleanPars ( char * s )

extern

Remove parenthesis (s), [s], or {s} from around the string s.

Only one set of pars are removed, and only if the same pars are as the first and last character of the string.

Precondition: If necessary, remove possible spaces or quotation marks before this.

See also: strInPars, strCleanSpaces

Author: Vesa Oikonen

Parameters

s	Pointer to the string.

Definition at line 476 of file stringext.c.

  {
  if(s==NULL || strInPars(s)==0) return;
  int len=strlen(s);
  s[--len]=(char)0;
  for(int i=1; i<=len; i++) s[i-1]=s[i];
  return;
}

Referenced by iftGetDoubleWithUnit(), tacRead4DM(), tacReadCarimasTxt(), and tacReadDFT().

◆ 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

See also: strClean, strncpyClean, strncpyCleanSpaces, strIsSpaceOnly

Parameters

s	Pointer to the string.

Definition at line 300 of file stringext.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 csvCleanSpaces(), dcmImgIsotope(), and strncpyClean().

◆ strDateRead()

int strDateRead	(	const char *	str,
		struct tm *	date )

extern

Reads date from a standard string representation of date.

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

See also: strDateTimeRead, strDateTimeValid, ctime_r_int

Author: Vesa Oikonen

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.
date	Pointer to allocated struct tm where the date is written; undefined if date is not valid.

Definition at line 378 of file datetime.c.

  {
  char buf[32];
  int ret, n, YYYY=0, MM=0, DD=0;
 
  if(strnlen(str, 10)<8) return 1;
  if((ret=strDateValid(str))<=0) {
    strlcpy(buf, str, 11); 
  } else {
    ret=strDateValid2(str, buf);
    if(ret>0) ret=strDateValid3(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 tacReadGEMS(), and tacReadScanditronics().

◆ strDateTimeAdd()

int strDateTimeAdd	(	int	s,
		char *	dt )

extern

Add given time in seconds to the date and time.

Returns: 0 if successful.

See also: tmAdd, tmDifference, time_to_tm, strDateTimeDifference

Parameters

s	Time to add in seconds; can be negative.
dt	Pointer to date and time string, edited here. It must be allocated for at least 20 characters.

Definition at line 537 of file datetime.c.

  {
  if(dt==NULL) return(1);
  struct tm tm0;
  if(strDateTimeRead(dt, &tm0)) return(1);
  if(s==0) return(0);
  tmAdd(s, &tm0);
  if(!strftime(dt, 20, "%Y-%m-%d %H:%M:%S", &tm0)) return(2);
  return(0);
}

◆ strDateTimeDifference()

double strDateTimeDifference	(	const char *	dt1,
		const char *	dt0 )

extern

Calculate the difference in seconds between two dates and times in international date and time string format (YYYY-MM-DD hh:mm:ss).

See also: tmDifference, strDateTimeAdd, strDateTimeRead, strDateTimeValid, ctime_r_int

Returns: Returns dt1-dt0 time difference in seconds.

Parameters

dt1	Pointer to date and time string.
dt0	Pointer to date and time string.

Definition at line 519 of file datetime.c.

  {
  if(dt1==NULL || dt0==NULL) return(0.0);
  struct tm tm1, tm0;
  if(strDateTimeRead(dt1, &tm1) || strDateTimeRead(dt0, &tm0)) return(0.0);
  return(tmDifference(&tm1, &tm0));
}

Referenced by dcmMListRead(), and imgReadDICOM().

◆ strDateTimeRead()

int strDateTimeRead	(	const char *	str,
		struct tm *	date )

extern

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

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

See also: strDateRead, strDateTimeValid, ctime_r_int, tmDifference, strClean

Author: Vesa Oikonen

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
date	Pointer to allocated struct tm where date and time is written. undefined if date or time is not valid.

Definition at line 350 of file datetime.c.

  {
  char buf[32];
  int ret, n, YYYY=0, MM=0, DD=0, hh=0, mm=0, ss=0;
 
  ret=strDateTimeValid(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;
  ret=strftime(buf, 32, "%Y-%m-%d %H:%M:%S", date);
  if(ret==0) return(50);
  return(0);
}

Referenced by abssWrite(), strDateTimeAdd(), strDateTimeDifference(), tacReadAllogg(), tacReadGEMS(), and tacReadScanditronics().

◆ strDateTimeValid()

int strDateTimeValid	(	const char *	str,
		char *	intdate )

extern

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

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

Author: Vesa Oikonen

See also: strDateTimeRead, strTimeValid, strDateValid, strClean

Parameters

str	String to be verified.
intdate	Pointer to preallocated string where date and time is written in international format (YYYY-MM-DD hh:mm:ss); enter NULL, if not needed. Undefined, if date is not valid.

Definition at line 308 of file datetime.c.

  {
  int ret1, ret2, len;
  char correct_date[20], *time_ptr=NULL, tmp[20];
 
  if(str==NULL) return 1;
  // 1/1/70 00:00:00 or 1970-01-01 00:00:00
  if(strnlen(str, 19)<15) return 1;
  /* Separate date and time */
  len=strcspn(str, " \t");
  if(len<6 || len>10) return 1;
  strlcpy(tmp, str, len+1); time_ptr=(char*)str+len+1;
  /* and check the date */
  correct_date[0]=(char)0;
  if((ret1=strDateValid(tmp))<=0) {strlcpy(correct_date, tmp, 20);}
  else if((ret1=strDateValid2(tmp, correct_date))<=0) {}
  else if((ret1=strDateValid3(tmp, correct_date))<=0) {}
  else {return 2;}
 
  /* the check the time */
  ret2=strTimeValid(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(), strDateTimeRead(), tacReadOldAllogg(), and tacReadSIF().

◆ strDateValid()

int strDateValid ( const char * str )

extern

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

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

Author: Vesa Oikonen

See also: ctime_r_int, strDateValid2, strDateValid3, strDateValid4, strDateTimeValid, strClean

Parameters

str	String to be verified.

Definition at line 155 of file datetime.c.

  {
  if(str==NULL || strnlen(str, 10)<10) return 1;
  if(str[4]!='-' || str[7]!='-') return 2;
  if(strncasecmp(str, "YYYY-MM-DD", 10)==0) return -1;
  int Y, M, D, n;
  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(), strDateRead(), and strDateTimeValid().

◆ strDateValid2()

int strDateValid2	(	const char *	str,
		char *	intdate )

extern

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

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

Author: Vesa Oikonen

See also: ctime_r_int, strDateValid, strDateValid3, strDateValid4, strDateTimeValid

Parameters

str	String to be verified.
intdate	Pointer to preallocated string where date is written in international format; enter NULL, if not needed. Undefined, if date is not valid.

Definition at line 179 of file datetime.c.

  {
  if(str==NULL || strnlen(str, 10)<8) return 1;
  int n; char *cptr;
  cptr=(char*)str; n=strcspn(cptr, "./"); if(n<1 || n>2) return 2;
  cptr+=n+1; n=strcspn(cptr, "./"); if(n<1 || n>2) return 2;
  cptr+=n+1; n=strcspn(cptr, "./- \t\n\r"); if(n!=4) return 2;
  if(strncasecmp(str, "DD.MM.YYYY", 10)==0 || strncasecmp(str, "DD/MM/YYYY", 10)==0) {
    return -1;
  }
  int Y, M, D;
  n=sscanf(str, "%d/%d/%4d", &D, &M, &Y);
  if(n<3) n=sscanf(str, "%d.%d.%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 strDateRead(), and strDateTimeValid().

◆ strDateValid3()

int strDateValid3	(	const char *	str,
		char *	intdate )

extern

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

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

Author: Vesa Oikonen

See also: ctime_r_int, strDateValid2, strDateValid, strDateValid4, strDateTimeValid

Parameters

str	String to be verified.
intdate	Pointer to allocated string where date is written in international format; enter NULL, if not needed. Undefined, if date is not valid.

Definition at line 214 of file datetime.c.

  {
  if(str==NULL || strnlen(str, 8)<6) return 1;
  int n; char *cptr;
  cptr=(char*)str; n=strcspn(cptr, "./"); if(n<1 || n>2) return 2;
  cptr+=n+1; n=strcspn(cptr, "./"); if(n<1 || n>2) return 2;
  cptr+=n+1; n=strcspn(cptr, "./- \t\n\r"); if(n!=2) return 2;
  if(strncasecmp(str, "DD.MM.YY", 8)==0 || strncasecmp(str, "DD/MM/YY", 8)==0) {
    return -1;
  }
  int Y, M, D;
  n=sscanf(str, "%d/%d/%2d", &D, &M, &Y);
  if(n<3) n=sscanf(str, "%d.%d.%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 strDateRead(), and strDateTimeValid().

◆ strDateValid4()

int strDateValid4	(	int	dateint,
		char *	intdate,
		int *	year,
		int *	month,
		int *	day )

extern

Verify that specified integer contains date in format YYYYMMDD.

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

Author: Vesa Oikonen

See also: ctime_r_int, strDateValid2, strDateValid3, strDateValid, strDateTimeValid

Parameters

dateint	Integer 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. Undefined, if date is not valid.
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 250 of file datetime.c.

  {
  if(dateint<18000101 || dateint>99991231) return 1;
  int Y, M, D, n;
  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;
  if(intdate!=NULL) {sprintf(intdate, "%04d-%02d-%02d", Y, M, D);}
  return 0;
}

◆ strdelstr()

char * strdelstr	(	char *	s1,
		const char *	s2 )

extern

Find (case-insensitive) s2 in s1, and delete it from s1.

Returns: Returns pointer to s1 at where the deletion stopped, or to the start of s1 in case s2 was not found.

See also: strReplaceChar, strClean, strTrimleft

Parameters

s1	Pointer to string in which sub-string s2 is searched; modified.
s2	Pointer to sub-string which is searched for in source string s1.

Definition at line 688 of file stringext.c.

  {
  if(s1==NULL || s2==NULL) return(s1); //printf("strdelstr('%s', '%s')\n", s1, s2);
  int n=strlen(s2); if(n<1) return(s1);
  char *cptr, *rptr;
  cptr=strcasestr(s1, s2);
  if(cptr==NULL || *cptr==(char)0) return(s1);
  rptr=cptr;
  while(*cptr) {*cptr=cptr[n]; cptr++;}
  return(rptr);
}

◆ strdup()

char * strdup ( const char * s )

extern

Allocates memory and copies into it the string addressed by s, including the terminating character.

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

Returns: pointer to allocated string, or NULL in case of error.

Parameters

s	String to be duplicated.

Definition at line 185 of file stringext.c.

  {
  if(s==NULL) return NULL;
  size_t length=strlen(s)+1;
  void *r;
  r=calloc(length, sizeof(char)); if(r==NULL) return NULL;
  return (char*)memcpy(r, s, length);
}

Referenced by atofCheck(), atofVerified(), csvCellReplace(), csvDuplicate(), csvPutString(), dcmFileList(), dcmMListRead(), iftPut(), iftPutFromString(), iftPutFromStringWithSpaceSeparator(), iftReplaceKey(), iftReplaceValue(), parCombineTACs(), parCompareParameterNames(), parCompareTacNames(), parSelectByAnother(), pathCreate(), pathExist(), roinameSubpart(), statusSet(), strCleanSpaces(), strTokenDup(), studynrFromFilename(), tacNameSplit(), tacRead(), tacReadCarimasTxt(), tacReadCSV(), tacReadDFT(), tacReadHRRTHC(), tacReadPMOD(), tpcHtmlUsage(), tpcPrintBuild(), tpcPrintUsage(), and tpcProgramName().

◆ 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: strCleanForXML

Parameters

s	Pointer to the string to be encoded.

Definition at line 731 of file stringext.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 parWriteXML(), and tacWriteXML().

◆ strHaveDecimalComma()

int strHaveDecimalComma ( const char * s )

extern

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

Author: Vesa Oikonen

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

See also: strHaveDecimalSeparator, atofVerified, strIsValidNumber

Parameters

s	Pointer to string.

Definition at line 115 of file decpoint.c.

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

Referenced by tacReadCSV().

◆ strHaveDecimalSeparator()

int strHaveDecimalSeparator ( const char * s )

extern

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

Author: Vesa Oikonen

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

See also: strHaveDecimalComma, atofVerified, strIsValidNumber

Parameters

s	Pointer to string.

Definition at line 132 of file decpoint.c.

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

◆ strInPars()

int strInPars ( char * s )

extern

Checks whether string is inside parenthesis (s), [s], or {s}, that is, the same pars as the first and last character of the string.

Precondition: If necessary, remove possible spaces or quotation marks before this.

Returns: 0 if not in pars, 1 if inside pars.

Author: Vesa Oikonen

See also: strCleanPars, strCleanSpaces, strIsValidNumber

Parameters

s	Pointer to the string.

Definition at line 453 of file stringext.c.

  {
  if(s==NULL) return 0;
  int len=strlen(s); if(len<2) return 0;
  if(s[0]=='(' && s[len-1]==')') return 1;
  if(s[0]=='{' && s[len-1]=='}') return 1;
  if(s[0]=='[' && s[len-1]==']') return 1;
  return 0;
}

Referenced by strCleanPars(), and tacReadDFT().

◆ strIsSpaceOnly()

int strIsSpaceOnly ( char * s )

extern

Checks whether string contains only space characters.

Text files saved in text editor or Excel may accidentally contain empty looking lines which still may contain tabs and spaces, leading to problems when reading them for example as CSV files.

Returns: 0 if lines contains something more than spaces, 1 if only space characters or string is empty.

Author: Vesa Oikonen

See also: strClean, strncpyCleanSpaces

Parameters

s	Pointer to the string.

Definition at line 671 of file stringext.c.

  {
  if(s==NULL) return(1);
  int len=strlen(s); if(len<1) return(1);
  for(int i=0; i<len; i++) if(!isspace(s[i])) return(0);
  return 1;
}

Referenced by csvRead().

◆ strIsValidNumber()

int strIsValidNumber ( const char * s )

extern

Verifies that given string seems like a valid representation of integer or floating point number in decimal or exponential format.

Returns: 1 if string is valid number, and 0 if not.

Author: Vesa Oikonen

See also: atofVerified, atoiCheck, tpcYesNo

Parameters

s	Pointer to the string to be verified; it must not contain any leading or trailing space characters etc, or it will never verify.

Definition at line 33 of file decpoint.c.

  {
  if(s==NULL || strnlen(s, 256)<1) return 0;
 
  char *p=(char*)s;
  int i, nn=0;
  /* It can start with + or minus */
  if(*p=='+' || *p=='-') p++;
  /* Then jump over all digits, counting them */
  i=strspn(p, "0123456789"); nn+=i; p+=i;
  if(!*p) {if(nn>0) return 1; else return 0;} // end of string
  /* Then there can be a decimal point or comma, but only one per string */
  if(*p=='.' || *p==',') p++;
  if(!*p) {if(nn>0) return 1; else return 0;} // end of string
  /* Again jump over all digits, counting them */
  i=strspn(p, "0123456789"); nn+=i; p+=i;
  /* At this point we must have got at least one digit */
  if(nn<1) return 0;
  if(!*p) return 1; // end of string
  /* String continues, the next character must be E or e */
  if(*p=='E' || *p=='e') p++; else return 0;
  /* next can be + or - but not necessarily */
  if(*p=='+' || *p=='-') p++;
  /* then we must have at least one digit */
  i=strspn(p, "0123456789"); if(i<1) return 0;
  nn+=i; p+=i;
  /* now we must be at the end of string */
  if(!*p) return 1; else return 0; 
}

Referenced by atofCheck(), and atofVerified().

◆ 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 592 of file stringext.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 fileBackup(), micropetExists(), parReadFIT(), parReadRES(), roinameAddField(), roinameEditByTemplate(), roinameSubpart(), and tpcProgramName().

◆ 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 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 632 of file stringext.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 abssWrite(), anaExists(), dcmAddItem(), dcmFileRead(), dcmMListRead(), dcmSameImage(), ecatReadMainheader(), ecatWriteMainheader(), fileBackup(), imgCopyHeader(), imgSetNiftiHeader(), imgWriteNifti(), micropetExists(), niftiCreateFNames(), niftiExists(), parReadCSV(), parReadFIT(), parReadRES(), parWriteFIT(), roinameAddField(), strDateRead(), strDateTimeValid(), strncpyClean(), strncpyCleanSpaces(), strTokenDup(), strTokenNCpy(), studynrFromFilename(), tacGetHeaderIsotope(), tacGetHeaderStudynr(), tacNameSplit(), tacRead4DM(), tacReadCarimasTxt(), tacReadCSV(), tacReadDFT(), tacReadInveonCSV(), tacReadMat(), tacReadOldAllogg(), tacReadPMOD(), tacReadQView(), tacReadSIF(), tacReadSimple(), and tpcProgramName().

◆ strncatDouble()

char * strncatDouble	(	char *	str1,
		const double	d,
		size_t	count )

extern

Concatenate not more than count characters of the string representation of double d to the string pointed to by str1. String is terminated with a null.

See also: strncatInt

Returns: Pointer to str1.

Parameters

str1	String to be catenated.
d	Double value to be catenated as string into str1.
count	No more than count characters are concatenated to str1. If str1 was allocated like char str1[MAX_LEN], then set count=MAX_LEN-strlen(str1)-1.

Definition at line 517 of file stringext.c.

  {
  char tmp[128]; sprintf(tmp, "%g", d);
  return strncat(str1, tmp, count);
}

◆ strncatInt()

char * strncatInt	(	char *	str1,
		const int	n,
		size_t	count )

extern

Concatenate not more than count characters of the string representation of integer n to the string pointed to by str1. String is terminated with a null.

See also: strncatDouble, strncatIntZ, strIsValidNumber

Returns: Pointer to str1.

Parameters

str1	String to be catenated.
n	Integer value to be catenated as string into str1.
count	No more than count characters are concatenated to str1. If str1 was allocated like char str1[MAX_LEN], then set count=MAX_LEN-strlen(str1)-1.

Definition at line 495 of file stringext.c.

  {
  char tmp[128]; sprintf(tmp, "%d", n);
  return strncat(str1, tmp, count);
}

◆ strncatIntZ()

char * strncatIntZ	(	char *	str1,
		const int	n,
		const int	maxn,
		size_t	count )

extern

Concatenate not more than count characters of the zero-padded string representation of integer n to the string pointed to by str1. String is terminated with a null.

See also: strncatInt, strncatDouble

Returns: Pointer to str1.

Parameters

str1	String to be catenated.
n	Integer value to be catenated as string into str1.
maxn	Integer value which is large enough that it would not need to be zero-padded.
count	No more than count characters are concatenated to str1. If str1 was allocated like char str1[MAX_LEN], then set count=MAX_LEN-strlen(str1)-1.

Definition at line 539 of file stringext.c.

  {
  if(n<0) strncat(str1, "-", count);
  char tmp[128]; sprintf(tmp, "%d", abs(maxn));
  int len=strlen(tmp);
  sprintf(tmp, "%0*d", len, abs(n));
  return strncat(str1, tmp, count);
}

◆ strncpyClean()

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

extern

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

Returns: the length of the new string s1.

Author: Vesa Oikonen

See also: strncpyCleanSpaces, strCleanSpaces, strClean, strIsSpaceOnly, strstrNoQuotation

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 321 of file stringext.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 *tmp;
  int i, sqn, dqn, m;
 
  /* Allocate temp string */
  tmp=calloc(strlen(s2)+1, sizeof(char));
 
  /* Trim space characters */
  i=strncpyCleanSpaces(tmp, s2, strlen(s2)+1);
  if(i<1) {strcpy(s1, ""); free(tmp); return(0);}
 
  /* Loop as long as there is something to trim */
  i=1;
  while(i) {
    m=0; // nothing modified so far  
    /* Count quotes */
    sqn=strChrCount(tmp, "\'");
    dqn=strChrCount(tmp, "\"");
    if(sqn==0 && dqn==0) {i=0; break;}
    /* Trim matching quotation marks */
    if(dqn>=2 && tmp[0]=='\"' && tmp[strlen(tmp)-1]=='\"') {
      tmp[0]=' '; tmp[strlen(tmp)-1]=(char)0; dqn-=2; m+=2;
    }
    if(sqn>=2 && tmp[0]=='\'' && tmp[strlen(tmp)-1]=='\'') {
      tmp[0]=' '; tmp[strlen(tmp)-1]=(char)0; sqn-=2; m+=2;
    }
    /* Trim orphan quotation marks */
    if(sqn%2) {
      if(tmp[0]=='\'') {tmp[0]=' '; sqn--; m++;}
      else if(tmp[strlen(tmp)-1]=='\'') {tmp[strlen(tmp)-1]=(char)0; sqn--; m++;}
    }
    if(dqn%2) {
      if(tmp[0]=='\"') {tmp[0]=' '; dqn--; m++;}
      else if(tmp[strlen(tmp)-1]=='\"') {tmp[strlen(tmp)-1]=(char)0; dqn--; m++;}
    }
    /* If nothing was modified, then do not continue */
    if(m==0) {i=0; break;}
    /* Trim spaces, quit if error encountered */
    if(strCleanSpaces(tmp)) i=0;
  }
 
  strlcpy(s1, tmp, maxlen);
  free(tmp);
  return(strlen(s1));
}

Referenced by imgReadDICOM(), parSelectParameters(), parSelectTACs(), strClean(), and tacSelectTACs().

◆ 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

See also: strCleanSpaces, strncpyClean, strIsSpaceOnly

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 265 of file stringext.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 niftiBasename(), parReadFIT(), strCleanSpaces(), strncpyClean(), and tacRead4DM().

◆ strndup()

char * strndup	(	const char *	s,
		size_t	n )

extern

Allocates memory and copies into it first n characters from string addressed by s, and adds the terminating character.

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

Returns: pointer to allocated string, or NULL in case of error.

Parameters

s	String to be duplicated.
n	Max nr of characters to copy, not including terminal zero.

Definition at line 205 of file stringext.c.

  {
  size_t length=strnlen(s, n);
  void *r; r=calloc(length+1, sizeof(char)); if(r==NULL) return NULL;
  return (char*)memcpy(r, s, length);
}

Referenced by csvPutLine(), csvPutLineWithSpaces(), iftPutFromString(), iftPutFromStringWithSpaceSeparator(), and tpcProgramName().

◆ 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, whicever 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 566 of file stringext.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 abssIdFromFName(), anaExists(), anaReadHeader(), atoiCheck(), dcmAddItem(), dcmDA2intl(), dcmDT2intl(), dcmFileList(), dcmFileRead(), dcmFileRemove(), dcmFileWrite(), dcmSOPIdentify(), dcmSOPUIDName(), dcmTM2intl(), dcmTrUID(), elementIdentify(), fileExist(), filenameGetExtension(), filenameGetExtensions(), filenameRmExtension(), filenameRmExtensions(), filenameRmFile(), filenameRmPath(), imgFormatFromFName(), imgRead(), imgReadDICOM(), imgReadNifti(), imgWrite(), imgWriteDICOM(), imgWriteNifti(), isotopeIdentify(), micropetExists(), micropetHeaderRead(), mtgaPlotSVG(), niftiExists(), niftiReadHeader(), niftiWriteHeader(), parCombineTACs(), parSelectByAnother(), parWriteFIT(), pathCreate(), pathFileList(), pathFileNr(), roinameMatch(), roinameSubpart(), strDateRead(), strDateTimeValid(), strDateValid(), strDateValid2(), strDateValid3(), strIsValidNumber(), strlcat(), strndup(), strTimeValid(), strTokenNr(), studynrFromFilename(), tacPlotFitSVG(), tacPlotHistogramSVG(), tacPlotLineSVG(), tacReadAllogg(), tacReadGEMS(), tacReadMat(), tacReadModelingData(), tacReadModelingInput(), tacReadOldAllogg(), tacReadQView(), tacReadReference(), tacReadScanditronics(), tpcHtmlUsage(), tpcProcessStdOptions(), tpcProgramName(), tpcYesNo(), unitIdentify(), and unitIdentifyFilename().

◆ strReplaceChar()

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

extern

Replace all characters c1 in string s with character c2.

See also: strChrCount, strdelstr, strTrimleft, strstrNoQuotation

Author: Vesa Oikonen

Parameters

s	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 134 of file stringext.c.

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

Referenced by csvPutDouble(), parCombineTACs(), parCompareParameterNames(), parCompareTacNames(), parSelectByAnother(), parToIFT(), strCleanForXML(), and tacWritePMOD().

◆ strstrNoQuotation()

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

extern

Version of strstr() which ignores any matches inside quotation marks "" or ''.

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

Author: Vesa Oikonen

Remarks: This function is an enhanced version of the one from libtpcmisc.

See also: strClean, strCleanSpaces

Parameters

haystack	Pointer to string to be searched.
needle	Pointer to string with quotation marks.

Definition at line 225 of file stringext.c.

  {
  if(haystack==NULL) return((char*)NULL);
  if(needle==NULL) return((char*)haystack);
 
  size_t i, test_len;
  int single_quotation=0;
  int double_quotation=0;
  char *cptr;
  
  test_len=strlen(needle); if(test_len<1) return((char*)haystack);
  for(i=0, cptr=(char*)haystack; i<strlen(haystack); i++, cptr++) {
    if(*cptr=='\'') {
      if(single_quotation==0 && strchr(cptr+1, '\'')!=NULL) single_quotation=1;
      else single_quotation=0;
      continue;
    }
    if(*cptr=='\"') {
      if(double_quotation==0 && strchr(cptr+1, '\"')!=NULL) double_quotation=1;
      else double_quotation=0;
      continue;
    }
    if(single_quotation==1 || double_quotation==1) continue;
    if(strncmp(cptr, needle, test_len)==0) return(cptr);
  }
  return((char*)NULL);
}

Referenced by csvRead(), iftPutFromString(), and iftPutFromStringWithSpaceSeparator().

◆ strTimeValid()

int strTimeValid ( const char * str )

extern

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

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

Author: Vesa Oikonen

See also: strDateTimeRead, strDateTimeValid, strDateValid

Parameters

str	String to be verified.

Definition at line 284 of file datetime.c.

  {
  if(str==NULL || strnlen(str, 8)<8) return 1;
  if(str[2]!=':' || str[5]!=':') return 2;
  if(strncasecmp(str, "hh:mm:ss", 8)==0) return -1;
  int h, m, s, n;
  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(), iftPutFromString(), and strDateTimeValid().

◆ 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

See also: strTokenNCpy, strstrNoQuotation

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 413 of file stringext.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 csvRead(), iftRead(), tacFormatDetermine(), and tacRead4DM().

◆ strTokenNCpy()

int strTokenNCpy	(	const char *	s1,
		const char *	s2,
		int	i,
		char *	s3,
		int	count )

extern

Copy the i'th token in the string s1 into string s3. The characters making up the string s2 are the delimiters that determine the tokens.

See also: strTokenNr, strTokenDup

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

Author: Vesa Oikonen

Parameters

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

Definition at line 53 of file stringext.c.

  {
  int j=0, n=0;
  char *cptr;
  if(s3!=NULL && count>0) s3[0]=(char)0;
  if(s1==NULL || s2==NULL || strlen(s1)==0 || strlen(s2)==0) return(0);
  if(i<1 || s3==NULL || count<2) return(0);
 
  cptr=(char*)s1;
  do {
    // pass delimiter characters
    j=strspn(cptr, s2); cptr+=j;
    // pass characters between delimiters
    j=strcspn(cptr, s2); 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(s3, cptr, j+1); //strncpy(s3, cptr, j); s3[j]=(char)0;
      break;
    }
    cptr+=j;
  } while(j>0);
  if(n>i) {s3[0]=(char)0; return(0);}
  return(j);
}

Referenced by atofList(), doubleGetWithUnit(), floatGetWithUnit(), iftGetDoubleWithUnit(), intlistAddFromString(), intlistExpandFromString(), parFromIFT(), parReadCSV(), parReadRES(), and tacReadSIF().

◆ strTokenNr()

int strTokenNr	(	const char *	s1,
		const char *	s2 )

extern

Calculate the number of tokens in the string s1. The characters making up the string s2 are the delimiters that determine the token.

See also: strTokenNCpy, strTokenDup

Returns: Returns the nr of tokens.

Author: Vesa Oikonen

Parameters

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

Definition at line 25 of file stringext.c.

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

Referenced by atofList(), doubleGetWithUnit(), floatGetWithUnit(), intlistAddFromString(), intlistExpandFromString(), parFromIFT(), parReadRES(), and tacNameSplit().

◆ strTrimLeft()

char * strTrimLeft	(	char *	s,
		const size_t	t )

extern

Strip string of t characters from the left.

Returns: Returns pointer to the string.

See also: filenameGetExtension

Parameters

s	Pointer to the string to be trimmed. Overwritten.
t	Number of characters to remove.

Definition at line 710 of file stringext.c.

  {
  if(s==NULL || t==0) return(s);
  size_t len=strlen(s); if(len<1) return(s);
  if(t>=len) {s[0]=(char)0; return(s);}
  memmove(s, s+t, 1+len-t);
  return(s);
}

◆ 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 118 of file stringext.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;
}

◆ studynrFromFilename()

int studynrFromFilename	(	const char *	fname,
		char *	studynr,
		int	force )

extern

Extract study number (max MAX_STUDYNR_LEN chars) from given file name.

If file name string includes path, the path is ignored. This function removes initial zeroes from the numerical part, and converts upper case letters to lower case, if necessary.

Returns: 0 if successful.

See also: studynrStandardize, studynrVerify, filenameRmExtensions, filenameRmPath

Parameters

fname	File name, which may include the path.
studynr	Pointer for the resulting study number. Memory must be allocated before calling this function (at least MAX_STUDYNR_LEN+1 characters).
force	Validity of studynr is verified (0) or not verified (1).

Definition at line 79 of file studynr.c.

  {
  if(studynr==NULL) return(1); 
  strcpy(studynr, "");
  if(fname==NULL || strnlen(fname, 2)<2) return(1);
  /* Work with copy of file name */
  char *temp=strdup(fname);
  /* Remove path and extensions */
  filenameRmPath(temp); filenameRmExtensions(temp);
  /* Check if any part of base file name could be a valid study number */
  int ret=0, j;
  size_t i, k, n, len;
  len=strlen(temp); if(len<1) {free(temp); return(2);}
  for(i=1; i<len; i++) if(isalpha(temp[i-1]) && isdigit(temp[i])) {
    j=i-1; while(j>0 && isalpha(temp[j-1])) j--;
    k=i; while(k<len-1 && isdigit(temp[k+1])) k++;
    n=k-j+1; if(n>=2 && n<=10) {
      strlcpy(studynr, temp+j, n+1);
      if(force==1) {ret=1; break;}
      if(studynrStandardize(studynr)==0 && studynrVerify(studynr, 0)==0) {
        ret=2; break;}
    } 
  }
  if(ret) {free(temp); return 0;}
  if(force==0) {free(temp); return 3;}
  /* ok we have to take something */
  for(j=0; j<(int)len; j++) if(isalpha(temp[j]) || isdigit(temp[j])) break;
  for(k=j; k<len-1; k++) if(!isalpha(temp[k+1]) && !isdigit(temp[k+1])) break;
  n=k-j+1; if(n<1) {free(temp); return 4;}
  if(n>10) n=10;
  strlcpy(studynr, temp+j, n+1);
  free(temp);
  return 0;
}

Referenced by tacRead().

◆ studynrStandardize()

int studynrStandardize ( char * s )

extern

Remove initial zeroes in the numerical part of the PET study number, and convert any capital letters to lower-case.

Returns: non-zero in case of failure.

See also: studynrVerify, studynrFromFilename

Parameters

s	Pointer to the study number string to be edited.

Definition at line 47 of file studynr.c.

  {
  if(s==NULL) return(1);
  size_t i, j, len;
  len=strlen(s); if(len<2) return(1);
  if(isdigit((int)s[0])) return(2);
  if(!isdigit((int)s[len-1])) return(3);
  /* Convert to lower-case */
  for(i=0; i<len; i++) {
    if(isdigit((int)s[i])) break;
    s[i]=(char)tolower((int)s[i]);
  }
  /* Remove initial zeroes */
  for(i=1; i<len; i++) if(isdigit((int)s[i])) break;
  for(j=i; j<len; j++) if(s[j]!='0') break;
  if(i==j) return(0); /* no initial zeroes */
  for(;j<=len; i++, j++) s[i]=s[j];
  return(0);
}

Referenced by studynrFromFilename().

◆ studynrVerify()

int studynrVerify	(	const char *	s,
		int	zero_ok )

extern

Check that the argument string is a valid TPC study number containing 1-5 letters (upper or lower case) followed by 1-5 digits.

Author: Vesa Oikonen

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

See also: studynrFromFilename, studynrStandardize

Parameters

s	Pointer to string to test.
zero_ok	Number part of study number may start with zero (1) or may not (0).

Definition at line 22 of file studynr.c.

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

Referenced by studynrFromFilename().

◆ swabip()

void swabip	(	void *	buf,
		int	size )

extern

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

See also: swawbip, swawip, swap, endianLittle

Parameters

buf	Pointer to memory
size	Size of buf in bytes

Definition at line 115 of file endian.c.

  {
  if(buf==NULL) return;
 
  int i;
  unsigned char c;
 
  for(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(), dcmFileWrite(), dcmWriteFileTag(), niftiReadHeader(), and niftiWriteHeader().

◆ swap()

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

extern

Swap the specified short int, int, long int, float, or double from little endian to big endian or vice versa.

Arguments are allowed to overlap.

See also: swawbip, swabip, swawip, endianLittle

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 69 of file endian.c.

  {
  if(from==NULL || to==NULL) return;
 
  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(), and swap().

◆ swap16ip()

void swap16ip	(	void *	buf,
		unsigned long long	size )

extern

In-place big-little endian swapping of an array of 16 bit (2 byte) integers.

See also: swap32ip, swap64ip, endianLittle

Parameters

buf	Pointer to memory.
size	Nr of 16 bit values in the array.

Definition at line 234 of file endian.c.

  {
  if(buf==NULL || size==0) return;
 
  unsigned long long int i;
  unsigned char c, *cptr;
 
  cptr=(unsigned char*)buf;
  for(i=0; i<size; i++, cptr+=2) {
    c=cptr[0]; cptr[0]=cptr[1]; cptr[1]=c;
  }
}

Referenced by imgReadNifti(), niftiReadHeader(), and niftiWriteHeader().

◆ swap32ip()

void swap32ip	(	void *	buf,
		unsigned long long	size )

extern

In-place big-little endian swapping of an array of 32 bit (4 byte) integers or floating point values.

See also: swap16ip, swap64ip, endianLittle

Parameters

buf	Pointer to memory.
size	Nr of 32 bit values in the array.

Definition at line 210 of file endian.c.

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

Referenced by imgReadNifti(), niftiReadHeader(), and niftiWriteHeader().

◆ swap64ip()

void swap64ip	(	void *	buf,
		unsigned long long	size )

extern

In-place big-little endian swapping of an array of 64 bit (8 byte) integers or floating point values.

See also: swap32ip, swap16ip, endianLittle

Parameters

buf	Pointer to memory.
size	Nr of 64 bit values in the array.

Definition at line 184 of file endian.c.

  {
  if(buf==NULL || size==0) return;
 
  unsigned long long int i;
  unsigned char c, *cptr;
 
  cptr=(unsigned char*)buf;
  for(i=0; i<size; i++, cptr+=8) {
    c=cptr[0]; cptr[0]=cptr[7]; cptr[7]=c;
    c=cptr[1]; cptr[1]=cptr[6]; cptr[6]=c;
    c=cptr[2]; cptr[2]=cptr[5]; cptr[5]=c;
    c=cptr[3]; cptr[3]=cptr[4]; cptr[4]=c;
  }
}

Referenced by imgReadNifti(), niftiReadHeader(), and niftiWriteHeader().

◆ swawbip()

void swawbip	(	void *	buf,
		int	size )

extern

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

See also: swawip, swabip, swap, endianLittle

Parameters

buf	Pointer to memory
size	Size of buf in bytes

Definition at line 138 of file endian.c.

  {
  if(buf==NULL) return;
 
  int i;
  unsigned char c, *cptr;
 
  cptr=(unsigned char*)buf;
  for(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(), ecatMListRead(), niftiReadHeader(), and niftiWriteHeader().

◆ swawip()

void swawip	(	void *	buf,
		int	size )

extern

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

See also: swawbip, swabip, swap, endianLittle

Parameters

buf	Pointer to memory
size	Size of buf in bytes

Definition at line 161 of file endian.c.

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

◆ 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.

See also: ctime_r_int, tmDifference, tmAdd

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 443 of file datetime.c.

  {
  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;
}

◆ 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 thread-safe in Windows.

See also: gmtime_r, localtime_r

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

Parameters

tm	Pointer to struct tm.

Definition at line 77 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 tmAdd().

◆ tmAdd()

void tmAdd	(	int	s,
		struct tm *	d )

extern

Add given time in seconds to the date and time.

See also: tmDifference, time_to_tm, strDateTimeAdd, strDateTimeDifference

Parameters

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

Definition at line 500 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
}

Referenced by abssWrite(), and strDateTimeAdd().

◆ tmDifference()

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

extern

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

See also: strDateTimeDifference, tmAdd, time_to_tm

Returns: Returns tm1-tm0 in seconds.

Parameters

tm1	Pointer to tm structure.
tm0	Pointer to tm structure.

Definition at line 486 of file datetime.c.

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

Referenced by strDateTimeDifference().

◆ 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, tpcProgramName, tpcPrintBuild

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 169 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 filename */
  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.

See also: tpcProgramName, tpcPrintUsage

Author: Vesa Oikonen

Parameters

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

Definition at line 339 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.
See also
tpcHtmlUsage, tpcProgramName, tpcPrintBuild

Author
Vesa Oikonen

Parameters

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

Definition at line 114 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: tpcYesNo, statusSet

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 47 of file proginfo.c.

  {
  int n;
 
  char *cptr;
  /* Check that string is option, starting with '-' or '--' */
  if(s==NULL || strnlen(s, 2)<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(strnlen(cptr, 1)<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 && strnlen(cptr, 2)>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, based on user-given string, usually argv[0].

See also: tpcPrintBuild, tpcPrintUsage

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	Buffer length of prname, max n-1 characters are written into it, plus nul char.

Definition at line 406 of file proginfo.c.

  {
  char *tmp;
 
  /* Check the input */
  if(prname==NULL || n<1) return;
  prname[0]=(char)0;
  
  /* Remove path and extension */
  if(strnlen(program, 1024)>0) tmp=strndup(program, 1024); 
  else tmp=strdup("unknown"); 
  filenameRmPath(tmp); filenameRmExtension(tmp);
  /* Copy it if possible */
  n-=strnlen(tmp, 1024); if(n>0) strlcpy(prname, tmp, n); 
  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 */
    strlcat(prname, v, n);
  }
  
  /* Add copyright, if required */
  if(copyright!=0) {
    /* Copy it if possible */
    strlcat(prname, " ", n); 
    strlcat(prname, tpcclib_COPYRIGHT, n); 
  }
  
  return;
}

◆ tpcYesNo()

int tpcYesNo ( const char * s )

extern

Check whether argument string is 'y(es)' or 'n(o)', or 'on' or 'off'.

Returns: Returns 0 when argument means no, 0 when it means no, and -1 if neither.

Author: Vesa Oikonen

See also: atofVerified, atofCheck, atoiCheck

Parameters

s	String to be tested.

Definition at line 459 of file proginfo.c.

  {
  if(s==NULL || strnlen(s, 3)<1) return(-1);
  if(strncasecmp(s, "YES", 1)==0) return(1);
  if(strncasecmp(s, "NO", 1)==0) return(0);
  if(strcasecmp(s, "ON")==0) return(1);
  if(strcasecmp(s, "OFF")==0) return(0);
  return(-1);
}

◆ unitCombination()

int unitCombination	(	const int	u1,
		const int	u2,
		const int	v1,
		const int	v2 )

extern

Check whether the combination of units is an identifiable combinatorial unit.

Returns: enum unit, or 0 (enum UNIT_UNKNOWN) if not supported.

See also: unitConversionFactor, unitIsCombinatorial, unitIdentify, unitName

Parameters

u1	Unit #1 of dividend
u2	Unit #2 of dividend
v1	Unit #1 of divider
v2	Unit #2 of divider

Definition at line 745 of file units.c.

  {
  /* Try if unit table has this combination */
  int i, n=0;
  while(strlen(tpc_unit[n].name)>0) n++;
  for(i=0; i<n; i++) {
    if(!((u1==tpc_unit[i].u1 && u2==tpc_unit[i].u2) || (u1==tpc_unit[i].u2 && u2==tpc_unit[i].u1)))
      continue;
    if(!((v1==tpc_unit[i].v1 && v2==tpc_unit[i].v2) || (v1==tpc_unit[i].v2 && v2==tpc_unit[i].v1)))
      continue;
    return i;
  }
  return(UNIT_UNKNOWN);
}

Referenced by unitDividerMassVolumeConversion(), and unitMultiply().

◆ unitConversionFactor()

double unitConversionFactor	(	const int	u1,
		const int	u2 )

extern

Calculates conversion factor between specified two units, if possible.

Multiply value with factor to get the value in unit2.

Returns: Conversion factor (or NaN, if conversion is not possible).

See also: unitIsCombinatorial, unitIdentify, unitName, unitInverse, unitCombination, unitDividerMassVolumeConversion, unitIsRAConc

Parameters

u1	Enum unit 1.
u2	Enum unit 2.

Definition at line 487 of file units.c.

  {
  //printf("unitConversionFactor(%s, %s)\n", unitName(u1), unitName(u2));
  double cf=nan("");
  if(u1==UNIT_UNKNOWN || u2==UNIT_UNKNOWN) return cf;
  if(u1==u2) return 1.0;
 
  /* If both are time units, then this is easy */
  if(unitIsTime(u1) && unitIsTime(u2)) {
    if(u1==UNIT_MSEC) cf=1.0E-3;
    else if(u1==UNIT_SEC) cf=1.0;
    else if(u1==UNIT_MIN) cf=60.0;
    else if(u1==UNIT_HOUR) cf=3600.0;
    else if(u1==UNIT_DAY) cf=24.0*36000.0;
    else if(u1==UNIT_MONTH) cf=30.0*24.0*36000.0;
    else if(u1==UNIT_YEAR) cf=365.0*30.0*24.0*36000.0;
    if(u2==UNIT_MSEC) cf/=1.0E-3;
    else if(u2==UNIT_SEC) cf/=1.0;
    else if(u2==UNIT_MIN) cf/=60.0;
    else if(u2==UNIT_HOUR) cf/=3600.0;
    else if(u2==UNIT_DAY) cf/=24.0*36000.0;
    else if(u2==UNIT_MONTH) cf/=30.0*24.0*36000.0;
    else if(u2==UNIT_YEAR) cf/=365.0*30.0*24.0*36000.0;
    else cf=nan("");
    return cf;
  }
 
  /* If both are distances */
  if(unitIsDistance(u1) && unitIsDistance(u2)) {
    if(u1==UNIT_UM) cf=1.0E-6;
    else if(u1==UNIT_MM) cf=1.0E-3;
    else if(u1==UNIT_CM) cf=1.0E-2;
    else if(u1==UNIT_M) cf=1.0;
    if(u2==UNIT_UM) cf/=1.0E-6;
    else if(u2==UNIT_MM) cf/=1.0E-3;
    else if(u2==UNIT_CM) cf/=1.0E-2;
    else if(u2==UNIT_M) cf/=1.0;
    else cf=nan("");
    return cf;
  }
 
  /* If both are volumes */
  if(unitIsVolume(u1) && unitIsVolume(u2)) {
    if(u1==UNIT_UL) cf=1.0E-6;
    else if(u1==UNIT_ML) cf=1.0E-3;
    else if(u1==UNIT_DL) cf=1.0E-1;
    else if(u1==UNIT_L) cf=1.0;
    if(u2==UNIT_UL) cf/=1.0E-6;
    else if(u2==UNIT_ML) cf/=1.0E-3;
    else if(u2==UNIT_DL) cf/=1.0E-1;
    else if(u2==UNIT_L) cf/=1.0;
    else cf=nan("");
    return cf;
  }
 
  /* If both are masses */
  if(unitIsMass(u1) && unitIsMass(u2)) {
    if(u1==UNIT_UG) cf=1.0E-6;
    else if(u1==UNIT_MG) cf=1.0E-3;
    else if(u1==UNIT_G) cf=1.0;
    else if(u1==UNIT_100G) cf=1.0E+2;
    else if(u1==UNIT_KG) cf=1.0E+3;
    if(u2==UNIT_UG) cf/=1.0E-6;
    else if(u2==UNIT_MG) cf/=1.0E-3;
    else if(u2==UNIT_G) cf/=1.0;
    else if(u2==UNIT_100G) cf/=1.0E+2;
    else if(u2==UNIT_KG) cf/=1.0E+3;
    else cf=nan("");
    return cf;
  }
 
  /* If both are moles */
  if(unitIsMole(u1) && unitIsMole(u2)) {
    if(u1==UNIT_PMOL)      cf=1.0E-12;
    else if(u1==UNIT_NMOL) cf=1.0E-9;
    else if(u1==UNIT_UMOL) cf=1.0E-6;
    else if(u1==UNIT_MMOL) cf=1.0E-3;
    else if(u1==UNIT_MOL)  cf=1.0;
    if(u2==UNIT_PMOL)      cf/=1.0E-12;
    else if(u2==UNIT_NMOL) cf/=1.0E-9;
    else if(u2==UNIT_UMOL) cf/=1.0E-6;
    else if(u2==UNIT_MMOL) cf/=1.0E-3;
    else if(u2==UNIT_MOL)  cf/=1.0;
    else cf=nan("");
    return cf;
  }
 
  /* If both are radioactivities */
  if(unitIsRadioactivity(u1) && unitIsRadioactivity(u2)) {
    if(u1==UNIT_BQ)       cf=1.0;
    else if(u1==UNIT_KBQ) cf=1.0E+3;
    else if(u1==UNIT_MBQ) cf=1.0E+6;
    else if(u1==UNIT_GBQ) cf=1.0E+9;
    else if(u1==UNIT_NCI) cf=37.0;
    else if(u1==UNIT_UCI) cf=3.7E+4;
    else if(u1==UNIT_MCI) cf=3.7E+7;
    else if(u1==UNIT_CI)  cf=3.7E+10;
    if(u2==UNIT_BQ)       cf/=1.0;
    else if(u2==UNIT_KBQ) cf/=1.0E+3;
    else if(u2==UNIT_MBQ) cf/=1.0E+6;
    else if(u2==UNIT_GBQ) cf/=1.0E+9;
    else if(u2==UNIT_NCI) cf/=37.0;
    else if(u2==UNIT_UCI) cf/=3.7E+4;
    else if(u2==UNIT_MCI) cf/=3.7E+7;
    else if(u2==UNIT_CI)  cf/=3.7E+10;
    else cf=nan("");
    return cf;
  }
  
  /* Fraction vs percentage */
  if(u1==UNIT_PERCENTAGE) {
    if(u2==UNIT_PERCENTAGE) cf=1.0;
    else if(u2==UNIT_UNITLESS) cf=0.01;
    return cf;
  }    
  if(u2==UNIT_PERCENTAGE) {
    if(u1==UNIT_PERCENTAGE) cf=1.0;
    else if(u1==UNIT_UNITLESS) cf=100.0;
    return cf;
  }    
 
  /* Counts */
  if(u1==UNIT_COUNTS) {
    if(u2==UNIT_COUNTS) cf=1.0;
    else if(u2==UNIT_KCOUNTS) cf=0.001;
    return cf;
  }    
  if(u1==UNIT_KCOUNTS) {
    if(u2==UNIT_KCOUNTS) cf=1.0;
    else if(u2==UNIT_COUNTS) cf=1000.0;
    return cf;
  }    
  
  
  /* And last, the combination units, recursively */
  if(!unitIsCombinatorial(u1) || !unitIsCombinatorial(u2)) return cf;
 
  if(tpc_unit[u1].u1==UNIT_UNITLESS && tpc_unit[u2].u1==UNIT_UNITLESS)
    cf=1.0;
  else
    cf=unitConversionFactor(tpc_unit[u1].u1, tpc_unit[u2].u1);
  if(!isnan(cf)) {
    if(tpc_unit[u1].u2!=UNIT_UNITLESS && tpc_unit[u2].u2!=UNIT_UNITLESS)
      cf*=unitConversionFactor(tpc_unit[u1].u2, tpc_unit[u2].u2);
    if(!isnan(cf)) {
      if(tpc_unit[u1].v1!=UNIT_UNITLESS && tpc_unit[u2].v1!=UNIT_UNITLESS)
        cf/=unitConversionFactor(tpc_unit[u1].v1, tpc_unit[u2].v1);
      if(!isnan(cf)) {
        if(tpc_unit[u1].v2!=UNIT_UNITLESS && tpc_unit[u2].v2!=UNIT_UNITLESS)
          cf/=unitConversionFactor(tpc_unit[u1].v2, tpc_unit[u2].v2);
      }
    }
  }
  
  return cf;
}

Referenced by parReadCSV(), tacDecayCorrection(), tacFittime(), tacXMatch(), tacXUnitConvert(), tacYUnitConvert(), and unitConversionFactor().

◆ unitDividendHasRadioactivity()

int unitDividendHasRadioactivity ( int u )

extern

Check whether the unit has radioactivity in the dividend.

Returns: 1 if radioactivity in dividend, otherwise 0.

See also: unitIsRadioactivity, unitIsRAConc

Parameters

u Enum unit

Definition at line 710 of file units.c.

  {
  if(unitIsRadioactivity(u)) return(1);
  if(!unitIsCombinatorial(u)) return(0);
  if(unitIsRadioactivity(tpc_unit[u].u1) || unitIsRadioactivity(tpc_unit[u].u2)) return(1);
  return(0);
}

Referenced by unitIsRAConc().

◆ unitDividerHasMass()

int unitDividerHasMass ( int u )

extern

Check whether the unit has mass in divider.

Returns: 1 if mass found in divider, otherwise 0.

See also: unitDividerHasVolume, unitIsMass, unitInverse, unitDividerMassVolumeConversion

Parameters

u Enum unit

Definition at line 694 of file units.c.

  {
  if(!unitIsCombinatorial(u)) return 0;
  if(unitIsMass(tpc_unit[u].v1)) return 1;
  if(unitIsMass(tpc_unit[u].v2)) return 1;
  return 0;
}

Referenced by tacYUnitMass2Volume(), unitDividerMassVolumeConversion(), and unitIsRAConc().

◆ unitDividerHasVolume()

int unitDividerHasVolume ( int u )

extern

Check whether the unit has volume in divider.

Returns: 1 if volume found in divider, otherwise 0.

See also: unitIsVolume, unitDividerHasMass, unitInverse, unitDividerMassVolumeConversion

Parameters

u Enum unit

Definition at line 678 of file units.c.

  {
  if(!unitIsCombinatorial(u)) return 0;
  if(unitIsVolume(tpc_unit[u].v1)) return 1;
  if(unitIsVolume(tpc_unit[u].v2)) return 1;
  return 0;
}

Referenced by tacYUnitVolume2Mass(), unitDividerMassVolumeConversion(), and unitIsRAConc().

◆ unitDividerMassVolumeConversion()

int unitDividerMassVolumeConversion ( int u )

extern

Convert unit that has either volume or mass unit in divider to corresponding unit (assuming density of 1 g/mL) with either mass or volume in divider.

Returns: enum unit, or 0 (enum UNIT_UNKNOWN) if not supported.

See also: unitConversionFactor, unitDividerHasVolume, unitDividerHasMass

Parameters

u	Unit (enum)

Definition at line 775 of file units.c.

  {
  if(unitDividerHasVolume(u)) {
    int uo=UNIT_UNKNOWN, uu=UNIT_UNKNOWN, uoo=UNIT_UNKNOWN;
    if(unitIsVolume(tpc_unit[u].v1)) {uo=tpc_unit[u].v1; uoo=tpc_unit[u].v2;} 
    else if(unitIsVolume(tpc_unit[u].v2)) {uo=tpc_unit[u].v2; uoo=tpc_unit[u].v1;} 
    if(uo==UNIT_ML) uu=UNIT_G;
    else if(uo==UNIT_DL) uu=UNIT_100G;
    else if(uo==UNIT_L) uu=UNIT_KG;
    else if(uo==UNIT_UL) uu=UNIT_MG;
    if(uu==UNIT_UNKNOWN) return(UNIT_UNKNOWN);
    return(unitCombination(tpc_unit[u].u1, tpc_unit[u].u2, uoo, uu));
  } else if(unitDividerHasMass(u)) {
    int uo=UNIT_UNKNOWN, uu=UNIT_UNKNOWN, uoo=UNIT_UNKNOWN;
    if(unitIsMass(tpc_unit[u].v1)) {uo=tpc_unit[u].v1; uoo=tpc_unit[u].v2;} 
    else if(unitIsMass(tpc_unit[u].v2)) {uo=tpc_unit[u].v2; uoo=tpc_unit[u].v1;} 
    if(uo==UNIT_G) uu=UNIT_ML;
    else if(uo==UNIT_100G) uu=UNIT_DL;
    else if(uo==UNIT_KG) uu=UNIT_L;
    else if(uo==UNIT_MG) uu=UNIT_UL;
    if(uu==UNIT_UNKNOWN) return(UNIT_UNKNOWN);
    return(unitCombination(tpc_unit[u].u1, tpc_unit[u].u2, uoo, uu));
  }
  return(UNIT_UNKNOWN);
}

Referenced by tacYUnitMass2Volume(), and tacYUnitVolume2Mass().

◆ unitIdentify()

int unitIdentify ( const char * s )

extern

Identify the unit based on string representation s of the unit.

Returns: enum unit, or 0 (enum UNIT_UNKNOWN) if not identified.

See also: unitName, unitIdentifyFilename, unitIsTime, unitDividerMassVolumeConversion

Parameters

s	Unit as a string.

Definition at line 162 of file units.c.

  {
  if(s==NULL || strnlen(s, 1)<1) return UNIT_UNKNOWN;
  /* Try if string can be found directly in the table */
  int i, n=0;
  while(strlen(tpc_unit[n].name)>0) n++;
  for(i=0; i<n; i++) {
    if(strcasecmp(tpc_unit[i].name, s)==0) return i;
  }
  /* Unit string is not following TPC standard, lets try something else */
  if(     strcasecmp(s, "ms")==0)               return UNIT_MSEC;
  else if(strcasecmp(s, "s")==0)                return UNIT_SEC;
  else if(strcasecmp(s, "second")==0)           return UNIT_SEC;
  else if(strcasecmp(s, "seconds")==0)          return UNIT_SEC;
  else if(strcasecmp(s, "mins")==0)             return UNIT_MIN;
  else if(strcasecmp(s, "minute")==0)           return UNIT_MIN;
  else if(strcasecmp(s, "minutes")==0)          return UNIT_MIN;
  else if(strcasecmp(s, "hours")==0)            return UNIT_HOUR;
  else if(strcasecmp(s, "days")==0)             return UNIT_DAY;
  else if(strcasecmp(s, "months")==0)           return UNIT_MONTH;
  else if(strcasecmp(s, "years")==0)            return UNIT_YEAR;
  else if(strcasecmp(s, "yrs")==0)              return UNIT_YEAR;
  else if(strcasecmp(s, "microm")==0)           return UNIT_UM;
  else if(strcasecmp(s, "millim")==0)           return UNIT_MM;
  else if(strcasecmp(s, "microl")==0)           return UNIT_UL;
  else if(strcasecmp(s, "mm^3")==0)             return UNIT_UL;
  else if(strcasecmp(s, "cc")==0)               return UNIT_ML;
  else if(strcasecmp(s, "liter")==0)            return UNIT_L;
  else if(strcasecmp(s, "microgram")==0)        return UNIT_UG;
  else if(strcasecmp(s, "milligram")==0)        return UNIT_MG;
  else if(strcasecmp(s, "gram")==0)             return UNIT_G;
  else if(strcasecmp(s, "100 g")==0)            return UNIT_100G;
  else if(strcasecmp(s, "kilogram")==0)         return UNIT_KG;
  else if(strcasecmp(s, "micromol")==0)         return UNIT_UMOL;
  else if(strcasecmp(s, "micromoles")==0)       return UNIT_UMOL;
  else if(strcasecmp(s, "millimol")==0)         return UNIT_MMOL;
  else if(strcasecmp(s, "millimoles")==0)       return UNIT_MMOL;
  else if(strcasecmp(s, "mole")==0)             return UNIT_MOL;
  else if(strcasecmp(s, "moles")==0)            return UNIT_MOL;
  else if(strcasecmp(s, "cnts")==0)             return UNIT_COUNTS;
  else if(strcasecmp(s, "kcnts")==0)            return UNIT_KCOUNTS;
  else if(strcasecmp(s, "becquerels")==0)       return UNIT_BQ;
  else if(strcasecmp(s, "kbecquerels")==0)      return UNIT_KBQ;
  else if(strcasecmp(s, "cnts/s")==0)           return UNIT_CPS;
  else if(strcasecmp(s, "counts/s")==0)         return UNIT_CPS;
  else if(strcasecmp(s, "counts/sec")==0)       return UNIT_CPS;
  else if(strcasecmp(s, "ECAT counts/sec")==0)  return UNIT_CPS;
  else if(strcasecmp(s, "kcounts/s")==0)        return UNIT_KCPS;
  else if(strcasecmp(s, "kcounts/sec")==0)      return UNIT_KCPS;
  else if(strcasecmp(s, "counts/min")==0)       return UNIT_CPM;
  else if(strcasecmp(s, "kcounts/min")==0)      return UNIT_KCPM;
  else if(strcasecmp(s, "Bq/cc")==0)            return UNIT_BQ_PER_ML;
  else if(strcasecmp(s, "Bqcc")==0)             return UNIT_BQ_PER_ML;
  else if(strcasecmp(s, "BqmL")==0)             return UNIT_BQ_PER_ML;
  else if(strcasecmp(s, "kBq/cc")==0)           return UNIT_KBQ_PER_ML;
  else if(strcasecmp(s, "kBqcc")==0)            return UNIT_KBQ_PER_ML;
  else if(strcasecmp(s, "kBqmL")==0)            return UNIT_KBQ_PER_ML;
  else if(strcasecmp(s, "MBq/cc")==0)           return UNIT_MBQ_PER_ML;
  else if(strcasecmp(s, "MBqcc")==0)            return UNIT_MBQ_PER_ML;
  else if(strcasecmp(s, "MBqmL")==0)            return UNIT_MBQ_PER_ML;
  else if(strcasecmp(s, "nCi/cc")==0)           return UNIT_NCI_PER_ML;
  else if(strcasecmp(s, "nCicc")==0)            return UNIT_NCI_PER_ML;
  else if(strcasecmp(s, "nCimL")==0)            return UNIT_NCI_PER_ML;
  else if(strcasecmp(s, "uCi/cc")==0)           return UNIT_UCI_PER_ML;
  else if(strcasecmp(s, "uCicc")==0)            return UNIT_UCI_PER_ML;
  else if(strcasecmp(s, "uCimL")==0)            return UNIT_UCI_PER_ML;
  else if(strcasecmp(s, "integral")==0)         return UNIT_SEC_KBQ_PER_ML;
  else if(strcasecmp(s, "percent")==0)          return UNIT_PERCENTAGE;
  else if(strcasecmp(s, "perc")==0)             return UNIT_PERCENTAGE;
  else if(strcasecmp(s, "cc/cc")==0)            return UNIT_ML_PER_ML;
  else if(strcasecmp(s, "mL/cc")==0)            return UNIT_ML_PER_ML;
  else if(strcasecmp(s, "cc/mL")==0)            return UNIT_ML_PER_ML;
  else if(strcasecmp(s, "g/cc")==0)             return UNIT_G_PER_ML;
  else if(strcasecmp(s, "SUV")==0)              return UNIT_G_PER_ML;
  else if(strcasecmp(s, "SUVbw")==0)            return UNIT_G_PER_ML;
  else if(strcasecmp(s, "SUV_bw")==0)           return UNIT_G_PER_ML;
  else if(strcasecmp(s, "cc/g")==0)             return UNIT_ML_PER_G;
  else if(strcasecmp(s, "1/s")==0)              return UNIT_PER_SEC;
  else if(strcasecmp(s, "s-1")==0)              return UNIT_PER_SEC;
  else if(strcasecmp(s, "min-1")==0)            return UNIT_PER_MIN;
  else if(strcasecmp(s, "h-1")==0)              return UNIT_PER_HOUR;
  else if(strcmp(s, "pM")==0)                   return UNIT_PMOL_PER_L;
  else if(strcmp(s, "nM")==0)                   return UNIT_NMOL_PER_L;
  else if(strcmp(s, "uM")==0)                   return UNIT_UMOL_PER_L;
  else if(strcmp(s, "mM")==0)                   return UNIT_MMOL_PER_L;
  else if(strcasecmp(s, "mL/(sec*mL)")==0)      return UNIT_ML_PER_ML_SEC;
  else if(strcasecmp(s, "mL/sec/mL")==0)        return UNIT_ML_PER_ML_SEC;
  else if(strcasecmp(s, "mL/mL/sec")==0)        return UNIT_ML_PER_ML_SEC;
  else if(strcasecmp(s, "cc/sec/cc")==0)        return UNIT_ML_PER_ML_SEC;
  else if(strcasecmp(s, "cc/cc/sec")==0)        return UNIT_ML_PER_ML_SEC;
  else if(strcasecmp(s, "mL/(sec*cc)")==0)      return UNIT_ML_PER_ML_SEC;
  else if(strcasecmp(s, "mL/(cc*sec)")==0)      return UNIT_ML_PER_ML_SEC;
  else if(strcasecmp(s, "cc/(cc*sec)")==0)      return UNIT_ML_PER_ML_SEC;
  else if(strcasecmp(s, "mL/(min*mL)")==0)      return UNIT_ML_PER_ML_MIN;
  else if(strcasecmp(s, "mL/min/mL")==0)        return UNIT_ML_PER_ML_MIN;
  else if(strcasecmp(s, "mL/mL/min")==0)        return UNIT_ML_PER_ML_MIN;
  else if(strcasecmp(s, "cc/min/cc")==0)        return UNIT_ML_PER_ML_MIN;
  else if(strcasecmp(s, "cc/cc/min")==0)        return UNIT_ML_PER_ML_MIN;
  else if(strcasecmp(s, "mL/(min*g)")==0)       return UNIT_ML_PER_G_MIN;
  else if(strcasecmp(s, "mL/(min*cc)")==0)      return UNIT_ML_PER_ML_MIN;
  else if(strcasecmp(s, "mL/(cc*min)")==0)      return UNIT_ML_PER_ML_MIN;
  else if(strcasecmp(s, "cc/(cc*min)")==0)      return UNIT_ML_PER_ML_MIN;
  else if(strcasecmp(s, "mL/(min*dL)")==0)      return UNIT_ML_PER_DL_MIN;
  else if(strcasecmp(s, "mL/min/dL")==0)        return UNIT_ML_PER_DL_MIN;
  else if(strcasecmp(s, "mL/dL/min")==0)        return UNIT_ML_PER_DL_MIN;
  else if(strcasecmp(s, "umol/(min*mL)")==0)    return UNIT_UMOL_PER_ML_MIN;
  else if(strcasecmp(s, "umol/(min*dL)")==0)    return UNIT_UMOL_PER_DL_MIN;
  else if(strcasecmp(s, "umol/dL/min)")==0)     return UNIT_UMOL_PER_DL_MIN;
  else if(strcasecmp(s, "umol/(min*g)")==0)     return UNIT_UMOL_PER_G_MIN;
  else if(strcasecmp(s, "umol/(min*100g)")==0)  return UNIT_UMOL_PER_100G_MIN;
  else if(strcasecmp(s, "mmol/(min*mL)")==0)    return UNIT_MMOL_PER_ML_MIN;
  else if(strcasecmp(s, "mmol/(min*dL)")==0)    return UNIT_MMOL_PER_DL_MIN;
  else if(strcasecmp(s, "mmol/dL/min)")==0)     return UNIT_MMOL_PER_DL_MIN;
  else if(strcasecmp(s, "mmol/(min*g)")==0)     return UNIT_MMOL_PER_G_MIN;
  else if(strcasecmp(s, "mmol/(min*100g)")==0)  return UNIT_MMOL_PER_100G_MIN;
  else if(strcasecmp(s, "% ID")==0)             return UNIT_PID;
  else if(strcasecmp(s, "%i.d.")==0)            return UNIT_PID;
  else if(strcasecmp(s, "PID")==0)              return UNIT_PID;
  else if(strcasecmp(s, "% ID/mL")==0)          return UNIT_PID_PER_ML;
  else if(strcasecmp(s, "%i.d./mL")==0)         return UNIT_PID_PER_ML;
  else if(strcasecmp(s, "PID/mL")==0)           return UNIT_PID_PER_ML;
  else if(strcasecmp(s, "%ID/cc")==0)           return UNIT_PID_PER_ML;
  else if(strcasecmp(s, "% ID/cc")==0)          return UNIT_PID_PER_ML;
  else if(strcasecmp(s, "PID/cc")==0)           return UNIT_PID_PER_ML;
  else if(strcasecmp(s, "% ID/L")==0)           return UNIT_PID_PER_L;
  else if(strcasecmp(s, "%i.d./L")==0)          return UNIT_PID_PER_L;
  else if(strcasecmp(s, "PID/L")==0)            return UNIT_PID_PER_L;
  else if(strcasecmp(s, "% ID/g")==0)           return UNIT_PID_PER_G;
  else if(strcasecmp(s, "%i.d./g")==0)          return UNIT_PID_PER_G;
  else if(strcasecmp(s, "PID/g")==0)            return UNIT_PID_PER_G;
  else if(strcasecmp(s, "% ID/kg")==0)          return UNIT_PID_PER_KG;
  else if(strcasecmp(s, "%i.d./kg")==0)         return UNIT_PID_PER_KG;
  else if(strcasecmp(s, "PID/kg")==0)           return UNIT_PID_PER_KG;
  else if(strcasecmp(s, "Hounsfield Unit")==0)  return UNIT_HU;
  else if(strcasecmp(s, "HU")==0)               return UNIT_HU;
  else if(strcasecmp(s, "1/1")==0)              return UNIT_UNITLESS;
 
  return UNIT_UNKNOWN;
}

Referenced by doubleGetWithUnit(), floatGetWithUnit(), iftGetDoubleWithUnit(), imgGetNiftiHeader(), imgReadDICOM(), parFromIFT(), parReadCSV(), parReadFIT(), parReadRES(), parWriteCSV(), tacAllocateWithPAR(), tacGetHeaderTimeunit(), tacGetHeaderUnit(), tacRead4DM(), tacReadCarimasTxt(), tacReadCSV(), tacReadDFT(), tacReadHRRTHC(), tacReadInveonCSV(), and tacReadPMOD().

◆ unitIdentifyFilename()

int unitIdentifyFilename ( const char * s )

extern

Identify calibration unit based on file name.

Returns: enum unit, or 0 (enum UNIT_UNKNOWN) if not identified.

See also: unitName, unitIdentify, unitDividerHasVolume, unitIsRAConc

Parameters

s	File name.

Definition at line 311 of file units.c.

  {
  if(s==NULL || strnlen(s, 3)<3) return UNIT_UNKNOWN;
  char *cptr;
  for(int i=0; i<2; i++) {
    if(i==0) /* On the first run, look in the extension */
      {cptr=strrchr((char*)s, '.'); if(cptr==NULL) {cptr=(char*)s; i++;}}
    else /* Then, look into whole filename */
      cptr=(char*)s;
    if(strcasestr(cptr, "KBQ")) return UNIT_KBQ_PER_ML;
    if(strcasestr(cptr, "MBQ")) return UNIT_MBQ_PER_ML;
    if(strcasestr(cptr, "BQ")) return UNIT_BQ_PER_ML;
    if(strcasestr(cptr, "NCI")) return UNIT_NCI_PER_ML;
    if(strcasestr(cptr, "KCPS")) return UNIT_KCPS;
    if(strcasestr(cptr, "CPS")) return UNIT_CPS;
    if(strcasestr(cptr, "SUV")) return UNIT_G_PER_ML;
  }
  return UNIT_UNKNOWN;
}

Referenced by tacRead().

◆ unitInverse()

int unitInverse ( int u )

extern

Inverse of given unit.

Returns: enum unit, or 0 (enum UNIT_UNKNOWN) if not supported.

See also: unitConversionFactor, unitIsCombinatorial, unitIsTime

Parameters

u	Unit (enum)

Definition at line 654 of file units.c.

  {
  if(u==UNIT_UNKNOWN) return(UNIT_UNKNOWN);
  else if(u==UNIT_UNITLESS) return(UNIT_UNITLESS);
  else if(u==UNIT_SEC) return(UNIT_PER_SEC);
  else if(u==UNIT_PER_SEC) return(UNIT_SEC);
  else if(u==UNIT_MIN) return(UNIT_PER_MIN);
  else if(u==UNIT_PER_MIN) return(UNIT_MIN);
  else if(u==UNIT_HOUR) return(UNIT_PER_HOUR);
  else if(u==UNIT_PER_HOUR) return(UNIT_HOUR);
  else if(u==UNIT_ML_PER_ML) return(UNIT_ML_PER_ML);
  else if(u==UNIT_G_PER_ML) return(UNIT_ML_PER_G);
  else if(u==UNIT_ML_PER_G) return(UNIT_G_PER_ML);
  else return(UNIT_UNKNOWN);
}

◆ unitIsCombinatorial()

int unitIsCombinatorial ( int u )

extern

Check whether the unit is a combinatorial unit.

Returns: 1 if combination unit, otherwise 0.

See also: unitConversionFactor, unitIsRadioactivity, unitIsVolume, unitIsMole, unitCombination

Parameters

u Enum unit

Definition at line 468 of file units.c.

  {
  if(u<1 || u>=UNIT_LAST) return 0;
  if(tpc_unit[u].u2==UNIT_UNITLESS && tpc_unit[u].v1==UNIT_UNITLESS &&
     tpc_unit[u].v2==UNIT_UNITLESS) return 0;
  return 1;
}

Referenced by unitConversionFactor(), unitDividendHasRadioactivity(), unitDividerHasMass(), unitDividerHasVolume(), and unitIsRAConc().

◆ unitIsDistance()

int unitIsDistance ( int u )

extern

Check whether the unit is a unit of distance.

Returns: 1 if distance unit, otherwise 0.

See also: unitIsTime, unitIdentify, unitName

Parameters

u	Enum unit.

Definition at line 339 of file units.c.

  {
  switch(u) {
    case UNIT_UM:
    case UNIT_MM:
    case UNIT_CM:
    case UNIT_M:
      return 1;  
  }
  return 0;
}

Referenced by tacWriteDFT(), and unitConversionFactor().

◆ unitIsMass()

int unitIsMass ( int u )

extern

Check whether the unit is a unit of mass.

Returns: 1 if mass unit, otherwise 0.

See also: unitDividerHasMass, unitIsVolume, unitIsMole

Parameters

u	Enum unit.

Definition at line 402 of file units.c.

  {
  switch(u) {
    case UNIT_UG:
    case UNIT_MG:
    case UNIT_G:
    case UNIT_100G:
    case UNIT_KG:
      return 1;    
  }
  return 0;
}

Referenced by unitConversionFactor(), unitDividerHasMass(), unitDividerMassVolumeConversion(), and unitIsRAConc().

◆ unitIsMole()

int unitIsMole ( int u )

extern

Check whether the unit is a unit of chemical mass.

Returns: 1 if mole unit, otherwise 0.

See also: unitIsRadioactivity, unitIsMass

Parameters

u	Enum unit.

Definition at line 423 of file units.c.

  {
  switch(u) {
    case UNIT_PMOL:
    case UNIT_NMOL:
    case UNIT_UMOL:
    case UNIT_MMOL:
    case UNIT_MOL:
      return 1;    
  }
  return 0;
}

Referenced by unitConversionFactor().

◆ unitIsRAConc()

int unitIsRAConc ( int u )

extern

Check whether the unit is a known unit of radioactivity per mass or volume.

Returns: 1 if radioactivity concentration unit, otherwise 0.

See also: unitDividendHasRadioactivity, unitDividerHasVolume, unitDividerHasMass

Parameters

u Enum unit

Definition at line 726 of file units.c.

  {
  if(!unitIsCombinatorial(u)) return(0);
  if(!unitDividendHasRadioactivity(u)) return(0);
  if(!unitDividerHasMass(u) && !unitDividerHasVolume(u)) return(0);
  if(!unitIsRadioactivity(tpc_unit[u].u1)) return(0);
  if(!unitIsMass(tpc_unit[u].v1) && !unitIsVolume(tpc_unit[u].v1)) return(0);
  if(tpc_unit[u].u2!=UNIT_UNITLESS || tpc_unit[u].v2!=UNIT_UNITLESS) return(0);
  return(1);
}

Referenced by tacSetWeights().

◆ unitIsRadioactivity()

int unitIsRadioactivity ( int u )

extern

Check whether the unit is a unit of radioactivity.

Returns: 1 if radioactivity unit, otherwise 0.

See also: unitIsRAConc, unitIsMole, unitIsMass

Parameters

u Enum unit

Definition at line 444 of file units.c.

  {
  switch(u) {
    case UNIT_BQ:
    case UNIT_KBQ:
    case UNIT_MBQ:
    case UNIT_GBQ:
    case UNIT_NCI:
    case UNIT_UCI:
    case UNIT_MCI:
    case UNIT_CI:
      return 1;    
  }
  return 0;
}

Referenced by tacSetWeights(), unitConversionFactor(), unitDividendHasRadioactivity(), and unitIsRAConc().

◆ unitIsTime()

int unitIsTime ( int u )

extern

Check whether the unit is a unit of time.

Returns: 1 if time unit, otherwise 0.

See also: unitIsDistance, unitConversionFactor, unitInverse, unitIsRAConc, unitIdentify

Parameters

u	Enum unit.

Definition at line 359 of file units.c.

  {
  switch(u) {
    case UNIT_MSEC:
    case UNIT_SEC:
    case UNIT_MIN:
    case UNIT_HOUR:
    case UNIT_DAY:
    case UNIT_MONTH:
    case UNIT_YEAR:
      return 1;    
  }
  return 0;
}

Referenced by bfm1TCM(), parWriteCSV(), tacDecayCorrection(), tacPlotFitSVG(), tacPlotLineSVG(), tacSetWeights(), and unitConversionFactor().

◆ unitIsVolume()

int unitIsVolume ( int u )

extern

Check whether the unit is a unit of volume.

Returns: 1 if volume unit, otherwise 0.

See also: unitDividerHasVolume, unitIsMass, unitName

Parameters

u	Enum unit.

Definition at line 382 of file units.c.

  {
  switch(u) {
    case UNIT_UL:
    case UNIT_ML:
    case UNIT_DL:
    case UNIT_L:
      return 1;    
  }
  return 0;
}

Referenced by unitConversionFactor(), unitDividerHasVolume(), unitDividerMassVolumeConversion(), and unitIsRAConc().

◆ unitMultiply()

int unitMultiply	(	int	ua,
		int	ub )

extern

Multiply two units.

Returns: enum unit, or 0 (enum UNIT_UNKNOWN) if combination is not supported.

See also: unitConversionFactor, unitCombination

Parameters

ua	Unit A
ub	Unit B

Definition at line 809 of file units.c.

  {
  /* Verify validity of input units */
  if(ua<1 || ub<1) return(UNIT_UNKNOWN);
  int n=0; while(strlen(tpc_unit[n].name)>0) n++;
  if(ua>n-1 || ub>n-1) return(UNIT_UNKNOWN);
 
  /* Collect actual components of units */
  int u[4], v[4], nu=0, nv=0;
  for(int i=0; i<4; i++) u[i]=v[i]=UNIT_UNITLESS;
  if(tpc_unit[ua].u1!=UNIT_UNITLESS) u[nu++]=tpc_unit[ua].u1;
  if(tpc_unit[ua].u2!=UNIT_UNITLESS) u[nu++]=tpc_unit[ua].u2;
  if(tpc_unit[ua].v1!=UNIT_UNITLESS) v[nv++]=tpc_unit[ua].v1;
  if(tpc_unit[ua].v2!=UNIT_UNITLESS) v[nv++]=tpc_unit[ua].v2;
  if(tpc_unit[ub].u1!=UNIT_UNITLESS) u[nu++]=tpc_unit[ub].u1;
  if(tpc_unit[ub].u2!=UNIT_UNITLESS) u[nu++]=tpc_unit[ub].u2;
  if(tpc_unit[ub].v1!=UNIT_UNITLESS) v[nv++]=tpc_unit[ub].v1;
  if(tpc_unit[ub].v2!=UNIT_UNITLESS) v[nv++]=tpc_unit[ub].v2;
#if(0)
  printf("units: (");
  for(int i=0; i<nu; i++) printf(" %s", unitName(u[i]));
  printf(" ) / (");
  for(int j=0; j<nv; j++) printf(" %s", unitName(v[j]));
  printf(" )\n");
#endif
 
  /* Cancel out units if possible */
  int cn=0;
  for(int i=0; i<nu; i++) for(int j=0; j<nv; j++) if(u[i]==v[j]) {u[i]=v[i]=UNIT_UNITLESS; cn++;}
  if(cn>0) {
#if(0)
    printf("units after cancelling: (");
    for(int i=0; i<nu; i++) printf(" %s", unitName(u[i]));
    printf(" ) / (");
    for(int j=0; j<nv; j++) printf(" %s", unitName(v[j]));
    printf(" )\n");
#endif
    for(int i=0; i<nu-1; i++) if(u[i]==UNIT_UNITLESS)
      for(int j=i+1; j<nu; i++) if(u[j]!=UNIT_UNITLESS) {
        int s=u[i]; u[i]=u[j]; u[j]=s;
      }
    for(int i=0; i<nv-1; i++) if(v[i]==UNIT_UNITLESS)
      for(int j=i+1; j<nv; i++) if(v[j]!=UNIT_UNITLESS) {
        int s=v[i]; v[i]=v[j]; v[j]=s;
      }
    nu-=cn; nv-=cn;
#if(0)
    printf("units after removing gaps: (");
    for(int i=0; i<nu; i++) printf(" %s", unitName(u[i]));
    printf(" ) / (");
    for(int j=0; j<nv; j++) printf(" %s", unitName(v[j]));
    printf(" )\n");
#endif
  }
 
  /* Is the new combination a valid unit? */
  if(nu>2 || nv>2) return(UNIT_UNKNOWN);
  return(unitCombination(u[0], u[1], v[0], v[1]));
}

◆ unitName()

char * unitName ( int unit_code )

extern

Get the string representation of the unit specified with its unit_code.

Returns: pointer to string representation of the unit, "unknown" if not identified.

See also: unitIdentify, unitIsTime, unitCombination

Parameters

unit_code unit_code as enum or the index of unit in the table.

Definition at line 143 of file units.c.

  {
  int n=0;
  
  while(strlen(tpc_unit[n].name)>0) n++;
  if(unit_code<1 || unit_code>n-1)
    return tpc_unit[UNIT_UNKNOWN].name;
  else
    return(tpc_unit[unit_code].name);
}

Referenced by imgCompareUnit(), imgContents(), imgFillOHeader(), imgReadDICOM(), imgSetNiftiHeader(), parAllocateWithTAC(), parCompareParameterUnits(), parExamplePerfectBolus(), parExampleTTACs(), parToIFT(), parWriteCSV(), parWriteRES(), parWriteXML(), tacCompareUnit(), tacPlotFitSVG(), tacPlotLineSVG(), tacRead(), tacRead4DM(), tacSetHeaderTimeunit(), tacSetHeaderUnit(), tacWriteCSV(), tacWriteDFT(), tacWritePMOD(), tacWriteSheetIntoXML(), tacWriteSimple(), tacWriteXML(), tacXUnitConvert(), tacYUnitConvert(), and unitMultiply().

Data Structures

Macros

Enumerations

Functions

Detailed Description

Macro Definition Documentation

◆ MAX_PARNAME_LEN

◆ MAX_STUDYNR_LEN

◆ MAX_TACNAME_LEN

◆ MAX_UNITS_LEN

Enumeration Type Documentation

◆ anonymous enum

◆ tactype

◆ tpcerror

◆ unit

◆ weights

Function Documentation

◆ asciiCommentLine()

◆ asciiFileRead()

◆ asciiFileSize()

◆ atofCheck()

◆ atofList()

◆ atofVerified()

◆ atoiCheck()

◆ ctime_r_int()

◆ doubleAbsMaxIndex()

◆ doubleAbsMinIndex()

◆ doubleArrayMatch()

◆ doubleCopy()

◆ doubleCopyFinite()

◆ doubleCSpanPositives()

◆ doubleGEIndex()

◆ doubleGetWithUnit()

◆ doubleGTIndex()

◆ doubleMachEps()

◆ doubleMatch()

◆ doubleMatchRel()

◆ doubleMaxIndex()

◆ doubleMean()

◆ doubleMinIndex()

◆ doubleNaNs()

◆ doubleNonzeroes()

◆ doubleRange()

◆ doubleSpanPositives()

◆ doubleSum()

◆ doubleWMean()

◆ endianLittle()

◆ endianness()

◆ errorMsg()

◆ factorial()

◆ filenameGetExtension()

◆ filenameGetExtensions()

◆ filenameRmExtension()

◆ filenameRmExtensions()

◆ filenameRmFile()

◆ filenameRmPath()

◆ floatCopy()

◆ floatCSpanPositives()

◆ floatGetWithUnit()

◆ floatMachEps()

◆ floatMatch()

◆ floatMatchRel()

◆ floatMaxIndex()

◆ floatMean()

◆ floatNonzeroes()

◆ floatSpanPositives()

◆ floatSum()

◆ fncasematch()

◆ fnmatch()

◆ gmtime_r()

◆ intlistAdd()

◆ intlistAddFromString()

◆ intlistExpandFromString()

◆ intlistFree()

◆ intlistInit()

◆ intlistSort()

◆ inverfc()

◆ lfactorial()

◆ localtime_r()

◆ roinameAddField()