enbid-ananke 0.3.0__cp310-cp310-manylinux_2_39_x86_64.whl

enbid_ananke/Enbid-2.0/parameterfiles/myparameterfile1

@@ -0,0 +1,37 @@
+%  Input and Output
+InitCondFile     Examples/hernquist1_small/snapshot_ici
+%InitCondFile      Examples/uniform_6d_box_4/snapshot_ici
+
+%-------------------------------------------------------
+ICFormat                  1     % O)ASCII 1)Gadget 2)User defined
+SnapshotFileBase        _ph1
+
+%-------------------------------------------------------
+% Tree related options
+SpatialScale            1.0 %  x->x/SpatialScale and v->v
+PartBoundary            7   % Min particles in a node to do boundary correction 
+NodeSplittingCriterion  1   % 0)Alternate 1) Min Entropy 
+CubicCells              1   % use 1 in spherically symmetric systems
+MedianSplittingOn       0   %
+
+%--------------------------------------------------------
+% Smoothing options  AM=adaptive metric Ker=Kernel Sp=Spherical Pr=Product 
+% 0) None 1)FiEstAS 2)Ker Sp Normal 3)Ker Sp AM 4)KerPr Normal 5)KerPr AM
+TypeOfSmoothing      0
+DesNumNgb            2   % 2-10 for Fiestas and 25-100 for Kernel
+VolCorr              1    %  0) Disbale 1) Enable 
+
+%--------------------------------------------------------
+% Kernel smoothing  related options
+% 0) B-Spline 1)top hat 2)Bi_weight (1-x^2)^2 3)Epanechikov 4)CIC 5)TSC 
+TypeOfKernel           3 
+KernelBiasCorrection   1    % 0)none 1)shift central particle 
+AnisotropicKernel      0    % 0) Isotropic 1) Anisotropic 
+Anisotropy             0    % fix minimum c/a minor/major axis ratio 
+DesNumNgbA             128   % Neighbors for cal covar metric for Anisotropic Ker 
+%--------------------------------------------------------
+% other miscellaneous option
+TypeListOn        0
+PeriodicBoundaryOn 0
+%--------------------------------------------------------
+

enbid_ananke/Enbid-2.0/parameterfiles/myparameterfile2

@@ -0,0 +1,37 @@
+%  Input and Output
+InitCondFile     Examples/hernquist1_small/snapshot_ici
+%InitCondFile      Examples/uniform_6d_box_4/snapshot_ici
+
+%-------------------------------------------------------
+ICFormat                  1     % O)ASCII 1)Gadget 2)User defined
+SnapshotFileBase        _ph2
+
+%-------------------------------------------------------
+% Tree related options
+SpatialScale            1.0 %  x->x/SpatialScale and v->v
+PartBoundary            7   % Min particles in a node to do boundary correction 
+NodeSplittingCriterion  1   % 0)Alternate 1) Min Entropy 
+CubicCells              1   % use 1 in spherically symmetric systems
+MedianSplittingOn       0   %
+
+%--------------------------------------------------------
+% Smoothing options  AM=adaptive metric Ker=Kernel Sp=Spherical Pr=Product 
+% 0) None 1)FiEstAS 2)Ker Sp Normal 3)Ker Sp AM 4)KerPr Normal 5)KerPr AM
+TypeOfSmoothing      1
+DesNumNgb            2   % 2-10 for Fiestas and 25-100 for Kernel 
+VolCorr              1    %  0) Disbale 1) Enable 
+
+%--------------------------------------------------------
+% Kernel smoothing  related options
+% 0) B-Spline 1)top hat 2)Bi_weight (1-x^2)^2 3)Epanechikov 4)CIC 5)TSC 
+TypeOfKernel           3 
+KernelBiasCorrection   1    % 0)none 1)shift central particle 
+AnisotropicKernel      0    % 0) Isotropic 1) Anisotropic 
+Anisotropy             0    % fix minimum c/a minor/major axis ratio 
+DesNumNgbA             128   % Neighbors for cal covar metric for Anisotropic Ker 
+%--------------------------------------------------------
+% other miscellaneous option
+TypeListOn        0
+PeriodicBoundaryOn 0
+%--------------------------------------------------------
+

enbid_ananke/Enbid-2.0/parameterfiles/myparameterfile3

@@ -0,0 +1,37 @@
+%  Input and Output
+InitCondFile     Examples/hernquist1_small/snapshot_ici
+%InitCondFile      Examples/uniform_6d_box_4/snapshot_ici
+
+%-------------------------------------------------------
+ICFormat                  1     % O)ASCII 1)Gadget 2)User defined
+SnapshotFileBase        _ph3
+
+%-------------------------------------------------------
+% Tree related options
+SpatialScale            1.0 %  x->x/SpatialScale and v->v
+PartBoundary            7   % Min particles in a node to do boundary correction 
+NodeSplittingCriterion  1   % 0)Alternate 1) Min Entropy 
+CubicCells              1   % use 1 in spherically symmetric systems
+MedianSplittingOn       0   %
+
+%--------------------------------------------------------
+% Smoothing options  AM=adaptive metric Ker=Kernel Sp=Spherical Pr=Product 
+% 0) None 1)FiEstAS 2)Ker Sp Normal 3)Ker Sp AM 4)KerPr Normal 5)KerPr AM
+TypeOfSmoothing      3
+DesNumNgb            40   % 2-10 for Fiestas and 25-100 for Kernel
+VolCorr              1    %  0) Disbale 1) Enable 
+
+%--------------------------------------------------------
+% Kernel smoothing  related options
+% 0) B-Spline 1)top hat 2)Bi_weight (1-x^2)^2 3)Epanechikov 4)CIC 5)TSC 
+TypeOfKernel           3 
+KernelBiasCorrection   1    % 0)none 1)shift central particle 
+AnisotropicKernel      0    % 0) Isotropic 1) Anisotropic 
+Anisotropy             0    % fix minimum c/a minor/major axis ratio 
+DesNumNgbA             128   % Neighbors for cal covar metric for Anisotropic Ker 
+%--------------------------------------------------------
+% other miscellaneous option
+TypeListOn        0
+PeriodicBoundaryOn 0
+%--------------------------------------------------------
+

enbid_ananke/Enbid-2.0/parameterfiles/myparameterfile4

@@ -0,0 +1,37 @@
+%  Input and Output
+InitCondFile     Examples/hernquist1_small/snapshot_ici
+%InitCondFile      Examples/uniform_6d_box_4/snapshot_ici
+
+%-------------------------------------------------------
+ICFormat                  1     % O)ASCII 1)Gadget 2)User defined
+SnapshotFileBase        _ph4
+
+%-------------------------------------------------------
+% Tree related options
+SpatialScale            1.0 %  x->x/SpatialScale and v->v
+PartBoundary            7   % Min particles in a node to do boundary correction 
+NodeSplittingCriterion  1   % 0)Alternate 1) Min Entropy 
+CubicCells              1   % use 1 in spherically symmetric systems
+MedianSplittingOn       0   %
+
+%--------------------------------------------------------
+% Smoothing options  AM=adaptive metric Ker=Kernel Sp=Spherical Pr=Product 
+% 0) None 1)FiEstAS 2)Ker Sp Normal 3)Ker Sp AM 4)KerPr Normal 5)KerPr AM
+TypeOfSmoothing      3
+DesNumNgb            40   % 2-10 for Fiestas and 25-100 for Kernel
+VolCorr              1    %  0) Disbale 1) Enable 
+
+%--------------------------------------------------------
+% Kernel smoothing  related options
+% 0) B-Spline 1)top hat 2)Bi_weight (1-x^2)^2 3)Epanechikov 4)CIC 5)TSC 
+TypeOfKernel           3 
+KernelBiasCorrection   1    % 0)none 1)shift central particle 
+AnisotropicKernel      1    % 0) Isotropic 1) Anisotropic 
+Anisotropy             0    % fix minimum c/a minor/major axis ratio 
+DesNumNgbA             64   % Neighbors for cal covar metric for Anisotropic Ker 
+%--------------------------------------------------------
+% other miscellaneous option
+TypeListOn        0
+PeriodicBoundaryOn 0
+%--------------------------------------------------------
+

enbid_ananke/Enbid-2.0/periodic_lengths.txt

@@ -0,0 +1,5 @@
+0 0 0 0 0 0 
+
+
+
+

enbid_ananke/Enbid-2.0/src/Makefile

@@ -0,0 +1,43 @@
+#  There are a number of code options that need to be switched on/off
+#  on compile-time.
+#
+#  From the following list, please select the desired options:
+OPT1     =  -DWARN        # to enable or disable warnings
+#OPT2     =  -DDIM6        # Does calculations in 6 dimensional space
+OPT2    =  -DDIM3         # Does calculations in 3 dimensional space
+#OPT2    =  -DDIMO         # specify dimension in allvars.h (#define ND 4) 
+#OPT3    =  -DMEDIAN       # reduces memory, if SPH not needed
+#OPT4    =  -DPERIODIC     # Does calculations in 6 dimensional space
+
+
+OPTIMIZE =    -O3 -Wall 
+# Extra for intel machines
+#OPTIMIZE =    -O3 -Wall -ffast-math -msse2
+# Extra for MAC OS X
+#OPTIMIZE =    -O3 -Wall -fast  
+
+OPTIONS =  $(OPTIMIZE) $(OPT1) $(OPT2) $(OPT3) $(OPT4) $(OPT5) $(OPT6)\
+           $(OPT7) $(OPT8) $(OPT9)  $(OPT10)
+
+CXX       = g++  #   
+CXXFLAGS =   $(OPTIONS)
+LIBS   =  -lm 
+
+OBJS   = main.o  begrun.o   init.o   io.o allvars.o functions.o  \
+	allocate.o  read_ic.o  nr.o    tree.o tree_search.o  ngb_search.o density_nd.o
+
+INCL   = allvars.h proto.h  nr.h functions.h ngb_search.h Makefile 
+
+
+EXEC   = Enbid
+
+$(EXEC): $(OBJS) 
+	$(CXX) $(OBJS) $(LIBS) -o  $(EXEC)  
+	mv $(EXEC) ..
+
+$(OBJS): $(INCL) 
+
+.PHONY : clean
+clean:
+	rm -f $(OBJS) 
+

enbid_ananke/Enbid-2.0/src/allocate.cpp

@@ -0,0 +1,47 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <string.h>
+#include <math.h>
+#include "allvars.h"
+#include "proto.h"
+#include "functions.h"
+
+
+/* This routine allocates memory for 
+ * particle storage.
+ */
+void allocate_memory(void)
+{
+  
+  
+    int bytes=0,bytes_tot=0;
+  
+    if(All.MaxPart>0)
+    {
+	if(!(P_data=new particle_data[All.MaxPart]))
+	{
+	    printf("failed to allocate memory for `P_data' (%d bytes).\n",bytes);
+	    endrun(1);
+	}
+	bytes_tot+=All.MaxPart*sizeof(struct particle_data);
+    
+	P= P_data-1;   /* start with offset 1 */
+	Part= &P[1];
+	printf("Allocated %g MByte for particle storage.\n",bytes_tot/(1024.0*1024.0));
+    }
+  
+}
+
+
+
+/* This routine frees the memory for the particle storage,
+ */
+void free_memory(void)
+{
+  
+    if(All.MaxPart>0)
+    {
+	free(P_data);
+    }
+  
+}

enbid_ananke/Enbid-2.0/src/allvars.cpp

@@ -0,0 +1,62 @@
+#include "allvars.h"
+//int    BatchFlag;
+int    NumPart;
+
+/* variables for input/output ,  usually only used on process 0 */
+char   ParameterFile[100];
+
+/* tabulated smoothing kernel */
+double  Kernel[KERNEL_TABLE+2],
+    KernelDer[KERNEL_TABLE+2],
+    KernelDer2[KERNEL_TABLE+2],
+    KernelRad[KERNEL_TABLE+2];
+
+/* this struct contains mostly code parameters read from the parameter file */
+struct global_data_all_processes All;
+
+vector<int> trees;
+FILE *treefile_fpt;
+//real BoxHalf, Box;
+
+/* The following structure holds all the information that is
+ * stored for each particle.
+ */
+
+vector<int> npart,npartc;
+
+struct io_header_1 header1;
+
+//struct param_data Var;
+
+struct particle_data  *P,*P_data,*Part;
+
+
+/* Header for the GADGET file format. */
+
+double t00,t11,t22;
+double xbmax[ND],xbmin[ND],kbuf[ND]; /* boundary calculation variables */
+
+
+real searchx[ND][2], searchcenter[ND],metric[ND],metric1[ND],gmatrix[ND][ND];
+int *list_kdt,*list_kd;
+float *listr;
+
+struct NODE
+*nodes,*nodes_base,*nodesC;
+struct SNODE
+*snodes,*snodes_base;
+struct PNODE
+*pnodes,*pnodes_base;
+
+
+struct XNODE
+*xnodes,*xnodes_base;
+
+struct linklist  *pqHead;
+struct pqueue *pqx1;
+int pnext,pnew;
+
+float **mrho,*d,**ve;
+
+
+

enbid_ananke/Enbid-2.0/src/allvars.h

@@ -0,0 +1,263 @@
+/* This file declares all global variables. Further variables should be added here,
+   and declared as 'extern'. The actual existence of these variables is provided by
+   the file 'allvars.cpp'. To produce 'allvars.cpp' from 'allvars.h', do the following:
+
+   1.) Erase all #define's
+   2.) add #include "allvars.h"
+   3.) delete all keywords 'extern'
+   4.) delete all struct definitions enclosed in {...}, e.g.
+   "extern struct global_data_all_processes {....} All;"
+   becomes "struct global_data_all_processes All;"
+   5.) delete extern struct xyz_data {...}; completely
+*/
+
+#include <stdio.h>
+#include<iostream>
+#include<vector>
+using namespace std;
+
+#ifdef MEDIAN
+#define ROOT		1
+#define LOWER(i,node)	(i<<1)
+#define UPPER(i,node)	((i<<1)+1)
+#define PARENT(i,node)	(i>>1)
+#define SIBLING(i,node) 	((i&1)?i-1:i+1)
+#define SETNEXT(i,node)\
+{\
+	while (i&1) i=i>>1;\
+	++i;\
+	}
+#else
+#define ROOT	1
+#define LOWER(i,node)	(node[i].left)
+#define UPPER(i,node)	(node[i].left+1)
+#define PARENT(i,node)	(node[i].parent)
+#define SIBLING(i,node) 	((i&1)?i-1:i+1)
+#define SETNEXT(i,node)\
+{\
+	while (i&1) i=node[i].parent;\
+	++i;\
+	}
+#endif
+
+
+#ifdef DIM6
+#define ND 6  //const int ND=6;
+#endif
+#ifdef DIM3
+#define ND 3  //const int ND=3;
+#endif
+#ifdef DIMO
+#define ND 4 //const int ND=5;
+#endif
+
+#ifdef T3E
+typedef short int int4byte;   /* Note: int has 8 Bytes on the T3E ! */
+#else
+typedef int int4byte;
+#endif
+
+#define real double
+#define real1 float
+
+/* ... often used constants (cgs units) */
+#define  MAX_REAL_NUMBER  1e37
+#define  MIN_REAL_NUMBER  1e-37
+#define  THIRD            (1.0/3.0)
+#ifndef  PI
+#define  PI               3.14159265358979323846
+#endif
+#define  PI_INV           (1/PI)
+#define  LN2              0.69314718
+#define  KERNEL_TABLE 1000
+#define  MAX_NGB  1000000  /* 20000defines maximum length of neighbour list */
+
+//extern int    BatchFlag;
+extern int    NumPart;
+
+/* variables for input/output ,  usually only used on process 0 */
+extern  char   ParameterFile[100];
+
+/* tabulated smoothing kernel */
+extern double  Kernel[KERNEL_TABLE+2],
+    KernelDer[KERNEL_TABLE+2],
+    KernelDer2[KERNEL_TABLE+2],
+    KernelRad[KERNEL_TABLE+2];
+
+/* this struct contains mostly code parameters read from the parameter file */
+extern struct global_data_all_processes
+{
+    /* Code options */
+    int    TypeOfSmoothing;
+    int    KernelBiasCorrection;
+    int    AnisotropicKernel;
+    int    Dimensions;
+    int    VolCorr;
+    int    TypeOfKernel;
+    double SpatialScale;
+    int    PartBoundary;
+    int    NodeSplittingCriterion;
+    int    CubicCells;
+    double Anisotropy;
+    int    DesNumNgb;
+    int    NumBucket;
+    int    DesNumNgbA;
+    int    NumBucketA;
+    /* Cosmology */
+//    double BoxSize, BoxHalf;
+    int    MedianSplittingOn;
+    float   hs[ND],hsv;
+    double  MassTable[6];
+    /* File options */
+    int   ICFormat;
+    double OmegaBaryon;
+    double Omega0;
+    char    InitCondFile[100],SnapshotFileBase[100];
+    char   InputDir[100],
+	InitFileBase[100];
+    /* Some other global parameters */
+    int   MaxPart,TypeListOn;
+    int flag_swap;
+    int order_flag;
+    int PeriodicBoundaryOn;
+#ifdef PERIODIC
+    double boxh[ND];
+#endif
+
+} All;
+
+
+
+
+extern struct particle_data
+{
+    int  ID,Type;//,NumNgb;           /* unique particle identifier */
+    float     Pos[ND];       // particle position
+#ifdef DIM3
+    float     Vel[3];       // particle velocity
+#endif
+    float     Mass;         /* particle mass */
+    float     Density;
+} *P,*P_data,*Part;//,*P1;
+
+
+extern vector<int> npart,npartc; 
+
+/* Header for the standard file format. */
+extern struct io_header_1
+{
+    int4byte npart[6];
+    double   mass[6];
+    double   time;
+    double   redshift;
+    int4byte flag_sfr;
+    int4byte flag_feedback;
+    int4byte npartTotal[6];
+    int4byte flag_cooling;
+    int4byte num_files;
+    double   BoxSize;
+    double   Omega0;
+    double   OmegaLambda;
+    double   HubbleParam;
+    /* extra flags other than gadget*/
+    int4byte flag_id;          /* if IDs needed in output */
+    int4byte flag_dim;         /* no of dimensions */
+    /* signifies that output file contains estimated density */
+    int4byte flag_density;
+    char     fill[256- 6*4- 6*8- 2*8- 2*4- 6*4- 2*4 - 4*8-3*4];  /* fills to 256 Bytes */
+} header1;
+
+extern double t00,t11,t22;
+extern double xbmax[ND],xbmin[ND],kbuf[ND]; /* boundary calculation variables */
+
+
+/* extern struct param_data */
+/* { */
+/*     /\* Code options *\/ */
+/*     int    BeginFileNum,EndFileNum,FileNum; */
+/* } Var; */
+
+
+
+
+
+extern int *list_kdt,*list_kd;//,*father;
+extern float *listr;
+extern real searchx[ND][2], searchcenter[ND],metric[ND],metric1[ND],gmatrix[ND][ND];//,h[ND];
+extern vector<int> trees;
+extern FILE *treefile_fpt;
+//extern real BoxHalf, Box;
+// count = no of particles
+// k1=dim to be split  ,kcut= split location
+// for leaf node k1 is index to its pnode
+// left and right left= index(pointer) to left and right nodes
+
+#ifdef MEDIAN
+extern struct NODE
+{
+    int count,k1,lid;  // 2*20=40 bytes,
+    real kcut;
+    struct XNODE * bnd;
+} *nodes,*nodes_base,*nodesC;
+#else
+extern struct NODE
+{
+    real kcut;
+    int count,k1,lid,left,parent;  // 2*20=40 bytes,
+    struct XNODE * bnd;
+} *nodes,*nodes_base,*nodesC;
+#endif
+
+//contains additional details of node
+// lid: first id of particle
+// k0:  last dimension split
+// bleft: 0 or 1 checks if the surface is derived from boundary or not
+extern struct SNODE
+{
+    //  real Density,h[ND];//,ex,ev,e[ND];  2*20=40bytes
+    int  lid,k0;
+    bool bleft[ND],bright[ND];
+} *snodes,*snodes_base;
+
+// contains co-ordinates (surfaces) of leaf node
+extern struct PNODE
+{
+    real x[ND][2];              // 48 bytes
+} *pnodes,*pnodes_base;
+
+
+// contains co-ordinates (surfaces) of all nodes
+extern struct XNODE
+{
+    real x[ND][2];            // 2*48=96 bytes +12bytes =
+} *xnodes,*xnodes_base;
+
+
+
+extern struct linklist 
+{
+    real r;
+    int p;
+    float x[ND];
+} *pqHead;
+
+extern struct pqueue 
+{
+    struct linklist* pq;
+    bool operator<( const struct pqueue &a ) const
+    {
+	return pq->r < a.pq->r;
+    }
+    bool operator>( const struct pqueue &a ) const
+    {
+	return pq->r > a.pq->r;
+    }
+}*pqx1;
+
+
+extern float **mrho,*d,**ve;
+extern int pnext,pnew;
+
+
+
+