latticerule.h

Go to the documentation of this file.

/*
 * Authors:
 *
 *   David Munger <mungerd@iro.umontreal.ca>        (2015)
 *   Pierre L'Ecuyer <lecuyer@iro.umontreal.ca>     (2015)
 *
 */

#pragma once
#ifndef CLQMC_LATTICERULE_H
#define CLQMC_LATTICERULE_H

#include <clQMC/clQMC.h>
#include <stdio.h>

struct clqmcLatticeRule_;

typedef struct clqmcLatticeRule_ clqmcLatticeRule;


struct clqmcLatticeRuleStream_;

typedef struct clqmcLatticeRuleStream_ clqmcLatticeRuleStream;


#ifdef __cplusplus
extern "C" {
#endif

cl_uint clqmcLatticeRuleNumPoints(const clqmcLatticeRule* lattice);

cl_uint clqmcLatticeRuleDimension(const clqmcLatticeRule* lattice);

#define clqmcLatticeRuleCreate _CLQMC_TAG_FPTYPE(clqmcLatticeRuleCreate)

clqmcLatticeRule* clqmcLatticeRuleCreate             (cl_uint numPoints, cl_uint dimension, const cl_int* genVec, size_t* objectSize, clqmcStatus* err);
clqmcLatticeRule* clqmcLatticeRuleCreate_clqmc_float (cl_uint numPoints, cl_uint dimension, const cl_int* genVec, size_t* objectSize, clqmcStatus* err);
clqmcLatticeRule* clqmcLatticeRuleCreate_clqmc_double(cl_uint numPoints, cl_uint dimension, const cl_int* genVec, size_t* objectSize, clqmcStatus* err);

clqmcLatticeRule* clqmcLatticeRuleCreateKorobov(cl_uint numPoints, cl_uint dimension, cl_int gen, size_t* objectSize, clqmcStatus* err);

clqmcStatus       clqmcLatticeRuleDestroy(clqmcLatticeRule* lattice);

clqmcStatus       clqmcLatticeRuleWriteInfo(const clqmcLatticeRule* lattice, FILE* file);

#define clqmcLatticeRuleCreateStream       _CLQMC_TAG_FPTYPE(clqmcLatticeRuleCreateStream)
#define clqmcLatticeRuleCreateOverStream   _CLQMC_TAG_FPTYPE(clqmcLatticeRuleCreateOverStream)
#define clqmcLatticeRuleNextCoordinate     _CLQMC_TAG_FPTYPE(clqmcLatticeRuleNextCoordinate)
#define clqmcLatticeRuleNextPoint          _CLQMC_TAG_FPTYPE(clqmcLatticeRuleNextPoint)

CLQMCAPI clqmcLatticeRuleStream* clqmcLatticeRuleCreateStream             (const clqmcLatticeRule* lattice, cl_uint partCount, cl_uint partIndex, const _CLQMC_FPTYPE* shift, clqmcStatus* err);
CLQMCAPI clqmcLatticeRuleStream* clqmcLatticeRuleCreateStream_clqmc_float (const clqmcLatticeRule* lattice, cl_uint partCount, cl_uint partIndex, const cl_float*      shift, clqmcStatus* err);
CLQMCAPI clqmcLatticeRuleStream* clqmcLatticeRuleCreateStream_clqmc_double(const clqmcLatticeRule* lattice, cl_uint partCount, cl_uint partIndex, const cl_double*     shift, clqmcStatus* err);

CLQMCAPI clqmcStatus clqmcLatticeRuleCreateOverStream             (clqmcLatticeRuleStream* stream, const clqmcLatticeRule* lattice, cl_uint partCount, cl_uint partIndex, const _CLQMC_FPTYPE* shift);
CLQMCAPI clqmcStatus clqmcLatticeRuleCreateOverStream_clqmc_float (clqmcLatticeRuleStream* stream, const clqmcLatticeRule* lattice, cl_uint partCount, cl_uint partIndex, const cl_float*      shift);
CLQMCAPI clqmcStatus clqmcLatticeRuleCreateOverStream_clqmc_double(clqmcLatticeRuleStream* stream, const clqmcLatticeRule* lattice, cl_uint partCount, cl_uint partIndex, const cl_double*     shift);

CLQMCAPI clqmcStatus clqmcLatticeRuleDestroyStream(clqmcLatticeRuleStream* stream);

CLQMCAPI _CLQMC_FPTYPE clqmcLatticeRuleNextCoordinate         (clqmcLatticeRuleStream* stream);
CLQMCAPI cl_float  clqmcLatticeRuleNextCoordinate_clqmc_float (clqmcLatticeRuleStream* stream);
CLQMCAPI cl_double clqmcLatticeRuleNextCoordinate_clqmc_double(clqmcLatticeRuleStream* stream);

CLQMCAPI cl_uint clqmcLatticeRuleNextPoint             (clqmcLatticeRuleStream* stream, clqmc_fptype* coords);
CLQMCAPI cl_uint clqmcLatticeRuleNextPoint_clqmc_float (clqmcLatticeRuleStream* stream, clqmc_fptype* coords);
CLQMCAPI cl_uint clqmcLatticeRuleNextPoint_clqmc_double(clqmcLatticeRuleStream* stream, clqmc_fptype* coords);

CLQMCAPI cl_uint clqmcLatticeRuleForwardToNextPoint(clqmcLatticeRuleStream* stream);

CLQMCAPI cl_uint clqmcLatticeRuleCurrentPointIndex(const clqmcLatticeRuleStream* stream);

CLQMCAPI cl_uint clqmcLatticeRuleCurrentCoordIndex(const clqmcLatticeRuleStream* stream);

#ifdef __cplusplus
}
#endif

#endif