tRemove unused code - cngf-pf - continuum model for granular flows with pore-pr… | |
git clone git://src.adamsgaard.dk/cngf-pf | |
Log | |
Files | |
Refs | |
README | |
LICENSE | |
--- | |
commit 349c6b2a79c906087e03ff4462c3210e20653c3f | |
parent de519ac795f069865beefb13be55713b3413382e | |
Author: Anders Damsgaard <[email protected]> | |
Date: Mon, 16 Mar 2020 07:49:13 +0100 | |
Remove unused code | |
Diffstat: | |
M fluid.c | 32 -----------------------------… | |
1 file changed, 0 insertions(+), 32 deletions(-) | |
--- | |
diff --git a/fluid.c b/fluid.c | |
t@@ -136,38 +136,6 @@ darcy_pressure_change_1d(const int i, | |
return 1.0/(beta_f*phi[i]*mu_f)*div_k_grad_p; | |
} | |
-/*void | |
-darcy_pressure_change_1d(double* dp_f_dt, | |
- const int nz, | |
- const double* p_f_ghost_in, | |
- const double* phi, | |
- const double* k, | |
- const double dz, | |
- const double beta_f, | |
- const double mu_f) | |
-{ | |
- int i; | |
- double k_, div_k_grad_p, k_zn, k_zp; | |
- | |
-#pragma omp parallel for private(i, k_, div_k_grad_p, k_zn, k_zp) if(nz>12) | |
- for (i=0; i<nz; ++i) { | |
- | |
- k_ = k[i]; | |
- if (i==0) k_zn = k_; else k_zn = k[i-1]; | |
- if (i==nz-1) k_zp = k_; else k_zp = k[i+1]; | |
- | |
- div_k_grad_p = (2.0*k_zp*k_/(k_zp + k_) | |
- * (p_f_ghost_in[i+2] | |
- - p_f_ghost_in[i+1])/dz | |
- - 2.0*k_zn*k_/(k_zn + k_) | |
- * (p_f_ghost_in[i+1] | |
- - p_f_ghost_in[i])/dz | |
- )/dz; | |
- | |
- dp_f_dt[i] = 1.0/(beta_f*phi[i]*mu_f)*div_k_grad_p; | |
- } | |
-}*/ | |
- | |
int | |
darcy_solver_1d(const int max_iter, | |
const double rel_tol) |