diff --git a/tests/queries/0_stateless/00087_math_functions.reference b/tests/queries/0_stateless/00087_math_functions.reference
index 3e709ddfe52..e02aac11faf 100644
--- a/tests/queries/0_stateless/00087_math_functions.reference
+++ b/tests/queries/0_stateless/00087_math_functions.reference
@@ -72,3 +72,43 @@
 1
 1
 1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
+1
diff --git a/tests/queries/0_stateless/00087_math_functions.sql b/tests/queries/0_stateless/00087_math_functions.sql
index 8de8a774ff8..8b4b84ead9b 100644
--- a/tests/queries/0_stateless/00087_math_functions.sql
+++ b/tests/queries/0_stateless/00087_math_functions.sql
@@ -56,6 +56,12 @@ select log2(2) = 1;
 select log2(4) = 2;
 select sum(abs(log2(exp2(x)) - x) < 1.0e-9) / count() from system.one array join range(1000) as x;
 
+select log1p(-1) = -inf;
+select log1p(0) = 0;
+select abs(log1p(exp(2) - 1) - 2) < 1e8;
+select abs(log1p(exp(3) - 1) - 3) < 1e8;
+select sum(abs(log1p(exp(x) - 1) - x) < 1e-8) / count() from system.one array join range(100) as x;
+
 select sin(0) = 0;
 select sin(pi() / 4) = 1 / sqrt(2);
 select sin(pi() / 2) = 1;
@@ -82,6 +88,48 @@ select acos(-1) = pi();
 select atan(0) = 0;
 select atan(1) = pi() / 4;
 
+select atan2(0, 1) = 0;
+select atan2(0, 2) = 0;
+select atan2(1, 0) = pi() / 2;
+select atan2(1, 1) = pi() / 4;
+select atan2(-1, -1) = -3 * pi() / 4);
+
+select hypot(0, 1) = 1;
+select hypot(1, 0) = 1;
+select hypot(1, 1) = sqrt(2);
+select hypot(-1, 1) = sqrt(2);
+select hypot(3, 4) = 5;
+
+select sinh(0) = 0;
+select sinh(1) = -sinh(-1);
+select abs(sinh(1) - 0.5 * (e() - exp(-1))) < 1e-9;
+select abs(sinh(2) - 0.5 * (exp(2) - exp(-2))) < 1e-9;
+select (sum(abs(sinh(x) - 0.5 * (exp(x) - exp(-x)))) < 1e-9) / count() from system.one array join range(1000) as x;
+
+select cosh(0) = 1;
+select cosh(1) = cosh(-1);
+select abs(cosh(1) - 0.5 * (e() + exp(-1))) < 1e-9;
+select abs(pow(cosh(1), 2) - pow(sinh(1), 2) - 1) < 1e-9;
+select (sum(abs(cosh(x) * cosh(x) - sinh(x) * sinh(x) - 1)) < 1e-9) / count() from system.one array join range(1000) as x;
+
+select asinh(0) = 0;
+select asinh(1) = -asinh(-1);
+select abs(asinh(1) - ln(1 + sqrt(2))) < 1e-9;
+select (asinh(sinh(1)) - 1) < 1e-9
+select sum(abs(asinh(sinh(x)) - x) < 1e-9) / count() from system.one array join range(1000) as x;
+
+select acosh(1) = 0;
+select acosh(2) = acosh(-2);
+select abs(acosh(2) - ln(2 + sqrt(3))) < 1e-9;
+select (acosh(cosh(2)) - 2) < 1e-9
+select sum(abs(acosh(cosh(x)) - x) < 1e-9) / count() from system.one array join range(1, 1001) as x;
+
+select atanh(0) = 0;
+select atanh(0.5) = -asinh(-0.5);
+select abs(atanh(0.9) - 0.5 * ln(19)) < 1e-9;
+select (atanh(tanh(1)) - 1) < 1e-9
+select sum(abs(asinh(sinh(x)) - x) < 1e-9) / count() from system.one array join range(1000) as x;
+
 select erf(0) = 0;
 select erf(-10) = -1;
 select erf(10) = 1;