volume_macros.h File Reference

#include "plmbase_config.h"
#include "plm_int.h"

Go to the source code of this file.

Macros
#define	COORDS_FROM_INDEX(ijk, idx, dim)
#define	LOOP_Z(ijk, fxyz, vol)
#define	LOOP_Z_OMP(k, vol)
#define	LOOP_Y(ijk, fxyz, vol)
#define	LOOP_X(ijk, fxyz, vol)
#define	PROJECT_Z(xyz, proj)    (xyz[0] * proj[2*3+0] + xyz[1] * proj[2*3+1] + xyz[2] * proj[2*3+2])
#define	PROJECT_Y(xyz, proj)    (xyz[0] * proj[1*3+0] + xyz[1] * proj[1*3+1] + xyz[2] * proj[1*3+2])
#define	PROJECT_X(xyz, proj)    (xyz[0] * proj[0*3+0] + xyz[1] * proj[0*3+1] + xyz[2] * proj[0*3+2])
#define	POSITION_FROM_COORDS(xyz, ijk, origin, step)

Functions
static plm_long	volume_index (const plm_long dims[3], plm_long i, plm_long j, plm_long k)
static plm_long	volume_index (const plm_long dims[3], const plm_long ijk[3])
static bool	index_in_volume (const plm_long dims[3], const plm_long ijk[3])

Macro Definition Documentation

COORDS_FROM_INDEX

#define COORDS_FROM_INDEX	(		ijk,
			idx,
			dim
	)

Value:

    ijk[2] = idx / (dim[0] * dim[1]);                                   \
    ijk[1] = (idx - (ijk[2] * dim[0] * dim[1])) / dim[0];               \
    ijk[0] = idx - ijk[2] * dim[0] * dim[1] - (ijk[1] * dim[0]);

LOOP_X

#define LOOP_X	(		ijk,
			fxyz,
			vol
	)

Value:

    for (                                                               \
    ijk[0] = 0,                                                         \
        fxyz[0] = vol->origin[0] + ijk[2]*vol->step[0*3+2]              \
        + ijk[1]*vol->step[0*3+1];                                      \
    ijk[0] < vol->dim[0];                                               \
    ++ijk[0],                                                           \
        fxyz[0] += vol->step[0*3+0],                                    \
        fxyz[1] += vol->step[1*3+0],                                    \
        fxyz[2] += vol->step[2*3+0]                                     \
        )

LOOP_Y

#define LOOP_Y	(		ijk,
			fxyz,
			vol
	)

Value:

    for (                                                               \
    ijk[1] = 0,                                                         \
        fxyz[1] = vol->origin[1] + ijk[2]*vol->step[1*3+2];             \
    ijk[1] < vol->dim[1];                                               \
    ++ijk[1],                                                           \
        fxyz[2] = vol->origin[2] + ijk[2]*vol->step[2*3+2]              \
        + ijk[1]*vol->step[2*3+1],                                      \
        fxyz[1] = vol->origin[1] + ijk[2]*vol->step[1*3+2]              \
        + ijk[1]*vol->step[1*3+1]                                       \
        )

LOOP_Z

#define LOOP_Z	(		ijk,
			fxyz,
			vol
	)

Value:

    for (                                                               \
    ijk[2] = 0,                                                         \
        fxyz[2] = vol->origin[2];                                       \
    ijk[2] < vol->dim[2];                                               \
    ++ijk[2],                                                           \
        fxyz[2] = vol->origin[2] + ijk[2]*vol->step[2*3+2]              \
        )

LOOP_Z_OMP

#define LOOP_Z_OMP	(		k,
			vol
	)

Value:

    for (                                                               \
    long k = 0;                                                         \
    k < vol->dim[2];                                                    \
    ++k                                                                 \
        )

POSITION_FROM_COORDS

#define POSITION_FROM_COORDS	(		xyz,
			ijk,
			origin,
			step
	)

Value:

    do {                                                                \
        xyz[0] = origin[0]                                              \
            + ijk[0]*step[3*0+0]                                        \
            + ijk[1]*step[3*0+1]                                        \
            + ijk[2]*step[3*0+2];                                       \
        xyz[1] = origin[1]                                              \
            + ijk[0]*step[3*1+0]                                        \
            + ijk[1]*step[3*1+1]                                        \
            + ijk[2]*step[3*1+2];                                       \
        xyz[2] = origin[2]                                              \
            + ijk[0]*step[3*2+0]                                        \
            + ijk[1]*step[3*2+1]                                        \
            + ijk[2]*step[3*2+2];                                       \
    } while (0);

PROJECT_X

#define PROJECT_X	(		xyz,
			proj
	)		(xyz[0] * proj[0*3+0] + xyz[1] * proj[0*3+1] + xyz[2] * proj[0*3+2])

PROJECT_Y

#define PROJECT_Y	(		xyz,
			proj
	)		(xyz[0] * proj[1*3+0] + xyz[1] * proj[1*3+1] + xyz[2] * proj[1*3+2])

PROJECT_Z

#define PROJECT_Z	(		xyz,
			proj
	)		(xyz[0] * proj[2*3+0] + xyz[1] * proj[2*3+1] + xyz[2] * proj[2*3+2])

Function Documentation

index_in_volume()

static bool index_in_volume ( const plm_long  dims[3], const plm_long  ijk[3]  )

inlinestatic

volume_index() [1/2]

static plm_long volume_index ( const plm_long  dims[3], const plm_long  ijk[3]  )

inlinestatic

volume_index() [2/2]

static plm_long volume_index ( const plm_long  dims[3], plm_long  i, plm_long  j, plm_long  k  )

inlinestatic