diff --git a/Project.toml b/Project.toml
index 07f7695..5495ee3 100644
--- a/Project.toml
+++ b/Project.toml
@@ -1,7 +1,7 @@
 name = "NaturalNeighbours"
 uuid = "f16ad982-4edb-46b1-8125-78e5a8b5a9e6"
 authors = ["Daniel VandenHeuvel <danj.vandenheuvel@gmail.com>"]
-version = "1.2.0"
+version = "1.2.1"
 
 [deps]
 ChunkSplitters = "ae650224-84b6-46f8-82ea-d812ca08434e"
diff --git a/src/interpolation/eval.jl b/src/interpolation/eval.jl
index b2ab538..4b4d936 100644
--- a/src/interpolation/eval.jl
+++ b/src/interpolation/eval.jl
@@ -47,7 +47,7 @@ end
 @inline function _eval_interp(method::Union{<:Farin,Sibson{1}}, itp::NaturalNeighboursInterpolant, p, cache; kwargs...)
     gradients = get_gradient(itp)
     if isnothing(gradients)
-        throw(ArgumentError("Gradients must be provided for Sibson-1 or Farin interpolation. Consider using e.g. interpolate(tri, z; derivatives = true)."))
+        throw(ArgumentError("Gradients must be provided for Sibson-1, Farin, or Hiyoshi-2 interpolation. Consider using e.g. interpolate(tri, z; derivatives = true)."))
     end
     tri = get_triangulation(itp)
     nc = compute_natural_coordinates(Sibson(), tri, p, cache; kwargs...)
diff --git a/test/interpolation/basic_tests.jl b/test/interpolation/basic_tests.jl
index ffe7c84..6d85c86 100644
--- a/test/interpolation/basic_tests.jl
+++ b/test/interpolation/basic_tests.jl
@@ -73,7 +73,7 @@ end
     f′′ = (x, y) -> [2+6x*y 3x^2; 3x^2 2]
     z = [f(x, y) for (x, y) in each_point(tri)]
     itp = interpolate(tri, z; derivatives=false)
-    @test_throws ArgumentError("Gradients must be provided for Sibson-1 or Farin interpolation. Consider using e.g. interpolate(tri, z; derivatives = true).") itp(0.5, 0.5; method=Sibson(1))
+    @test_throws ArgumentError("Gradients must be provided for Sibson-1, Farin, or Hiyoshi-2 interpolation. Consider using e.g. interpolate(tri, z; derivatives = true).") itp(0.5, 0.5; method=Sibson(1))
 end
 
 @testset "Hiyoshi(2) errors without gradients and Hessians" begin