mirror of
https://bitbucket.org/cosmicvoids/vide_public.git
synced 2025-07-04 07:11:12 +00:00
58 lines
1.6 KiB
C
58 lines
1.6 KiB
C
#include <stdio.h>
|
|
#include <stddef.h>
|
|
#include <stdlib.h>
|
|
#define NR_END 1
|
|
#define FREE_ARG char*
|
|
|
|
|
|
int ***i3tensor_2(long nrl, long nrh, long ncl, long nch, long ndl, long ndh)
|
|
/* allocate a float 3tensor with range t[nrl..nrh][ncl..nch][ndl..ndh] */
|
|
{
|
|
long i,j,nrow=nrh-nrl+1,ncol=nch-ncl+1,ndep=ndh-ndl+1;
|
|
int ***t;
|
|
|
|
/* allocate pointers to pointers to rows */
|
|
t=(int ***) malloc((size_t)((nrow+NR_END)*sizeof(int**)));
|
|
if (!t)exit(0);
|
|
t += NR_END;
|
|
t -= nrl;
|
|
|
|
/* allocate pointers to rows and set pointers to them */
|
|
t[nrl]=(int **) malloc((size_t)((nrow*ncol+NR_END)*sizeof(int*)));
|
|
if (!t[nrl])exit(0);
|
|
t[nrl] += NR_END;
|
|
t[nrl] -= ncl;
|
|
|
|
/* allocate rows and set pointers to them */
|
|
t[nrl][ncl]=(int *) malloc((size_t)((nrow*ncol*ndep+NR_END)*sizeof(int)));
|
|
if (!t[nrl][ncl])exit(0);
|
|
t[nrl][ncl] += NR_END;
|
|
t[nrl][ncl] -= ndl;
|
|
|
|
for(j=ncl+1;j<=nch;j++) t[nrl][j]=t[nrl][j-1]+ndep;
|
|
for(i=nrl+1;i<=nrh;i++) {
|
|
t[i]=t[i-1]+ncol;
|
|
t[i][ncl]=t[i-1][ncl]+ncol*ndep;
|
|
for(j=ncl+1;j<=nch;j++) t[i][j]=t[i][j-1]+ndep;
|
|
}
|
|
|
|
/* return pointer to array of pointers to rows */
|
|
return t;
|
|
}
|
|
|
|
|
|
|
|
void free_i3tensor(int ***t, long nrl, long nrh, long ncl, long nch,
|
|
long ndl, long ndh)
|
|
/* free a int i3tensor allocated by f3tensor() */
|
|
{
|
|
free((FREE_ARG) (t[nrl][ncl]+ndl-NR_END));
|
|
free((FREE_ARG) (t[nrl]+ncl-NR_END));
|
|
free((FREE_ARG) (t+nrl-NR_END));
|
|
}
|
|
|
|
|
|
|
|
|
|
#undef NR_END
|
|
#undef FREE_ARG
|