Callable tridiagonal coefficients produces catastrophic loss of performance

[glwagner]

Currently BatchedTridiagonalSolver supports Callable tridiagonal coefficients via

Oceananigans.jl/src/Solvers/batched_tridiagonal_solver.jl

Lines 58 to 59 in bcc34f0

	@inline get_coefficient(a::Base.Callable, i, j, k, grid, p, args...) = a(i, j, k, grid, p, args...)
	@inline get_coefficient(a::Base.Callable, i, j, k, grid, ::Nothing, args...) = a(i, j, k, grid, args...)

However, as noted on #3030 this can produce a catastrophic loss of performance. Therefore, this support should be discontinued.

cc @simone-silvestri

[navidcy]

Any idea if this affects GPU or CPU or both?

[glwagner]

It would seem just CPU