tRemoved porosity forcing at iteration=20 - sphere - GPU-based 3D discrete element method algorithm with optional fluid coupling
 (HTM) git clone git://src.adamsgaard.dk/sphere
 (DIR) Log
 (DIR) Files
 (DIR) Refs
 (DIR) LICENSE
       ---
 (DIR) commit df3e897f18c146bcf635a58d25a144b4ffc154fa
 (DIR) parent b68e846905de544f6a0d61356ed4c1b4b0c7cad5
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Tue,  1 Apr 2014 13:39:48 +0200
       
       Removed porosity forcing at iteration=20
       
       Diffstat:
         M src/navierstokes.cuh                |      12 ++++--------
       
       1 file changed, 4 insertions(+), 8 deletions(-)
       ---
 (DIR) diff --git a/src/navierstokes.cuh b/src/navierstokes.cuh
       t@@ -1018,13 +1018,8 @@ __global__ void findPorositiesVelocitiesDiametersSpherical(
                    //dev_ns_phi[cellidx]  = 1.0;
                    //dev_ns_dphi[cellidx] = 0.0;
        
       -            if (x == nx/2 && y == ny/2 && z == nz/2 && iteration == 20) {
       -                dev_ns_phi[cellidx]  = 0.9;
       -                dev_ns_dphi[cellidx] = -0.1;
       -            } else {
       -                dev_ns_phi[cellidx]  = 1.0;
       -                dev_ns_dphi[cellidx] = 0.0;
       -            }
       +            dev_ns_phi[cellidx]  = 1.0;
       +            dev_ns_dphi[cellidx] = 0.0;
        
                    dev_ns_vp_avg[cellidx] = MAKE_FLOAT3(0.0, 0.0, 0.0);
                    dev_ns_d_avg[cellidx]  = 0.0;
       t@@ -2427,7 +2422,8 @@ __global__ void applyParticleInteractionForce(
                }
            }
        }
       - 
       +
       +// Find the effective pressure that the top wall should exert
        
        // Print final heads and free memory
        void DEM::endNSdev()