fixed aliasing of matrices in mul! call.

Author: thomvet

test/gpu_all.jl

@@ -35,15 +35,15 @@ end
 chunk_size = 10
 u0 = cu(rand(10,10)) #example kept small for test purposes.
 A  = cu(-randn(10,10))                  
-cache = dualcache(cu(zeros(10,10)), chunk_size)
-prob = ODEProblem(foo, u0, (0.0f0,1.0f0), (A, cache))
+cache = dualcache(cu(zeros(10, 10)), chunk_size)
+prob = ODEProblem(foo, u0, (0.0f0, 1.0f0), (A, cache))
 sol = solve(prob, TRBDF2(chunk_size = chunk_size))
 @test sol.retcode == :Success
 
 #with auto-detected chunk_size
 u0 = cu(rand(10,10)) #example kept small for test purposes.
 A  = cu(-randn(10,10))                  
-cache = dualcache(A)
+cache = dualcache(cu(zeros(10, 10)))
 prob = ODEProblem(foo, u0, (0.0f0,1.0f0), (A, cache))
 sol = solve(prob, TRBDF2())
 @test sol.retcode == :Success