tupdate pressure update debug output - 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 fe438fd43b79359a6b4f1877f366fb1b0e4f066d
 (DIR) parent 8bf65a86fbf64159696822fd4bce622ca670f7a6
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Sun,  9 Nov 2014 10:55:27 +0100
       
       update pressure update debug output
       
       Diffstat:
         M src/darcy.cuh                       |      19 +++++++++----------
       
       1 file changed, 9 insertions(+), 10 deletions(-)
       ---
 (DIR) diff --git a/src/darcy.cuh b/src/darcy.cuh
       t@@ -903,20 +903,19 @@ __global__ void updateDarcySolution(
                const Float res_norm = (p_new - p)/(p + 1.0e-16);
        
        #ifdef REPORT_FORCING_TERMS
       +            const Float dp_diff = (ndem*devC_dt)/(beta_f*phi*mu)
       +                *(k*laplace_p + dot(grad_k, grad_p));
       +            const Float dp_forc = -dphi/(beta_f*phi*(1.0 - phi));
                printf("\n%d,%d,%d updateDarcySolution\n"
       -                "dpdt        = %e\n"
       -                "p           = %e\n"
                        "p_new       = %e\n"
       -                "f           = %e\n"
       -                "f_transient = %e\n"
       -                "f_forcing   = %e\n"
       -                "f_diff      = %e\n"
       +                "p           = %e\n"
       +                "p_old       = %e\n"
       +                "dp_diff     = %e\n"
       +                "dp_forc     = %e\n"
                        "res_norm    = %e\n",
                        x,y,z,
       -                dpdt,
       -                p, p_new,
       -                f,
       -                f_transient, f_forcing, f_diff,
       +                p_new, p, p_old,
       +                dp_diff, dp_forc,
                        res_norm);
        #endif