diff --git a/README.md b/README.md
index 6927ac1..6f900de 100644
--- a/README.md
+++ b/README.md
@@ -11,7 +11,7 @@ jacobi_pd
 ## Description
 
 This repository contains a small C++
-[header file](include/jacobi.hpp)
+[header file](include/jacobi_pd.hpp)
 that implements the
 [Jacobi eigenvalue algorithm](https://en.wikipedia.org/wiki/Jacobi_eigenvalue_algorithm).
 It is [free of copyright](https://creativecommons.org/publicdomain/zero/1.0/).
@@ -46,7 +46,7 @@ No lines of code were borrowed or adapted from other sources.
 *Caveats:* The code in this repository does not run in parallel,
 and only works on dense square real symmetric matrices.
 However it is reasonably
-[short, simple](include/jacobi.hpp), 
+[short, simple](include/jacobi_pd.hpp), 
 [fast](benchmarks/README.md) and
 [reliable](.travis.yml).
 You can do anything you like with this code.
@@ -55,8 +55,8 @@ You can do anything you like with this code.
 ##  Example usage
 
 ```cpp
-#include "jacobi.hpp"
-using namespace jacobi_public_domain;
+#include "jacobi_pd.hpp"
+using namespace jacobi_pd;
 
 // ...
 
diff --git a/include/jacobi_pd.hpp b/include/jacobi_pd.hpp
index 30faa00..9ceb2c8 100644
--- a/include/jacobi_pd.hpp
+++ b/include/jacobi_pd.hpp
@@ -7,11 +7,12 @@
 #ifndef _JACOBI_HPP
 #define _JACOBI_HPP
 
+#include <algorithm>
 #include <cmath>
 #include <cassert>
 #include "matrix_alloc.hpp"
 
-namespace jacobi_public_domain {
+namespace jacobi_pd {
 
 using namespace matrix_alloc;
 
diff --git a/tests/test_jacobi.cpp b/tests/test_jacobi.cpp
index 1712527..897a24c 100644
--- a/tests/test_jacobi.cpp
+++ b/tests/test_jacobi.cpp
@@ -11,7 +11,7 @@
 #include <vector>
 #include <array>
 #include "matrix_alloc.hpp"
-#include "jacobi.hpp"
+#include "jacobi_pd.hpp"
 
 using std::cout;
 using std::cerr;
@@ -20,7 +20,7 @@ using std::setprecision;
 using std::vector;
 using std::array;
 using namespace matrix_alloc;
-using namespace jacobi_public_domain;
+using namespace jacobi_pd;
 
 
 // This code works with various types of C++ matrices (for example,