Adapt to MultiFloats.jl v2: exp(::MultiFloat) has been removed

src/SparseIR.jl

@@ -13,14 +13,19 @@ export AugmentedBasis, TauConst, TauLinear, MatsubaraConst
 export TauSampling, MatsubaraSampling, evaluate, fit, evaluate!, fit!,
     MatsubaraSampling64F, MatsubaraSampling64B, TauSampling64, sampling_points
 
-using MultiFloats: Float64x2
+using MultiFloats: MultiFloats, Float64x2, _call_big
+# If we call MultiFloats.use_bigfloat_transcendentals() like MultiFloats
+# recommends, we get an error during precompilation:
+for name in (:exp, :sinh, :cosh)
+    eval(:(function Base.$name(x::Float64x2)
+        Float64x2(_call_big($name, x, precision(Float64x2) + 20))
+    end))
+end
 using LinearAlgebra: LinearAlgebra, cond, dot, svd, SVD, QRIteration, mul!
 using QuadGK: gauss, quadgk
 using Bessels: sphericalbesselj
 using PrecompileTools
 
-include("_multifloat_funcs.jl")
-
 include("_linalg.jl")
 include("_roots.jl")
 include("_specfuncs.jl")

test/_multifloat_funcs.jl

@@ -18,12 +18,6 @@ logrange(x1, x2, length) = (exp10(y) for y in range(log10(x1), log10(x2); length
         @test cosh(x) isa Float64x2
     end
 
-    @testset "allocations" begin
-        x = rand(xx)
-        @test iszero(@allocated sinh(x))
-        @test iszero(@allocated cosh(x))
-    end
-
     @testset "sinh(x) where x = $x" for x in xx
         @test sinh(x)≈sinh(big(x)) rtol=eps(Float64x2)
     end