tTemporarily replaced contactLinear with contactLinearViscous in shearmodel2 - 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 beba9aac6c98f3bd322dd6f6da14f91d96699318
(DIR) parent 82d698d30f1988547ebf70a50e0fb9a7b295d99f
(HTM) Author: Anders Damsgaard <adc@geo.au.dk>
Date: Tue, 4 Sep 2012 09:22:51 +0200
Temporarily replaced contactLinear with contactLinearViscous in shearmodel2
Diffstat:
M src/contactsearch.cuh | 15 ++++++---------
1 file changed, 6 insertions(+), 9 deletions(-)
---
(DIR) diff --git a/src/contactsearch.cuh b/src/contactsearch.cuh
t@@ -411,9 +411,9 @@ __global__ void interact(unsigned int* dev_gridParticleIndex, // Input: Unsorted
Float delta_n, x_ab_length, radius_b;
Float3 x_ab;
Float4 x_b, distmod;
- Float4 vel_a = dev_vel_sorted[idx_a];
- Float4 angvel4_a = dev_angvel_sorted[idx_a];
- Float3 angvel_a = MAKE_FLOAT3(angvel4_a.x, angvel4_a.y, angvel4_a.z);
+ //Float4 vel_a = dev_vel_sorted[idx_a];
+ //Float4 angvel4_a = dev_angvel_sorted[idx_a];
+ //Float3 angvel_a = MAKE_FLOAT3(angvel4_a.x, angvel4_a.y, angvel4_a.z);
// Loop over all possible contacts, and remove contacts
// that no longer are valid (delta_n > 0.0)
t@@ -439,17 +439,14 @@ __global__ void interact(unsigned int* dev_gridParticleIndex, // Input: Unsorted
// Process collision if the particles are overlapping
if (delta_n < 0.0f) {
- //cuPrintf("\nProcessing contact, idx_a_orig = %u, idx_b_orig = %u, contact = %d, delta_n = %f\n",
- // idx_a_orig, idx_b_orig, i, delta_n);
- /*contactLinearViscous(&F, &T, &es_dot, &p,
+ contactLinearViscous(&F, &T, &es_dot, &p,
idx_a_orig, idx_b_orig,
dev_vel,
dev_angvel,
radius_a, radius_b,
x_ab, x_ab_length,
delta_n, devC_kappa);
- dev_delta_t[mempos] = MAKE_FLOAT4(0.0f, 0.0f, 0.0f, 0.0f);*/
- contactLinear(&F, &T, &es_dot, &p,
+ /*contactLinear(&F, &T, &es_dot, &p,
idx_a_orig,
idx_b_orig,
vel_a,
t@@ -459,7 +456,7 @@ __global__ void interact(unsigned int* dev_gridParticleIndex, // Input: Unsorted
radius_a, radius_b,
x_ab, x_ab_length,
delta_n, dev_delta_t,
- mempos);
+ mempos);*/
} else {
__syncthreads();
// Remove this contact (there is no particle with index=np)