def verify_using_ga(g_a, g_b, g_c):
g_v = ga.duplicate(g_c)
ga.gemm(False,False,N,N,N,1,g_a,g_b,0,g_v)
c = ga.access(g_c)
v = ga.access(g_v)
if c is not None:
val = int(np.abs(np.sum(c-v))>0.0001)
else:
val = 0
val = ga.gop_add(val)
ga.destroy(g_v)
return val == 0
def verify_using_ga(g_a, g_b, g_c):
g_v = ga.duplicate(g_c)
ga.gemm(False,False,N,N,N,1,g_a,g_b,0,g_v)
c = ga.access(g_c)
v = ga.access(g_v)
if c is not None:
val = int(np.abs(np.sum(c-v))>0.0001)
else:
val = 0
val = ga.gop_add(val)
ga.destroy(g_v)
return val == 0
def verify(g_a, g_b, g_c):
g_chk = ga.duplicate(g_a, "array check")
if not g_chk: ga.error("duplicate failed")
ga.sync()
ga.gemm(False, False, TOTALELEMS, TOTALELEMS, TOTALELEMS, 1.0, g_a, g_b,
0.0, g_chk)
ga.sync()
ga.add(g_c, g_chk, g_chk, 1.0, -1.0)
rchk = ga.dot(g_chk, g_chk)
if not me:
print "Normed difference in matrices: %12.4f" % rchk
if not (-TOLERANCE < rchk < TOLERANCE):
ga.error("Matrix multiply verify failed")
else:
print "Matrix Multiply OK"
ga.destroy(g_chk)
def verify(g_a, g_b, g_c):
g_chk = ga.duplicate(g_a, "array check")
if not g_chk: ga.error("duplicate failed")
ga.sync()
ga.gemm(False, False, TOTALELEMS, TOTALELEMS, TOTALELEMS, 1.0, g_a, g_b,
0.0, g_chk);
ga.sync()
ga.add(g_c, g_chk, g_chk, 1.0, -1.0)
rchk = ga.dot(g_chk, g_chk)
if not me:
print "Normed difference in matrices: %12.4f" % rchk
if not (-TOLERANCE < rchk < TOLERANCE):
ga.error("Matrix multiply verify failed")
else:
print "Matrix Multiply OK"
ga.destroy(g_chk)