GPU-SPH源代码

这是三维的SPH方法模拟代码，具有半隐式开放边界哦！绝对的经典，而且是市面上难得一见的粒子法开放边界代码！

#!/usr/bin/env python

“““
This script converts an ASCII GRID Digital Elevation Model into a VTK Structured Points
dataset （legacy format） simply by replacing the ASCII GRID header with the appropriate VTK
metadata. 

The VTK dataset is two-dimensional （flat） so a Warp By Scalar filter should be applied in
ParaView to convert the Height field into actual height information.
“““

import sys
import string
import os

def usage（）:
    print “%s filename -- convert filename （ASCII GRID） into a VTK dataset“ % os.path.basename（sys.argv[0]）

if len（sys.argv） < 2:
    usage（）
    sys.exit（0）

grid_fname = sys.argv[1]
vts_fname = os.path.splitext（sys.argv[1]）[0] + ‘.vtk‘

grid = open（grid_fname ‘r‘）
vts = None


md = {}

# ASCII GRID stores lines north to south but in the VTK we want them south

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2017-07-05 14:29  gpusph-4.1\
     文件         577  2017-07-05 14:29  gpusph-4.1\.gitignore
     文件      106271  2017-07-05 14:29  gpusph-4.1\Doxygen_settings
     文件       35147  2017-07-05 14:29  gpusph-4.1\GPLv3.txt
     文件       43369  2017-07-05 14:29  gpusph-4.1\Makefile
     目录           0  2017-07-05 14:29  gpusph-4.1\asy\
     文件        1451  2017-07-05 14:29  gpusph-4.1\asy\euler_angles1.tex
     文件        2315  2017-07-05 14:29  gpusph-4.1\asy\euler_angles2.tex
     目录           0  2017-07-05 14:29  gpusph-4.1\aux\
     文件         838  2017-07-05 14:29  gpusph-4.1\aux\Codingstyle.txt
     文件       17476  2017-07-05 14:29  gpusph-4.1\aux\GPUSPH.xml
     文件       11712  2017-07-05 14:29  gpusph-4.1\aux\UsingGit.txt
     目录           0  2017-07-05 14:29  gpusph-4.1\docs\
     目录           0  2017-07-05 14:29  gpusph-4.1\docs\dev-guide\
     目录           0  2017-07-05 14:29  gpusph-4.1\docs\dev-guide\src\
     文件        1382  2017-07-05 14:29  gpusph-4.1\docs\dev-guide\src\group.dox
     文件         560  2017-07-05 14:29  gpusph-4.1\docs\dev-guide\src\main.dox
     目录           0  2017-07-05 14:29  gpusph-4.1\docs\user-guide\
     文件        2839  2017-07-05 14:29  gpusph-4.1\docs\user-guide\GPUSPHtemplate.cls
     文件         491  2017-07-05 14:29  gpusph-4.1\docs\user-guide\Makefile
     目录           0  2017-07-05 14:29  gpusph-4.1\docs\user-guide\fig\
     文件       32730  2017-07-05 14:29  gpusph-4.1\docs\user-guide\fig\vertex_and_boundary_elements.svg
     目录           0  2017-07-05 14:29  gpusph-4.1\docs\user-guide\gpusph-install\
     文件        1491  2017-07-05 14:29  gpusph-4.1\docs\user-guide\gpusph-install\Makefile
     文件       33595  2017-07-05 14:29  gpusph-4.1\docs\user-guide\gpusph-install\gpusph-install.tex
     文件       35248  2017-07-05 14:29  gpusph-4.1\docs\user-guide\gpusph-install\gpusph-license.tex
     文件       26172  2017-07-05 14:29  gpusph-4.1\docs\user-guide\gpusph-manual.bib
     目录           0  2017-07-05 14:29  gpusph-4.1\docs\user-guide\gpusph-setup\
     文件        1414  2017-07-05 14:29  gpusph-4.1\docs\user-guide\gpusph-setup\Makefile
     文件       46307  2017-07-05 14:29  gpusph-4.1\docs\user-guide\gpusph-setup\gpusph-setup.tex
     目录           0  2017-07-05 14:29  gpusph-4.1\docs\user-guide\gpusph-theory\
............此处省略214个文件信息