Lemina
A molecular dynamics package for network, granular material and point particles with a range of interaction potential.
 
Loading...
Searching...
No Matches
AccumVacf.c File Reference
#include <stdio.h>
#include "global.h"
+ Include dependency graph for AccumVacf.c:

Go to the source code of this file.

Functions

double Integrate (double *, int)
 
void PrintVacf ()
 
void ZeroVacf ()
 
void AccumVacf ()
 

Function Documentation

◆ AccumVacf()

void AccumVacf ( )

Definition at line 27 of file AccumVacf.c.

27 {
28 double fac;
29 int j, nb;
30 for(nb = 1 ; nb <= nBuffAcf ; nb ++){
31 if(indexAcf[nb] == nValAcf){
32 for(j = 1 ; j <= nValAcf; j ++){
33 viscAcfAv[j] += viscAcf[nb][j];
34 }
35 indexAcf[nb] = 0;
36 countAcfAv ++;
38 fac = 1./(kinEnergy*region[1]*region[2]*limitAcfAv);
40 PrintVacf();
41 ZeroVacf();
42} } } }
void PrintVacf()
Definition PrintVacf.c:25
double Integrate(double *, int)
Definition Integrate.c:25
void ZeroVacf()
Definition ZeroVacf.c:25
double * viscAcfAv
Definition global.h:99
double viscAcfInt
Definition global.h:99
double region[2+1]
double kinEnergy
Definition global.h:20
int nValAcf
int nBuffAcf
Definition global.h:100
int stepAcf
Definition global.h:100
double * indexAcf
double deltaT
Definition global.h:20
double ** viscAcf
Definition global.h:99
int limitAcfAv
Definition global.h:100
int countAcfAv
Definition global.h:100

References countAcfAv, deltaT, indexAcf, Integrate(), kinEnergy, limitAcfAv, nBuffAcf, nValAcf, PrintVacf(), region, stepAcf, viscAcf, viscAcfAv, viscAcfInt, and ZeroVacf().

Referenced by EvalVacf().

+ Here is the call graph for this function:
+ Here is the caller graph for this function:

◆ Integrate()

double Integrate ( double * f,
int nf )

Definition at line 25 of file Integrate.c.

25 {
26 double s;
27 int i;
28 s = 0.5*(f[1] + f[nf]);
29 for(i = 2 ; i <= nf - 1 ; i ++)
30 s += f[i];
31 return(s);
32}

Referenced by AccumVacf().

+ Here is the caller graph for this function:

◆ PrintVacf()

void PrintVacf ( )

Definition at line 25 of file PrintVacf.c.

25 {
26 double tVal;
27 int j;
28 fprintf(fpvisc,"viscosity acf\n");
29 for(j = 1 ; j <= nValAcf ; j ++){
30 tVal = (j-1)*stepAcf*deltaT;
31 fprintf(fpvisc, "%lf\t %lf\t %lf\n", tVal, viscAcfAv[j], viscAcfAv[j]/viscAcfAv[1]);
32 }
33 fprintf(fpvisc, "viscosity acf integral : %lf\n", viscAcfInt);
34}
FILE * fpvisc

References deltaT, fpvisc, nValAcf, stepAcf, viscAcfAv, and viscAcfInt.

Referenced by AccumVacf().

+ Here is the caller graph for this function:

◆ ZeroVacf()

void ZeroVacf ( )

Definition at line 25 of file ZeroVacf.c.

25 {
26 int j;
27 countAcfAv= 0 ;
28 for(j = 1 ; j <= nValAcf ; j ++)
29 viscAcfAv[j] = 0.;
30}

References countAcfAv, nValAcf, and viscAcfAv.

Referenced by AccumVacf().

+ Here is the caller graph for this function: