densitykernel.h File Reference

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Classes
struct	DensityKernel

Macros
#define	NUMDIMS   3
#define	NORM_COEFF   4.188790204786

Enumerations
enum	DensityKernelType { DENSITY_KERNEL_CUBIC_SPLINE = 1 , DENSITY_KERNEL_QUINTIC_SPLINE = 2 , DENSITY_KERNEL_QUARTIC_SPLINE = 4 }

Functions
double	density_kernel_desnumngb (DensityKernel *kernel, double eta)
void	density_kernel_init (DensityKernel *kernel, double H, enum DensityKernelType type)
double	density_kernel_wk (DensityKernel *kernel, double u)
double	density_kernel_dwk (DensityKernel *kernel, double u)
double	density_kernel_volume (DensityKernel *kernel)
static double	density_kernel_dW (DensityKernel *kernel, double u, double wk, double dwk)
static double	dotproduct (const double v1[3], const double v2[3])
static void	crossproduct (const double v1[3], const double v2[3], double out[3])

Macro Definition Documentation

NORM_COEFF

#define NORM_COEFF   4.188790204786

Coefficient for kernel normalization. Note: 4.0/3 * PI = 4.188790204786

Definition at line 6 of file densitykernel.h.

NUMDIMS

#define NUMDIMS   3

For 3D-normalized kernel

Definition at line 5 of file densitykernel.h.

Enumeration Type Documentation

DensityKernelType

enum DensityKernelType

Enumerator
DENSITY_KERNEL_CUBIC_SPLINE
DENSITY_KERNEL_QUINTIC_SPLINE
DENSITY_KERNEL_QUARTIC_SPLINE

Definition at line 17 of file densitykernel.h.

                       {
    DENSITY_KERNEL_CUBIC_SPLINE = 1,
    DENSITY_KERNEL_QUINTIC_SPLINE = 2,
    DENSITY_KERNEL_QUARTIC_SPLINE = 4,
};

Function Documentation

crossproduct()

static void crossproduct ( const double  v1[3], const double  v2[3], double  out[3]  )

inlinestatic

Definition at line 63 of file densitykernel.h.

{
    static const int D2[3] = {1, 2, 0};
    static const int D3[3] = {2, 0, 1};
 
    int d1;
    for(d1 = 0; d1 < 3; d1++)
    {
        const int d2 = D2[d1];
        const int d3 = D3[d1];
        out[d1] = (v1[d2] * v2[d3] -  v2[d2] * v1[d3]);
    }
}

Referenced by add_particle_to_group(), density_ngbiter(), and fof_finish_group_properties().

Here is the caller graph for this function:

density_kernel_desnumngb()

double density_kernel_desnumngb	(	DensityKernel *	kernel,
		double	eta
	)

Definition at line 124 of file densitykernel.c.

{
    /* this returns the expected number of ngb in for given sph resolution
     * deseta */
    /* See Price: arxiv 1012.1885. eq 12 */
    double support = kernel->support;
    return NORM_COEFF * pow(support * eta, NUMDIMS);
}

References kernel(), NORM_COEFF, NUMDIMS, and support.

Referenced by GetDensityKernelType().

Here is the call graph for this function:

Here is the caller graph for this function:

density_kernel_dW()

static double density_kernel_dW ( DensityKernel *  kernel, double  u, double  wk, double  dwk  )

inlinestatic

Definition at line 47 of file densitykernel.h.

{
    return - (NUMDIMS * kernel->Hinv * wk + u * dwk);
}

References dwk, kernel(), NUMDIMS, and wk.

Referenced by density_ngbiter().

Here is the call graph for this function:

Here is the caller graph for this function:

density_kernel_dwk()

double density_kernel_dwk	(	DensityKernel *	kernel,
		double	u
	)

Definition at line 108 of file densitykernel.c.

{
    double support = KERNELS[kernel->type].support;
    return kernel->dWknorm *
        KERNELS[kernel->type].dwk(kernel, u * support);
}

References kernel(), KERNELS, and support.

Referenced by density_ngbiter(), and hydro_ngbiter().

Here is the call graph for this function:

Here is the caller graph for this function:

density_kernel_init()

void density_kernel_init	(	DensityKernel *	kernel,
		double	H,
		enum DensityKernelType	type
	)

Definition at line 157 of file densitykernel.c.

{
    int t = -1;
    if(type == DENSITY_KERNEL_CUBIC_SPLINE) {
        t = 0;
    } else
    if(type == DENSITY_KERNEL_QUINTIC_SPLINE) {
        t = 1;
    } else
    if(type == DENSITY_KERNEL_QUARTIC_SPLINE) {
        t = 2;
    } else {
        endrun(1, "Density Kernel type is unknown\n");
    }
    density_kernel_init_with_type(kernel, t, H);
}

References DENSITY_KERNEL_CUBIC_SPLINE, density_kernel_init_with_type(), DENSITY_KERNEL_QUARTIC_SPLINE, DENSITY_KERNEL_QUINTIC_SPLINE, endrun(), and kernel().

Referenced by blackhole_accretion_ngbiter(), blackhole_dynfric_ngbiter(), blackhole_feedback_ngbiter(), density_ngbiter(), GetDensityKernelType(), hydro_ngbiter(), metal_return_ngbiter(), set_density_params(), and stellar_density_ngbiter().

Here is the call graph for this function:

Here is the caller graph for this function:

density_kernel_volume()

double density_kernel_volume

(

DensityKernel *

kernel

)

Definition at line 134 of file densitykernel.c.

{
    return NORM_COEFF * pow(kernel->H, NUMDIMS);
}

References kernel(), NORM_COEFF, and NUMDIMS.

Referenced by density_ngbiter(), and stellar_density_ngbiter().

Here is the call graph for this function:

Here is the caller graph for this function:

density_kernel_wk()

double density_kernel_wk	(	DensityKernel *	kernel,
		double	u
	)

Definition at line 116 of file densitykernel.c.

{
    double support = KERNELS[kernel->type].support;
    return kernel->Wknorm *
        KERNELS[kernel->type].wk(kernel, u * support);
}

References kernel(), KERNELS, and support.

Referenced by blackhole_accretion_ngbiter(), blackhole_dynfric_ngbiter(), blackhole_feedback_ngbiter(), density_ngbiter(), metal_return_ngbiter(), and stellar_density_ngbiter().

Here is the call graph for this function:

Here is the caller graph for this function:

dotproduct()

static double dotproduct ( const double  v1[3], const double  v2[3]  )

inlinestatic

Definition at line 53 of file densitykernel.h.

{
    double r =0;
    int d;
    for(d = 0; d < 3; d ++) {
        r += v1[d] * v2[d];
    }
    return r;
}

Referenced by density_ngbiter(), and hydro_ngbiter().

Here is the caller graph for this function: