tfix vtk files - sphere - GPU-based 3D discrete element method algorithm with o… | |
git clone git://src.adamsgaard.dk/sphere | |
Log | |
Files | |
Refs | |
LICENSE | |
--- | |
commit a46ed4ed9cb15cae01a5897578a35aeace71c035 | |
parent db378368459b79a351971445588510f7479c9455 | |
Author: Anders Damsgaard <[email protected]> | |
Date: Tue, 14 Oct 2014 13:49:48 +0200 | |
fix vtk files | |
Diffstat: | |
M python/sphere.py | 13 ++++++++----- | |
1 file changed, 8 insertions(+), 5 deletions(-) | |
--- | |
diff --git a/python/sphere.py b/python/sphere.py | |
t@@ -5091,19 +5091,19 @@ class sim: | |
# First iteration: Allocate arrays and find constant values | |
if (i == 0): | |
# Shear displacement | |
- self.xdisp = numpy.zeros(lastfile+1, dtype=numpy.float… | |
+ self.xdisp = numpy.zeros(lastfile+1, dtype=numpy.float6… | |
# Normal stress | |
- self.sigma_eff = numpy.zeros(lastfile+1, dtype=numpy.float… | |
+ self.sigma_eff= numpy.zeros(lastfile+1, dtype=numpy.float6… | |
# Normal stress | |
- self.sigma_def = numpy.zeros(lastfile+1, dtype=numpy.float… | |
+ self.sigma_def= numpy.zeros(lastfile+1, dtype=numpy.float6… | |
# Shear stress | |
- self.tau = numpy.zeros(lastfile+1, dtype=numpy.float… | |
+ self.tau = numpy.zeros(lastfile+1, dtype=numpy.float6… | |
# Upper wall position | |
- self.dilation = numpy.zeros(lastfile+1, dtype=numpy.float… | |
+ self.dilation = numpy.zeros(lastfile+1, dtype=numpy.float6… | |
# Upper wall position | |
self.tau_p = 0.0 # Peak shear stress | |
t@@ -5116,6 +5116,9 @@ class sim: | |
w_x0 = sb.w_x[0] # Original height | |
A = sb.L[0] * sb.L[1] # Upper surface area | |
+ if (i == 1): | |
+ w_x0 = sb.w_x[0] # Original height | |
+ | |
# Summation of shear stress contributions | |
for j in fixvel[0]: | |
if (sb.vel[j,0] > 0.0): |