diff --git a/Gui/mainwindow.cpp b/Gui/mainwindow.cpp
index 6f93923da3772c54540f27d7a4f6b1413542b028..9d4383f3391ceb29f8019f9fed864b775bd03323 100644
--- a/Gui/mainwindow.cpp
+++ b/Gui/mainwindow.cpp
@@ -14,6 +14,10 @@
 
 #include "mainwindow.h"
 
+#include <numeric>
+#include "quicksort.h"
+#include "BoundingSphere.h"
+
 // ThirdParty/logog
 #include "logog/include/logog.hpp"
 
@@ -1298,6 +1302,15 @@ void MainWindow::showDataExplorerSettingsDialog()
 
 void MainWindow::FEMTestStart()
 {
+	std::vector<GeoLib::Point*> *pnts = const_cast<std::vector<GeoLib::Point*>*>(this->_project.getGEOObjects()->getPointVec("testgeo"));
+	//GeoLib::BoundingSphere s(*pnts);
+	//std::vector<size_t> perm(4);
+	//std::iota(perm.begin(), perm.end(), 0);
+	//BaseLib::Quicksort<GeoLib::Point*>(*pnts, 0, pnts->size(), perm);
+	GeoLib::BoundingSphere s(*(pnts->at(0)), *(pnts->at(1)), *(pnts->at(2)));//, *(pnts->at(3)));
+	GeoLib::Point c = s.getCenter();
+	double r = s.getRadius();
+	std::cout << "Center: (" << c[0] << ", " << c[1] << ", " << c[2] << "), Radius: " << r << std::endl;
 }
 
 void MainWindow::showTrackingSettingsDialog()