[ESPResSo-users] lbfluid GPU vs CPU

[Markus Gusenbauer]

Hi,

there is an issue with lbfluid comparing GPU and CPU. I have a 
rectangular channel with wall velocity inflow from the left. With CPU 
the parabolic flow looks fine, but with GPU there is a standing wave 
(see pictures attached). I have a minimal working example below.

What could be the problem? I tried version 
"ESPResSo-3.4-dev-172-g5e97dc1-dirty-git" and the latest master from today.

Markus




setmd time_step 1.0
setmd skin 0.2
thermostat off

set boxX 30
set boxY 10
set boxZ 10

setmd box_l $boxX $boxY $boxZ

lbfluid cpu grid 1. dens 0.008 visc 0.032 tau 1.0 friction 0.5

# fluid inflow
lbboundary wall normal 1 0 0 dist 1 type 501 velocity 2.5e-4 0. 0.

lbboundary wall normal 0 1 0 dist 1 type 501
lbboundary wall normal 0 -1 0 dist [expr -$boxY+1] type 501
lbboundary wall normal 0 0 1 dist 1 type 501
lbboundary wall normal 0 0 -1 dist [expr -$boxZ+1] type 501

############################################################################
# main iteration loop
set cycle 0
set endCycle 2000

while { $cycle<$endCycle } {
    puts "$cycle"
    if { $cycle == 1999 } { lbfluid print vtk velocity "fluid.vtk" }
    integrate 10
    incr cycle
}

cpu.jpg
Description: JPEG image

gpu.jpg
Description: JPEG image