152 #include <config_auto.h>
157 #include "allheaders.h"
161 static const l_uint32 MaxDoubleArraySize = 100000000;
162 static const l_uint32 MaxPtrArraySize = 1000000;
183 if (n <= 0 || n > MaxDoubleArraySize)
187 if ((da->
array = (l_float64 *)LEPT_CALLOC(n,
sizeof(l_float64))) == NULL) {
189 return (
L_DNA *)ERROR_PTR(
"double array not made", __func__, NULL);
225 return (
L_DNA *)ERROR_PTR(
"iarray not defined", __func__, NULL);
227 return (
L_DNA *)ERROR_PTR(
"size must be > 0", __func__, NULL);
230 for (i = 0; i < size; i++)
261 return (
L_DNA *)ERROR_PTR(
"darray not defined", __func__, NULL);
263 return (
L_DNA *)ERROR_PTR(
"size must be > 0", __func__, NULL);
265 return (
L_DNA *)ERROR_PTR(
"invalid copyflag", __func__, NULL);
273 for (i = 0; i < size; i++)
299 return (
L_DNA *)ERROR_PTR(
"da not made", __func__, NULL);
301 for (i = 0; i < size; i++) {
302 val = startval + i * increment;
328 L_WARNING(
"ptr address is NULL\n", __func__);
332 if ((da = *pda) == NULL)
338 LEPT_FREE(da->
array);
363 return (
L_DNA *)ERROR_PTR(
"da not defined", __func__, NULL);
366 return (
L_DNA *)ERROR_PTR(
"dac not made", __func__, NULL);
370 for (i = 0; i < da->
n; i++)
387 return (
L_DNA *)ERROR_PTR(
"da not defined", __func__, NULL);
411 return ERROR_INT(
"da not defined", __func__, 1);
436 return ERROR_INT(
"da not defined", __func__, 1);
441 return ERROR_INT(
"extension failed", __func__, 1);
464 size_t oldsize, newsize;
467 return ERROR_INT(
"da not defined", __func__, 1);
468 if (da->
nalloc > MaxDoubleArraySize)
469 return ERROR_INT(
"da at maximum size; can't extend", __func__, 1);
470 oldsize = da->
nalloc *
sizeof(l_float64);
471 if (da->
nalloc > MaxDoubleArraySize / 2) {
472 newsize = MaxDoubleArraySize *
sizeof(l_float64);
473 da->
nalloc = MaxDoubleArraySize;
475 newsize = 2 * oldsize;
479 oldsize, newsize)) == NULL)
480 return ERROR_INT(
"new ptr array not returned", __func__, 1);
511 return ERROR_INT(
"da not defined", __func__, 1);
513 if (index < 0 || index > n) {
514 L_ERROR(
"index %d not in [0,...,%d]\n", __func__, index, n);
520 return ERROR_INT(
"extension failed", __func__, 1);
522 for (i = n; i > index; i--)
524 da->
array[index] = val;
551 return ERROR_INT(
"da not defined", __func__, 1);
553 if (index < 0 || index >= n) {
554 L_ERROR(
"index %d not in [0,...,%d]\n", __func__, index, n - 1);
558 for (i = index + 1; i < n; i++)
581 return ERROR_INT(
"da not defined", __func__, 1);
583 if (index < 0 || index >= n) {
584 L_ERROR(
"index %d not in [0,...,%d]\n", __func__, index, n - 1);
588 da->
array[index] = val;
606 return ERROR_INT(
"da not defined", __func__, 0);
632 return ERROR_INT(
"da not defined", __func__, 1);
633 if (newcount > da->
nalloc) {
635 sizeof(l_float64) * da->
nalloc,
636 sizeof(l_float64) * newcount)) == NULL)
637 return ERROR_INT(
"new ptr array not returned", __func__, 1);
665 return ERROR_INT(
"&val not defined", __func__, 1);
668 return ERROR_INT(
"da not defined", __func__, 1);
670 if (index < 0 || index >= da->
n)
671 return ERROR_INT(
"index not valid", __func__, 1);
673 *pval = da->
array[index];
700 return ERROR_INT(
"&ival not defined", __func__, 1);
703 return ERROR_INT(
"da not defined", __func__, 1);
705 if (index < 0 || index >= da->
n)
706 return ERROR_INT(
"index not valid", __func__, 1);
708 val = da->
array[index];
709 *pival = (l_int32)(val + L_SIGN(val) * 0.5);
728 return ERROR_INT(
"da not defined", __func__, 1);
729 if (index < 0 || index >= da->
n)
730 return ERROR_INT(
"index not valid", __func__, 1);
732 da->
array[index] = val;
751 return ERROR_INT(
"da not defined", __func__, 1);
752 if (index < 0 || index >= da->
n)
753 return ERROR_INT(
"index not valid", __func__, 1);
755 da->
array[index] += diff;
786 return (l_int32 *)ERROR_PTR(
"da not defined", __func__, NULL);
789 if ((array = (l_int32 *)LEPT_CALLOC(n,
sizeof(l_int32))) == NULL)
790 return (l_int32 *)ERROR_PTR(
"array not made", __func__, NULL);
791 for (i = 0; i < n; i++) {
829 return (l_float64 *)ERROR_PTR(
"da not defined", __func__, NULL);
835 if ((array = (l_float64 *)LEPT_CALLOC(n,
sizeof(l_float64))) == NULL)
836 return (l_float64 *)ERROR_PTR(
"array not made", __func__, NULL);
837 for (i = 0; i < n; i++)
838 array[i] = da->
array[i];
858 if (pstartx) *pstartx = 0.0;
859 if (pdelx) *pdelx = 1.0;
860 if (!pstartx && !pdelx)
861 return ERROR_INT(
"neither &startx nor &delx are defined", __func__, 1);
863 return ERROR_INT(
"da not defined", __func__, 1);
865 if (pstartx) *pstartx = da->
startx;
866 if (pdelx) *pdelx = da->
delx;
887 return ERROR_INT(
"da not defined", __func__, 1);
906 l_float64 start, binsize;
909 return ERROR_INT(
"das and dad not both defined", __func__, 1);
933 return (
L_DNA *)ERROR_PTR(
"filename not defined", __func__, NULL);
936 return (
L_DNA *)ERROR_PTR(
"stream not opened", __func__, NULL);
940 return (
L_DNA *)ERROR_PTR(
"da not read", __func__, NULL);
960 l_int32 i, n, index, ret, version;
961 l_float64 val, startx, delx;
965 return (
L_DNA *)ERROR_PTR(
"stream not defined", __func__, NULL);
967 ret = fscanf(fp,
"\nL_Dna Version %d\n", &version);
969 return (
L_DNA *)ERROR_PTR(
"not a l_dna file", __func__, NULL);
971 return (
L_DNA *)ERROR_PTR(
"invalid l_dna version", __func__, NULL);
972 if (fscanf(fp,
"Number of numbers = %d\n", &n) != 1)
973 return (
L_DNA *)ERROR_PTR(
"invalid number of numbers", __func__, NULL);
975 return (
L_DNA *)ERROR_PTR(
"num doubles < 0", __func__, NULL);
976 if (n > MaxDoubleArraySize)
977 return (
L_DNA *)ERROR_PTR(
"too many doubles", __func__, NULL);
978 if (n == 0) L_INFO(
"the dna is empty\n", __func__);
981 return (
L_DNA *)ERROR_PTR(
"da not made", __func__, NULL);
982 for (i = 0; i < n; i++) {
983 if (fscanf(fp,
" [%d] = %lf\n", &index, &val) != 2) {
985 return (
L_DNA *)ERROR_PTR(
"bad input data", __func__, NULL);
991 if (fscanf(fp,
"startx = %lf, delx = %lf\n", &startx, &delx) == 2)
1012 return (
L_DNA *)ERROR_PTR(
"data not defined", __func__, NULL);
1014 return (
L_DNA *)ERROR_PTR(
"stream not opened", __func__, NULL);
1018 if (!da) L_ERROR(
"dna not read\n", __func__);
1038 return ERROR_INT(
"filename not defined", __func__, 1);
1040 return ERROR_INT(
"da not defined", __func__, 1);
1043 return ERROR_INT(
"stream not opened", __func__, 1);
1047 return ERROR_INT(
"da not written to stream", __func__, 1);
1064 l_float64 startx, delx;
1067 return ERROR_INT(
"da not defined", __func__, 1);
1073 fprintf(fp,
"Number of numbers = %d\n", n);
1074 for (i = 0; i < n; i++)
1075 fprintf(fp,
" [%d] = %f\n", i, da->
array[i]);
1080 if (startx != 0.0 || delx != 1.0)
1081 fprintf(fp,
"startx = %f, delx = %f\n", startx, delx);
1097 l_float64 startx, delx;
1100 return ERROR_INT(
"da not defined", __func__, 1);
1105 for (i = 0; i < n; i++)
1111 if (startx != 0.0 || delx != 1.0)
1112 lept_stderr(
"startx = %f, delx = %f\n", startx, delx);
1139 if (pdata) *pdata = NULL;
1140 if (psize) *psize = 0;
1142 return ERROR_INT(
"&data not defined", __func__, 1);
1144 return ERROR_INT(
"&size not defined", __func__, 1);
1146 return ERROR_INT(
"da not defined", __func__, 1);
1149 if ((fp = open_memstream((
char **)pdata, psize)) == NULL)
1150 return ERROR_INT(
"stream not opened", __func__, 1);
1154 *psize = *psize - 1;
1156 L_INFO(
"work-around: writing to a temp file\n", __func__);
1159 return ERROR_INT(
"tmpfile stream not opened", __func__, 1);
1161 if ((fp = tmpfile()) == NULL)
1162 return ERROR_INT(
"tmpfile stream not opened", __func__, 1);
1188 if (n <= 0 || n > MaxPtrArraySize)
1192 if ((daa->
dna = (
L_DNA **)LEPT_CALLOC(n,
sizeof(
L_DNA *))) == NULL) {
1194 return (
L_DNAA *)ERROR_PTR(
"l_dna ptr array not made", __func__, NULL);
1226 for (i = 0; i < nptr; i++) {
1255 return ERROR_INT(
"daa not defined", __func__, 1);
1258 for (i = n - 1; i >= 0; i--) {
1287 L_WARNING(
"ptr address is NULL!\n", __func__);
1291 if ((daa = *pdaa) == NULL)
1294 for (i = 0; i < daa->
n; i++)
1296 LEPT_FREE(daa->
dna);
1322 return ERROR_INT(
"daa not defined", __func__, 1);
1324 return ERROR_INT(
"da not defined", __func__, 1);
1328 }
else if (copyflag ==
L_COPY) {
1330 return ERROR_INT(
"dac not made", __func__, 1);
1331 }
else if (copyflag ==
L_CLONE) {
1334 return ERROR_INT(
"invalid copyflag", __func__, 1);
1342 return ERROR_INT(
"extension failed", __func__, 1);
1366 size_t oldsize, newsize;
1369 return ERROR_INT(
"daa not defined", __func__, 1);
1370 if (daa->
nalloc > MaxPtrArraySize)
1371 return ERROR_INT(
"daa has too many ptrs", __func__, 1);
1373 newsize = 2 * oldsize;
1374 if (newsize > 8 * MaxPtrArraySize)
1375 return ERROR_INT(
"newsize > 8 MB; too large", __func__, 1);
1378 oldsize, newsize)) == NULL)
1379 return ERROR_INT(
"new ptr array not returned", __func__, 1);
1399 return ERROR_INT(
"daa not defined", __func__, 0);
1416 return ERROR_INT(
"daa not defined", __func__, 0);
1417 if (index < 0 || index >= daa->
n)
1418 return ERROR_INT(
"invalid index into daa", __func__, 0);
1437 return ERROR_INT(
"daa not defined", __func__, 0);
1440 for (sum = 0, i = 0; i < n; i++) {
1464 return (
L_DNA *)ERROR_PTR(
"daa not defined", __func__, NULL);
1465 if (index < 0 || index >= daa->
n)
1466 return (
L_DNA *)ERROR_PTR(
"index not valid", __func__, NULL);
1468 if (accessflag ==
L_COPY)
1470 else if (accessflag ==
L_CLONE)
1473 return (
L_DNA *)ERROR_PTR(
"invalid accessflag", __func__, NULL);
1500 return ERROR_INT(
"daa not defined", __func__, 1);
1502 return ERROR_INT(
"da not defined", __func__, 1);
1504 if (index < 0 || index >= n)
1505 return ERROR_INT(
"index not valid", __func__, 1);
1508 daa->
dna[index] = da;
1532 return ERROR_INT(
"&val not defined", __func__, 1);
1535 return ERROR_INT(
"daa not defined", __func__, 1);
1537 if (i < 0 || i >= n)
1538 return ERROR_INT(
"invalid index into daa", __func__, 1);
1540 if (j < 0 || j >= da->
n)
1541 return ERROR_INT(
"invalid index into da", __func__, 1);
1542 *pval = da->
array[j];
1569 return ERROR_INT(
"daa not defined", __func__, 1);
1571 if (index < 0 || index >= n)
1572 return ERROR_INT(
"invalid index in daa", __func__, 1);
1597 return (
L_DNAA *)ERROR_PTR(
"filename not defined", __func__, NULL);
1600 return (
L_DNAA *)ERROR_PTR(
"stream not opened", __func__, NULL);
1604 return (
L_DNAA *)ERROR_PTR(
"daa not read", __func__, NULL);
1623 l_int32 i, n, index, ret, version;
1628 return (
L_DNAA *)ERROR_PTR(
"stream not defined", __func__, NULL);
1630 ret = fscanf(fp,
"\nL_Dnaa Version %d\n", &version);
1632 return (
L_DNAA *)ERROR_PTR(
"not a l_dna file", __func__, NULL);
1634 return (
L_DNAA *)ERROR_PTR(
"invalid l_dnaa version", __func__, NULL);
1635 if (fscanf(fp,
"Number of L_Dna = %d\n\n", &n) != 1)
1636 return (
L_DNAA *)ERROR_PTR(
"invalid number of l_dna", __func__, NULL);
1638 return (
L_DNAA *)ERROR_PTR(
"num l_dna <= 0", __func__, NULL);
1639 if (n > MaxPtrArraySize)
1640 return (
L_DNAA *)ERROR_PTR(
"too many l_dna", __func__, NULL);
1641 if (n == 0) L_INFO(
"the dnaa is empty\n", __func__);
1644 return (
L_DNAA *)ERROR_PTR(
"daa not made", __func__, NULL);
1645 for (i = 0; i < n; i++) {
1646 if (fscanf(fp,
"L_Dna[%d]:", &index) != 1) {
1648 return (
L_DNAA *)ERROR_PTR(
"invalid l_dna header", __func__, NULL);
1652 return (
L_DNAA *)ERROR_PTR(
"da not made", __func__, NULL);
1676 return (
L_DNAA *)ERROR_PTR(
"data not defined", __func__, NULL);
1678 return (
L_DNAA *)ERROR_PTR(
"stream not opened", __func__, NULL);
1682 if (!daa) L_ERROR(
"daa not read\n", __func__);
1702 return ERROR_INT(
"filename not defined", __func__, 1);
1704 return ERROR_INT(
"daa not defined", __func__, 1);
1707 return ERROR_INT(
"stream not opened", __func__, 1);
1711 return ERROR_INT(
"daa not written to stream", __func__, 1);
1731 return ERROR_INT(
"stream not defined", __func__, 1);
1733 return ERROR_INT(
"daa not defined", __func__, 1);
1737 fprintf(fp,
"Number of L_Dna = %d\n\n", n);
1738 for (i = 0; i < n; i++) {
1740 return ERROR_INT(
"da not found", __func__, 1);
1741 fprintf(fp,
"L_Dna[%d]:", i);
1771 if (pdata) *pdata = NULL;
1772 if (psize) *psize = 0;
1774 return ERROR_INT(
"&data not defined", __func__, 1);
1776 return ERROR_INT(
"&size not defined", __func__, 1);
1778 return ERROR_INT(
"daa not defined", __func__, 1);
1781 if ((fp = open_memstream((
char **)pdata, psize)) == NULL)
1782 return ERROR_INT(
"stream not opened", __func__, 1);
1786 *psize = *psize - 1;
1788 L_INFO(
"work-around: writing to a temp file\n", __func__);
1791 return ERROR_INT(
"tmpfile stream not opened", __func__, 1);
1793 if ((fp = tmpfile()) == NULL)
1794 return ERROR_INT(
"tmpfile stream not opened", __func__, 1);
#define DNA_VERSION_NUMBER
L_DNAA * l_dnaaRead(const char *filename)
l_dnaaRead()
L_DNA * l_dnaReadStream(FILE *fp)
l_dnaReadStream()
l_ok l_dnaWriteStream(FILE *fp, L_DNA *da)
l_dnaWriteStream()
l_ok l_dnaaReplaceDna(L_DNAA *daa, l_int32 index, L_DNA *da)
l_dnaaReplaceDna()
L_DNA * l_dnaCreate(l_int32 n)
l_dnaCreate()
l_int32 l_dnaaGetCount(L_DNAA *daa)
l_dnaaGetCount()
l_int32 l_dnaaGetDnaCount(L_DNAA *daa, l_int32 index)
l_dnaaGetDnaCount()
L_DNAA * l_dnaaReadMem(const l_uint8 *data, size_t size)
l_dnaaReadMem()
L_DNAA * l_dnaaCreateFull(l_int32 nptr, l_int32 n)
l_dnaaCreateFull()
l_ok l_dnaGetIValue(L_DNA *da, l_int32 index, l_int32 *pival)
l_dnaGetIValue()
L_DNAA * l_dnaaReadStream(FILE *fp)
l_dnaaReadStream()
l_ok l_dnaAddNumber(L_DNA *da, l_float64 val)
l_dnaAddNumber()
l_ok l_dnaCopyParameters(L_DNA *dad, L_DNA *das)
l_dnaCopyParameters()
l_ok l_dnaaTruncate(L_DNAA *daa)
l_dnaaTruncate()
L_DNA * l_dnaRead(const char *filename)
l_dnaRead()
L_DNA * l_dnaMakeSequence(l_float64 startval, l_float64 increment, l_int32 size)
l_dnaMakeSequence()
l_ok l_dnaSetValue(L_DNA *da, l_int32 index, l_float64 val)
l_dnaSetValue()
static l_int32 l_dnaaExtendArray(L_DNAA *daa)
l_dnaaExtendArray()
L_DNA * l_dnaCreateFromIArray(l_int32 *iarray, l_int32 size)
l_dnaCreateFromIArray()
L_DNA * l_dnaaGetDna(L_DNAA *daa, l_int32 index, l_int32 accessflag)
l_dnaaGetDna()
l_ok l_dnaaAddNumber(L_DNAA *daa, l_int32 index, l_float64 val)
l_dnaaAddNumber()
static l_int32 l_dnaExtendArray(L_DNA *da)
l_dnaExtendArray()
l_ok l_dnaaWriteStream(FILE *fp, L_DNAA *daa)
l_dnaaWriteStream()
l_ok l_dnaShiftValue(L_DNA *da, l_int32 index, l_float64 diff)
l_dnaShiftValue()
l_ok l_dnaRemoveNumber(L_DNA *da, l_int32 index)
l_dnaRemoveNumber()
l_ok l_dnaEmpty(L_DNA *da)
l_dnaEmpty()
l_float64 * l_dnaGetDArray(L_DNA *da, l_int32 copyflag)
l_dnaGetDArray()
L_DNAA * l_dnaaCreate(l_int32 n)
l_dnaaCreate()
l_ok l_dnaSetParameters(L_DNA *da, l_float64 startx, l_float64 delx)
l_dnaSetParameters()
l_ok l_dnaaWrite(const char *filename, L_DNAA *daa)
l_dnaaWrite()
l_ok l_dnaGetDValue(L_DNA *da, l_int32 index, l_float64 *pval)
l_dnaGetDValue()
static const l_int32 InitialArraySize
l_ok l_dnaaWriteMem(l_uint8 **pdata, size_t *psize, L_DNAA *daa)
l_dnaaWriteMem()
L_DNA * l_dnaCreateFromDArray(l_float64 *darray, l_int32 size, l_int32 copyflag)
l_dnaCreateFromDArray()
l_ok l_dnaInsertNumber(L_DNA *da, l_int32 index, l_float64 val)
l_dnaInsertNumber()
l_int32 l_dnaaGetNumberCount(L_DNAA *daa)
l_dnaaGetNumberCount()
L_DNA * l_dnaClone(L_DNA *da)
l_dnaClone()
l_ok l_dnaSetCount(L_DNA *da, l_int32 newcount)
l_dnaSetCount()
l_ok l_dnaaAddDna(L_DNAA *daa, L_DNA *da, l_int32 copyflag)
l_dnaaAddDna()
l_ok l_dnaaGetValue(L_DNAA *daa, l_int32 i, l_int32 j, l_float64 *pval)
l_dnaaGetValue()
l_ok l_dnaWriteStderr(L_DNA *da)
l_dnaWriteStrderr()
l_ok l_dnaReplaceNumber(L_DNA *da, l_int32 index, l_float64 val)
l_dnaReplaceNumber()
void l_dnaDestroy(L_DNA **pda)
l_dnaDestroy()
l_ok l_dnaGetParameters(L_DNA *da, l_float64 *pstartx, l_float64 *pdelx)
l_dnaGetParameters()
l_ok l_dnaWriteMem(l_uint8 **pdata, size_t *psize, L_DNA *da)
l_dnaWriteMem()
l_ok l_dnaWrite(const char *filename, L_DNA *da)
l_dnaWrite()
l_int32 * l_dnaGetIArray(L_DNA *da)
l_dnaGetIArray()
L_DNA * l_dnaCopy(L_DNA *da)
l_dnaCopy()
l_int32 l_dnaGetCount(L_DNA *da)
l_dnaGetCount()
L_DNA * l_dnaReadMem(const l_uint8 *data, size_t size)
l_dnaReadMem()
void l_dnaaDestroy(L_DNAA **pdaa)
l_dnaaDestroy()
void lept_stderr(const char *fmt,...)
lept_stderr()
l_uint8 * l_binaryReadStream(FILE *fp, size_t *pnbytes)
l_binaryReadStream()
FILE * fopenReadFromMemory(const l_uint8 *data, size_t size)
fopenReadFromMemory()
FILE * fopenWriteStream(const char *filename, const char *modestring)
fopenWriteStream()
void * reallocNew(void **pindata, size_t oldsize, size_t newsize)
reallocNew()
FILE * fopenWriteWinTempfile(void)
fopenWriteWinTempfile()
FILE * fopenReadStream(const char *filename)
fopenReadStream()