ICODE-ADC/lammps
4
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

more code formatting consistency changes for loops and conditionals

Browse Source
This commit is contained in:
Axel Kohlmeyer
2020-12-30 09:27:47 -05:00
parent 28de7e351a
commit 51a20175c9
171 changed files with 744 additions and 747 deletions
Download Patch File Download Diff File Expand all files Collapse all files

14
src/USER-MISC/pair_local_density.cpp
View File

@ -594,7 +594,7 @@ void PairLocalDensity::interpolate_cbspl(int n, double delta,
fpa = 0.;
fpb = 0.;
for ( i = 0; i < n-1; i++ ) {
for (i = 0; i < n-1; i++) {
dl[i] = du[i] = delta;
}
@ -602,32 +602,32 @@ void PairLocalDensity::interpolate_cbspl(int n, double delta,
dd[n-1] = 2.0 * delta;
coeff_c[0] = ( 3.0 / delta ) * ( f[1] - f[0] ) - 3.0 * fpa;
coeff_c[n-1] = 3.0 * fpb - ( 3.0 / delta ) * ( f[n-1] - f[n-2] );
for ( i = 0; i < n-2; i++ ) {
for (i = 0; i < n-2; i++) {
dd[i+1] = 4.0 * delta;
coeff_c[i+1] = ( 3.0 / delta ) * ( f[i+2] - f[i+1] ) -
( 3.0 / delta ) * ( f[i+1] - f[i] );
}
// tridiagonal solver
for ( i = 0; i < n-1; i++ ) {
for (i = 0; i < n-1; i++) {
du[i] /= dd[i];
dd[i+1] -= dl[i]*du[i];
}
coeff_c[0] /= dd[0];
for ( i = 1; i < n; i++ )
for (i = 1; i < n; i++)
coeff_c[i] = ( coeff_c[i] - dl[i-1] * coeff_c[i-1] ) / dd[i];
for ( i = n-2; i >= 0; i-- )
for (i = n-2; i >= 0; i--)
coeff_c[i] -= coeff_c[i+1] * du[i];
for ( i = 0; i < n-1; i++ ) {
for (i = 0; i < n-1; i++) {
coeff_d[i] = ( coeff_c[i+1] - coeff_c[i] ) / ( 3.0 * delta );
coeff_b[i] = ( f[i+1] - f[i] ) / delta - delta * ( coeff_c[i+1] + 2.0*coeff_c[i] ) / 3.0;
}
// normalize
for ( i = 0; i < n-1; i++ ) {
for (i = 0; i < n-1; i++) {
coeff_b[i] = coeff_b[i] * delta ;
coeff_c[i] = coeff_c[i] * delta*delta ;
coeff_d[i] = coeff_d[i] * delta*delta*delta;