From 61b42ec69321dbd03fdc36eee29da950d475b900 Mon Sep 17 00:00:00 2001
From: Michael Hoffmann <hoffmann@inf.ethz.ch>
Date: Mon, 24 Apr 2023 17:37:12 +0200
Subject: [PATCH] add test case to check issue

---
 .../test/Bounding_volumes/issue_2836.cpp      | 39 +++++++++++++++++++
 1 file changed, 39 insertions(+)
 create mode 100644 Bounding_volumes/test/Bounding_volumes/issue_2836.cpp

diff --git a/Bounding_volumes/test/Bounding_volumes/issue_2836.cpp b/Bounding_volumes/test/Bounding_volumes/issue_2836.cpp
new file mode 100644
index 00000000000..9d47f1b9b39
--- /dev/null
+++ b/Bounding_volumes/test/Bounding_volumes/issue_2836.cpp
@@ -0,0 +1,39 @@
+#include <CGAL/Simple_cartesian.h>
+#include <CGAL/point_generators_2.h>
+#include <CGAL/rectangular_p_center_2.h>
+#include <CGAL/IO/Ostream_iterator.h>
+#include <CGAL/algorithm.h>
+#include <iostream>
+#include <algorithm>
+#include <vector>
+
+typedef double                                      FT;
+
+typedef CGAL::Simple_cartesian<FT>                  Kernel;
+
+typedef Kernel::Point_2                             Point;
+typedef std::vector<Point>                          Cont;
+typedef CGAL::Random_points_in_square_2<Point>      Generator;
+typedef CGAL::Ostream_iterator<Point,std::ostream>  OIterator;
+
+int main()
+{
+  CGAL::get_default_random() = CGAL::Random(1518508913);
+  int n = 10;
+  int p = 3;
+  OIterator cout_ip(std::cout);
+  CGAL::IO::set_pretty_mode(std::cout);
+
+  Cont points;
+  std::copy_n(Generator(1), n, std::back_inserter(points));
+  std::cout << "Generated Point Set:\n";
+  std::copy(points.begin(), points.end(), cout_ip);
+
+  FT p_radius;
+  std::cout << "\n\n" << p << "-centers:\n";
+  CGAL::rectangular_p_center_2(
+    points.begin(), points.end(), cout_ip, p_radius, p);
+  std::cout << "\n\n" << p << "-radius = " << p_radius << std::endl;
+
+  return 0;
+}