a = np.ones((10,1))
b = np.ones((10,1))
c = np.memmap('zeros.mat', dtype=np.float64, mode='w+', shape=a.dot(b.T).shape, order='C')
a.dot(b.T, out=c)
print(a.dtype == b.dtype == c.dtype)
print(np.dot(a, b.T).shape == c.shape)
print(c.flags['C_CONTIGUOUS'])