--- a/ChangeLog
+++ b/ChangeLog
@@ -1,5 +1,9 @@
 2011-10-09  Bruno Haible  <bruno@clisp.org>
 
+	Tests for module 'copysignf'.
+	* modules/copysignf-tests: New file.
+	* tests/test-copysignf.c: New file.
+
 	New module 'copysignf'.
 	* lib/math.in.h (copysignf): New declaration.
 	* lib/copysignf.c: New file.

new file mode 100644
--- /dev/null
+++ b/modules/copysignf-tests
@@ -0,0 +1,14 @@
+Files:
+tests/test-copysignf.c
+tests/signature.h
+tests/minus-zero.h
+tests/macros.h
+
+Depends-on:
+
+configure.ac:
+
+Makefile.am:
+TESTS += test-copysignf
+check_PROGRAMS += test-copysignf
+test_copysignf_LDADD = $(LDADD) @COPYSIGNF_LIBM@

new file mode 100644
--- /dev/null
+++ b/tests/test-copysignf.c
@@ -0,0 +1,111 @@
+/* Test of copysignf() function.
+   Copyright (C) 2010-2011 Free Software Foundation, Inc.
+
+   This program is free software: you can redistribute it and/or modify
+   it under the terms of the GNU General Public License as published by
+   the Free Software Foundation; either version 3 of the License, or
+   (at your option) any later version.
+
+   This program is distributed in the hope that it will be useful,
+   but WITHOUT ANY WARRANTY; without even the implied warranty of
+   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+   GNU General Public License for more details.
+
+   You should have received a copy of the GNU General Public License
+   along with this program.  If not, see <http://www.gnu.org/licenses/>.  */
+
+/* Written by Bruno Haible <bruno@clisp.org>, 2010-2011.  */
+
+#include <config.h>
+
+#include <math.h>
+
+#include "signature.h"
+SIGNATURE_CHECK (copysignf, double, (double, double));
+
+#include "macros.h"
+#include "minus-zero.h"
+
+#include <string.h>
+
+volatile float x;
+volatile float y;
+float z;
+float zero = 0.0f;
+
+int
+main ()
+{
+  /* A particular value in the first quadrant.  */
+  x = 0.6f;
+  y = 0.8f;
+  z = copysignf (x, y);
+  ASSERT (z == 0.6f);
+
+  /* A particular value in the second quadrant.  */
+  x = -0.6f;
+  y = 0.8f;
+  z = copysignf (x, y);
+  ASSERT (z == 0.6f);
+
+  /* A particular value in the third quadrant.  */
+  x = -0.6f;
+  y = -0.8f;
+  z = copysignf (x, y);
+  ASSERT (z == -0.6f);
+
+  /* A particular value in the fourth quadrant.  */
+  x = 0.6f;
+  y = -0.8f;
+  z = copysignf (x, y);
+  ASSERT (z == -0.6f);
+
+  /* From signed zero.  */
+  x = 1.0f;
+  y = 0.0f;
+  z = copysignf (x, y);
+  ASSERT (z == 1.0f);
+
+  x = 1.0f;
+  y = minus_zerof;
+  z = copysignf (x, y);
+  /* Assume all gnulib targets support -0.0f, until proven otherwise.  */
+  ASSERT (z == -1.0f);
+
+  x = -1.0f;
+  y = 0.0f;
+  z = copysignf (x, y);
+  ASSERT (z == 1.0f);
+
+  x = -1.0f;
+  y = minus_zerof;
+  z = copysignf (x, y);
+  ASSERT (z == -1.0f);
+
+  /* To signed zero.  */
+  x = 0.0f;
+  y = 1.0f;
+  z = copysignf (x, y);
+  ASSERT (z == 0.0f);
+  ASSERT (memcmp (&z, &zero, sizeof z) == 0);
+
+  x = 0.0f;
+  y = -1.0f;
+  z = copysignf (x, y);
+  ASSERT (z == 0.0f);
+  ASSERT (memcmp (&z, &zero, sizeof z) != 0);
+
+  x = minus_zerof;
+  y = 1.0f;
+  z = copysignf (x, y);
+  ASSERT (z == 0.0f);
+  ASSERT (memcmp (&z, &zero, sizeof z) == 0);
+
+  x = minus_zerof;
+  y = -1.0f;
+  z = copysignf (x, y);
+  ASSERT (z == 0.0f);
+  ASSERT (memcmp (&z, &zero, sizeof z) != 0);
+
+  return 0;
+}