@@ -295,10 +295,8 @@ function algebraicmultigrid2(W, du, u, p, t, newW, Plprev, Prprev, solverdata)
295295 if newW === nothing || newW
296296 A = convert(AbstractMatrix, W)
297297 Pl = AlgebraicMultigrid.aspreconditioner(AlgebraicMultigrid.ruge_stuben(A,
298- presmoother = AlgebraicMultigrid.Jacobi(rand(size(A,
299- 1))),
300- postsmoother = AlgebraicMultigrid.Jacobi(rand(size(A,
301- 1)))))
298+ presmoother = AlgebraicMultigrid.Jacobi(rand(size(A, 1))),
299+ postsmoother = AlgebraicMultigrid.Jacobi(rand(size(A, 1)))))
302300 else
303301 Pl = Plprev
304302 end
@@ -415,10 +413,8 @@ And similarly for algebraic multigrid:
415413
416414``` julia
417415prectmp2 = AlgebraicMultigrid. aspreconditioner (AlgebraicMultigrid. ruge_stuben (W,
418- presmoother = AlgebraicMultigrid. Jacobi (rand (size (W,
419- 1 ))),
420- postsmoother = AlgebraicMultigrid. Jacobi (rand (size (W,
421- 1 )))))
416+ presmoother = AlgebraicMultigrid. Jacobi (rand (size (W, 1 ))),
417+ postsmoother = AlgebraicMultigrid. Jacobi (rand (size (W, 1 )))))
422418const preccache2 = Ref (prectmp2)
423419function psetupamg (p, t, u, du, jok, jcurPtr, gamma)
424420 if jok
@@ -433,10 +429,8 @@ function psetupamg(p, t, u, du, jok, jcurPtr, gamma)
433429 # Build preconditioner on W
434430 preccache2[] = AlgebraicMultigrid. aspreconditioner (AlgebraicMultigrid. ruge_stuben (
435431 W,
436- presmoother = AlgebraicMultigrid. Jacobi (rand (size (W,
437- 1 ))),
438- postsmoother = AlgebraicMultigrid. Jacobi (rand (size (W,
439- 1 )))))
432+ presmoother = AlgebraicMultigrid. Jacobi (rand (size (W, 1 ))),
433+ postsmoother = AlgebraicMultigrid. Jacobi (rand (size (W, 1 )))))
440434 end
441435end
442436
0 commit comments