mccc_gc_milstein.c File Reference

Milstein integrator for collision operator in GC picture. More...

#include <math.h>
#include <float.h>
#include "../../ascot5.h"
#include "../../consts.h"
#include "../../physlib.h"
#include "../../error.h"
#include "../../particle.h"
#include "../../B_field.h"
#include "../../plasma.h"
#include "../../random.h"
#include "mccc_wiener.h"
#include "mccc_coefs.h"
#include "mccc.h"

Go to the source code of this file.

Functions
void	mccc_gc_milstein (particle_simd_gc *p, real *hin, real *hout, real tol, mccc_wienarr *w, B_field_data *Bdata, plasma_data *pdata, mccc_data *mdata, real *rnd)
	Integrate collisions for one time-step.

Detailed Description

Milstein integrator for collision operator in GC picture.

Definition in file mccc_gc_milstein.c.

Function Documentation

mccc_gc_milstein()

void mccc_gc_milstein	(	particle_simd_gc *	p,
		real *	hin,
		real *	hout,
		real	tol,
		mccc_wienarr *	w,
		B_field_data *	Bdata,
		plasma_data *	pdata,
		mccc_data *	mdata,
		real *	rnd )

Integrate collisions for one time-step.

Parameters

p	pointer to gc simd struct
hin	time-steps for NSIMD markers
hout	suggestions for the next timesteps for NSIMD markers
tol	relative error tolerance
w	array holding wiener structs for NSIMD markers
Bdata	pointer to magnetic field data
pdata	pointer to plasma data
mdata	pointer to collision data struct
rnd	array of normally distributed random numbers used to resolve collisions. Values for marker i are rnd[i*NSIMD + j]

Definition at line 34 of file mccc_gc_milstein.c.