Fokker_Planck_DKE_loop

Header source

0001 #include "tensor.h"
0002 
0003 void integral_dke_jd(struct Tensor A, struct Tensor B,int d,struct Tensor *C);
0004 /*Performs contraction of Tensor A with Tensor (vector B) along dimension d
0005 (dimensions start at 1) and returns C*/
0006 
0007 void Tensor_Create(int ndims,int* sizedim,struct Tensor *A);
0008 /*Initializes and allocates memory for Tensor with data specified by ndims,
0009 sizedim*/
0010 
0011 void Scalar_Tensor_product(DOUBLE s,struct Tensor A,struct Tensor* B);
0012 /*Performs *B=s*A*/
0013 
0014 int Tensor_Sum(struct Tensor A,struct Tensor B,DOUBLE m, struct Tensor *C);
0015 /*Performs *C=A+m*B*/
0016 
0017 void Tensor_Transpose(struct Tensor *A);
0018 /*Transposes *A if it's a vector or a matrix*/
0019 
0020 void SparseMat_Tensorproduct(struct SparseMatrix A,struct Tensor B,struct Tensor
0021 *C);
0022 /*Performs *C=A*B with A a sparsematrix as written by
0023 writeBinarySparseMat_mex_c.mex and B,*C Tensors (vectors)*/
0024 
0025 void Tensor_Fill(struct Tensor *T, DOUBLE a);
0026 /*Fills Tensor with DOUBLE value a*/
0027 
0028 void IntTensor_Fill(struct IntTensor *T, int a);
0029 /*Fills IntTensor with int value a)*/
0030 
0031 void ReadBinary_DKE(char* file, struct Tensor *A);
0032 /*READS Binary file containing Tensor as written by
0033 writeBinaryTensor_mex_c.mex*/
0034 
0035 void ReadBinaryArray_DKE(char* file,struct Tensor ** Array);
0036 /*READS Binary file containing Cell Array*/
0037 
0038 void ReadBinaryIntArray_DKE(char* file,struct IntTensor **Array);
0039 
0040 void ReadBinarySparseMatrix_DKE(char* file,struct SparseMatrix *A);
0041 /*READS SparseMatrix as written by writeBinarySparseMat_mex_c.mex*/
0042 
0043 void Tensor_Save(char * file, struct Tensor A);
0044 /*Writes Tensor A to file*/
0045 
0046 void IntTensor_Save(char * file, struct IntTensor A);
0047 /*Writes IntTensor A to file*/
0048 
0049 int firstordercollop_dke_ft_c2(struct Tensor f1,
0050 struct Tensor betath_ref,
0051 struct Tensor masku,struct Tensor maskl, 
0052 struct Tensor dpn, 
0053 struct Tensor xpn, struct Tensor xpn2, struct Tensor xpn3, struct Tensor xsigma, struct Tensor xgamma,
0054 struct Tensor xz, struct Tensor xz2,
0055 struct Tensor pn, struct Tensor pn2, struct Tensor pn3,  
0056 struct Tensor sigma, struct Tensor gamma,
0057 struct Tensor z, struct Tensor z2, 
0058 struct Tensor xJ1, struct Tensor xJ2, struct Tensor xJ3, struct Tensor xJ4,
0059 struct Tensor J1, struct Tensor J2, struct Tensor J3, struct Tensor J4,struct
0060 Tensor xTe_norm,struct Tensor *I1,struct Tensor *I2);
0061 /*C transalation of function firstordercollop_dke*/
0062 
0063 DOUBLE absolute(DOUBLE a);
0064 /*Returns absolute value of a*/
0065 
0066 int product(struct Tensor A,struct Tensor B,struct Tensor* C);
0067 int division(struct Tensor A,struct Tensor B,struct Tensor* C);