JULIA语言与并行计算

JULIA语言与并行计张常hangou@张先轶traits.zhang@gmail.om中国科学院软件研究所并行实验认识Helloprintln(“HelloLibrary-Friendly：C/Fortran/…shared顺序结 A=B=C=A+println(“C=,C选择结 x=1;ifx<println("xislessthany")elseifx>yprintln("xisgreaterthany")println("xisequaltoy")循环结 whilei<=5print(i+=1)endfori=1:5print(i)end;println(";")foriin[1,4,0]print(i)end;println(";")forsin["foo","bar","baz"]✓88✓88✓✓✓✓8

a=[1//3,3//5,7//9,8//9,2//7,9]foriinaprintln("doubleai-",2i)ca=1+2im;cb=2println(ca+cb)s="Julia,aprettylanguage!"println(typePoint{T}p1=Point{Int32}(123,p2=Point{Int32}(456,println("p1+p2=($(p1.x+p2.x),$(p1.y+p2.y))")Julia强调运行效functionsum(array)s=0foriinarrays+=ic=println("sum:",sum(c))f=+;调用C/FORTRAN调用C/FORTRAN语言函数制作或购买共享的动态库文WindowsLinux调用过获得库句柄查找需要的函数（dlsym(libc,获得结果CC语言库的例子-C程//C语言源码文件名#includedeclspec(dllexport)intsum(inta,intb){returna+b;}gcc-shared-otestdll4julia.dllCC语言库的例子-Julia程//Julia语言源码文件名libc=dlopen("testdll4julia")println("sumofthetwonumberis",rs)mean(v[,region]std(v[,region]stdm(v,mvar(v[,region]varm(v,mmedian(vhist(v[,n])→e,hist(v,e)→e,histrange(v,n

midpoints(equantile(v,pquantile(vcov(v1[,v2]cor(v1[,v2]CPU之间的通信使用内部消息机两个基础原Remotereferences：anobjectrefertoanotherRemotecalls:returnaremoteJulia可以管Julia可以管理多个CPU和CPU中的多计算核addprocs

x-线路，取y-站点，取

四惠四惠靛厂新西靛厂新靛厂新西靛厂新靛厂新四惠四惠靛厂新四惠西z-发车，取西

1 Y 2 3 4

杨杨西西国防大西西国防大前温都水XX计计算站间走行时计计算站间走行时长的平均busMeanZ=fori=forj=fork=busMeanZ[i,j]+=busMeanZ=Julia并行Julia并行的实现步添加处理器（核心juliajuliaconstX=4000;constY=40;constconstNP=4;添加处理器（核心添加处理器（核心juliajulia远程调用（并行远程调用（并行juliajulia拼合结果（时长拼合结果（时长juliajulia拼合结果（平均值远程调用（并行拼合结果（平均值远程调用（并行CPU：AMDPhenom(羿龙)IIX4960T -4核内存4GB宇瞻DDR31333MHz 计算耗时（串行计算耗时（串行processornumber:=计算公交线路站间行走时长(subinY)耗时elapsedtime:11.666052409=计算路段平均行走时长(meaninZ)耗时elapsedtime:11.810904218计算耗时（并行增加核心耗时：2.701658239processornumber:数据划分耗时：0.340269018走行时间耗时：2.015043338平均时间耗时：2.4313143meanDArrayCUDAmeanDArrayCUDA类核划平合串--并优拼合结果（平均值mean()优拼合结果（平均值functionzm(a::Array,x::Int32,y::Int32,z::Int32)forfor#fork=1:zsum[i,j]+=b[i,j,k]

mean()函mean()函数优类核划平合串--并优--UserJulia编辑器

MemCachedWebUserSession

OpenBLAS：开源高性能BLASIntelMKL，AMDACML，IBM开发者KazushigeGoto2010基于GotoBLAS21.13BSDOpenBLAS简目标：成为全球最好的BLAS开源实BSD协议张先轶,王茜,Zaheer完成龙芯3A完成IntelSandyBridgeBLAS3AMDBulldozerS/DGEMM改增强MacOSX,Windows,FreeBSD上的编译、安修正各种bug:SEGFAULT,计算结果错等OpenBLASOpenBLAS性能结龙芯3ABLAS34线程，OpenBLAS超过GotoBLAS120%，ATLASZhangXianyi,WangQian,ZhangYunquan,Model-drivenLevel3BLASPerformanceOptimizationonLoongson3AProcessor,2012IEEE18thInternationalConferenceonParallelandDistributedSystems