diff --git a/docs.html b/docs.html
index 7de14e0f..781def2e 100644
--- a/docs.html
+++ b/docs.html
@@ -2004,7 +2004,7 @@ All elements are treated as zero by default (ie. the matrix is initialised to co
 </li>
 <br>
 <li>
-Non-zero elements are stored in <a href="http://en.wikipedia.org/wiki/Sparse_matrix">compressed sparse column (CSC) format</a> (ie. column-major ordering);
+Non-zero elements are stored in <a href="https://en.wikipedia.org/wiki/Sparse_matrix">compressed sparse column (CSC) format</a> (ie. column-major ordering);
 zero-valued elements are never stored
 </li>
 <br>
@@ -2036,7 +2036,7 @@ For form&nbsp;4,
 <i>rowind</i> is a dense column vector of type <i>uvec</i> containing the row indices of the values to be inserted,
 and
 <i>colptr</i> is a dense column vector of type <i>uvec</i> (with length <i>n_cols&nbsp;+&nbsp;1</i>) containing indices of <i>values</i> corresponding to the start of new columns;
-the vectors correspond to the arrays used by the <a href="http://en.wikipedia.org/wiki/Sparse_matrix">compressed sparse column format</a>;
+the vectors correspond to the arrays used by the <a href="https://en.wikipedia.org/wiki/Sparse_matrix">compressed sparse column format</a>;
 this form is useful for copying data from other CSC sparse matrix containers
 </li>
 <br>
@@ -2131,7 +2131,7 @@ See also:
 <li><a href="#SpCol">SpCol class</a> (TODO: add to documentation)</li>
 <li><a href="#SpRow">SpRow class</a> (TODO: add to documentation)</li>
 -->
-<li><a href="http://en.wikipedia.org/wiki/Sparse_matrix">Sparse Matrix in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Sparse_matrix">Sparse Matrix in Wikipedia</a></li>
 <li><a href="#Mat">Mat class</a> (dense matrix)</li>
 </ul>
 </li>
@@ -3017,8 +3017,8 @@ See also:
 <li><a href="#ones_member">.ones()</a></li>
 <li><a href="#zeros_member">.zeros()</a></li>
 <li><a href="#size">size()</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Uniform_distribution_(continuous)">uniform distribution in Wikipedia</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Normal_distribution">normal distribution in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Uniform_distribution_(continuous)">uniform distribution in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Normal_distribution">normal distribution in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -3114,8 +3114,8 @@ See also:
 <li><a href="#fill">.fill()</a></li>
 <li><a href="#transform">.transform()</a></li>
 <li><a href="#element_access">element access</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Function_object">function object</a> at Wikipedia</li>
-<li><a href="http://en.wikipedia.org/wiki/C%2B%2B11#Lambda_functions_and_expressions">C++11 lambda functions</a> at Wikipedia</li>
+<li><a href="https://en.wikipedia.org/wiki/Function_object">function object</a> at Wikipedia</li>
+<li><a href="https://en.wikipedia.org/wiki/C%2B%2B11#Lambda_functions_and_expressions">C++11 lambda functions</a> at Wikipedia</li>
 <li><a href="http://www.cprogramming.com/c++11/c++11-lambda-closures.html">lambda function</a> at cprogramming.com</li>
 </ul>
 </li>
@@ -3318,8 +3318,8 @@ See also:
 <li><a href="#element_access">element access</a></li>
 <li><a href="#operators">overloaded operators</a></li>
 <li><a href="#misc_fns">miscellaneous element-wise functions</a> &nbsp; (exp, log, sqrt, square, round, ...)</li>
-<li><a href="http://en.wikipedia.org/wiki/Function_object">function object</a> at Wikipedia</li>
-<li><a href="http://en.wikipedia.org/wiki/C%2B%2B11#Lambda_functions_and_expressions">C++11 lambda functions</a> at Wikipedia</li>
+<li><a href="https://en.wikipedia.org/wiki/Function_object">function object</a> at Wikipedia</li>
+<li><a href="https://en.wikipedia.org/wiki/C%2B%2B11#Lambda_functions_and_expressions">C++11 lambda functions</a> at Wikipedia</li>
 <li><a href="http://www.cprogramming.com/c++11/c++11-lambda-closures.html">lambda function</a> at cprogramming.com</li>
 </ul>
 </li>
@@ -3389,8 +3389,8 @@ See also:
 <li><a href="#each_slice">.each_slice()</a></li>
 <li><a href="#element_access">element access</a></li>
 <li><a href="#misc_fns">miscellaneous element-wise functions</a> &nbsp; (exp, log, sqrt, square, round, ...)</li>
-<li><a href="http://en.wikipedia.org/wiki/Function_object">function object</a> at Wikipedia</li>
-<li><a href="http://en.wikipedia.org/wiki/C%2B%2B11#Lambda_functions_and_expressions">C++11 lambda functions</a> at Wikipedia</li>
+<li><a href="https://en.wikipedia.org/wiki/Function_object">function object</a> at Wikipedia</li>
+<li><a href="https://en.wikipedia.org/wiki/C%2B%2B11#Lambda_functions_and_expressions">C++11 lambda functions</a> at Wikipedia</li>
 <li><a href="http://www.cprogramming.com/c++11/c++11-lambda-closures.html">lambda function</a> at cprogramming.com</li>
 </ul>
 </li>
@@ -6355,7 +6355,7 @@ See also:
 <li><a href="#is_symmetric">.is_symmetric()</a></li>
 <li><a href="#is_diagmat">.is_diagmat()</a></li>
 <li><a href="http://mathworld.wolfram.com/TriangularMatrix.html">Triangular matrix in MathWorld</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Triangular_matrix">Triangular matrix in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Triangular_matrix">Triangular matrix in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -7792,7 +7792,7 @@ arma_rng::set_seed_random();  // set the seed to a random value
 <li><a href="#shuffle">shuffle()</a></li>
 <li><a href="#sprandu_sprandn">sprandu()</a></li>
 <li><a href="#size">size()</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Uniform_distribution_(continuous)">uniform distribution in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Uniform_distribution_(continuous)">uniform distribution in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -7877,7 +7877,7 @@ arma_rng::set_seed_random();  // set the seed to a random value
 <li><a href="#imbue">.imbue()</a></li>
 <li><a href="#sprandu_sprandn">sprandn()</a></li>
 <li><a href="#size">size()</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Normal_distribution">normal distribution in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Normal_distribution">normal distribution in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -7970,7 +7970,7 @@ fmat B2 = randg&lt;fmat&gt;(10, 10, distr_param(2,1));
 <li><a href="#chi2rnd">chi2rnd()</a></li>
 <li><a href="#imbue">.imbue()</a></li>
 <li><a href="#size">size()</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Gamma_distribution">gamma distribution in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Gamma_distribution">gamma distribution in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -8177,8 +8177,8 @@ sp_mat A = sprandu&lt;sp_mat&gt;(100, 200, 0.1);
 <li><a href="#randu">randu()</a></li>
 <li><a href="#randn">randn()</a></li>
 <li><a href="#size">size()</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Uniform_distribution_(continuous)">uniform distribution in Wikipedia</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Normal_distribution">normal distribution in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Uniform_distribution_(continuous)">uniform distribution in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Normal_distribution">normal distribution in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -8216,8 +8216,8 @@ mat Y = circ_toeplitz(A);
 <li>See also:
 <ul>
 <li><a href="http://mathworld.wolfram.com/ToeplitzMatrix.html">Toeplitz matrix in MathWorld</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Toeplitz_matrix">Toeplitz matrix in Wikipedia</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Circulant_matrix">Circulant matrix in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Toeplitz_matrix">Toeplitz matrix in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Circulant_matrix">Circulant matrix in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -8733,7 +8733,7 @@ See also:
 <li><a href="#inv">inv()</a></li>
 <li><a href="#solve">solve()</a></li>
 <li><a href="http://mathworld.wolfram.com/ConditionNumber.html">condition number in MathWorld</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Condition_number">condition number in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Condition_number">condition number in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -8840,7 +8840,7 @@ vec c = cross(a,b);
 See also:
 <ul>
 <li><a href="#dot">dot()</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Cross_product">Cross product in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Cross_product">Cross product in Wikipedia</a></li>
 <li><a href="http://mathworld.wolfram.com/CrossProduct.html">Cross product in MathWorld</a></li>
 </ul>
 </li>
@@ -8979,7 +8979,7 @@ See also:
 <li><a href="#log_det">log_det()</a></li>
 <li><a href="#rcond">rcond()</a></li>
 <li><a href="http://mathworld.wolfram.com/Determinant.html">determinant in MathWorld</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Determinant">determinant in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Determinant">determinant in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -9134,7 +9134,7 @@ vec b = diff(a);
 <li>See also:
 <ul>
 <li><a href="#trapz">trapz()</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Numerical_differentiation">numerical differentiation in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Numerical_differentiation">numerical differentiation in Wikipedia</a></li>
 <li><a href="http://mathworld.wolfram.com/NumericalDifferentiation.html">numerical differentiation in MathWorld</a></li>
 </ul>
 </li>
@@ -9213,7 +9213,7 @@ See also:
 <ul>
 <li><a href="#constants">datum::eps</a></li>
 <li><a href="http://mathworld.wolfram.com/Floating-PointArithmetic.html">Floating-Point Arithmetic in MathWorld</a></li>
-<li><a href="http://en.wikipedia.org/wiki/IEEE_754-2008">IEEE Standard for Floating-Point Arithmetic in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/IEEE_754-2008">IEEE Standard for Floating-Point Arithmetic in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -9257,7 +9257,7 @@ See also:
 <li><a href="#logmat">logmat()</a></li>
 <li><a href="#sqrtmat">sqrtmat()</a></li>
 <li><a href="#misc_fns">miscellaneous element-wise functions</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Matrix_exponential">matrix exponential in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Matrix_exponential">matrix exponential in Wikipedia</a></li>
 <li><a href="http://mathworld.wolfram.com/MatrixExponential.html">matrix exponential in MathWorld</a></li>
 <li><a href="http://doi.org/10.1137/S00361445024180">Nineteen Dubious Ways to Compute the Exponential of a Matrix, Twenty-Five Years Later</a></li>
 </ul>
@@ -9311,7 +9311,7 @@ See also:
 <li><a href="#is_symmetric">.is_symmetric()</a></li>
 <li><a href="#is_hermitian">.is_hermitian()</a></li>
 <li><a href="#misc_fns">miscellaneous element-wise functions</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Matrix_exponential">matrix exponential in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Matrix_exponential">matrix exponential in Wikipedia</a></li>
 <li><a href="http://mathworld.wolfram.com/MatrixExponential.html">matrix exponential in MathWorld</a></li>
 <li><a href="http://doi.org/10.1137/S00361445024180">Nineteen Dubious Ways to Compute the Exponential of a Matrix, Twenty-Five Years Later</a></li>
 </ul>
@@ -9532,7 +9532,7 @@ See also:
 <li><a href="#replace">.replace()</a></li>
 <li><a href="#has_nan">.has_nan()</a></li>
 <li><a href="#constants">constants (pi, nan, inf, ...)</a></li>
-<li><a href="http://en.wikipedia.org/wiki/NaN">NaN</a> in Wikipedia</li>
+<li><a href="https://en.wikipedia.org/wiki/NaN">NaN</a> in Wikipedia</li>
 </ul>
 </li>
 <br>
@@ -10096,7 +10096,7 @@ See also:
 <li><a href="#repmat">repmat()</a></li>
 <li><a href="#repelem">repelem()</a></li>
 <li><a href="http://mathworld.wolfram.com/KroneckerProduct.html">Kronecker product in MathWorld</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Kronecker_product">Kronecker product in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Kronecker_product">Kronecker product in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -10168,7 +10168,7 @@ See also:
 <li><a href="#rcond">rcond()</a></li>
 <li><a href="#cx_double">cx_double</a></li>
 <li><a href="http://mathworld.wolfram.com/Determinant.html">determinant in MathWorld</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Determinant">determinant in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Determinant">determinant in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -10220,7 +10220,7 @@ See also:
 <li><a href="#log_det">log_det()</a></li>
 <li><a href="#rcond">rcond()</a></li>
 <li><a href="http://mathworld.wolfram.com/Determinant.html">determinant in MathWorld</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Determinant">determinant in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Determinant">determinant in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -10514,9 +10514,9 @@ See also:
 <li><a href="#normalise">normalise()</a></li>
 <li><a href="#vectorise">vectorise()</a></li>
 <li><a href="#dot">dot()</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Norm_(mathematics)">vector norm in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Norm_(mathematics)">vector norm in Wikipedia</a></li>
 <li><a href="http://mathworld.wolfram.com/VectorNorm.html">vector norm in MathWorld</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Matrix_norm">matrix norm in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Matrix_norm">matrix norm in Wikipedia</a></li>
 <li><a href="http://mathworld.wolfram.com/MatrixNorm.html">matrix norm in MathWorld</a></li>
 </ul>
 </li>
@@ -10566,7 +10566,7 @@ See also:
 <li><a href="#norm">norm()</a></li>
 <li><a href="#dot">norm_dot()</a></li>
 <li><a href="http://mathworld.wolfram.com/NormalizedVector.html">Normalised vector in MathWorld</a>
-<li><a href="http://en.wikipedia.org/wiki/Unit_vector">Unit vector in Wikipedia</a>
+<li><a href="https://en.wikipedia.org/wiki/Unit_vector">Unit vector in Wikipedia</a>
 </ul>
 </li>
 <br>
@@ -10803,7 +10803,7 @@ See also:
 <li><a href="#orth">orth()</a></li>
 <li><a href="#constants">datum::eps</a></li>
 <li><a href="http://mathworld.wolfram.com/MatrixRank.html">Rank in MathWorld</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Rank_(linear_algebra)">Rank in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Rank_(linear_algebra)">Rank in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -11656,7 +11656,7 @@ mat C = symmatl(A);
 <li><a href="#trimat">trimatu() / trimatl()</a></li>
 <li><a href="#is_symmetric">.is_symmetric()</a></li>
 <li><a href="#is_hermitian">.is_hermitian()</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Symmetric_matrix">Symmetric matrix in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Symmetric_matrix">Symmetric matrix in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -11854,7 +11854,7 @@ mat LL = trimatl(A, -1);  // omit the main diagonal
 <li><a href="#diagmat">diagmat()</a></li>
 <li><a href="#nonzeros">nonzeros()</a></li>
 <li><a href="http://mathworld.wolfram.com/TriangularMatrix.html">Triangular matrix in MathWorld</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Triangular_matrix">Triangular matrix in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Triangular_matrix">Triangular matrix in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -12582,8 +12582,8 @@ See also:
 <li><a href="#is_symmetric">.is_symmetric()</a></li>
 <li><a href="#is_hermitian">.is_hermitian()</a></li>
 <li><a href="http://mathworld.wolfram.com/EigenDecomposition.html">eigen decomposition in MathWorld</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Eigenvalues_and_eigenvectors">eigenvalues &amp; eigenvectors in Wikipedia</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Divide-and-conquer_eigenvalue_algorithm">divide &amp; conquer eigenvalue algorithm in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Eigenvalues_and_eigenvectors">eigenvalues &amp; eigenvectors in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Divide-and-conquer_eigenvalue_algorithm">divide &amp; conquer eigenvalue algorithm in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -12659,7 +12659,7 @@ See also:
 <li><a href="#schur">schur()</a></li>
 <li><a href="#eigs_gen">eigs_gen()</a></li>
 <li><a href="http://mathworld.wolfram.com/EigenDecomposition.html">eigen decomposition in MathWorld</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Eigenvalues_and_eigenvectors">eigenvalues &amp; eigenvectors in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Eigenvalues_and_eigenvectors">eigenvalues &amp; eigenvectors in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -12719,7 +12719,7 @@ See also:
 <li><a href="#eig_sym">eig_sym()</a></li>
 <li><a href="#qz">qz()</a></li>
 <li><a href="http://mathworld.wolfram.com/EigenDecomposition.html">eigen decomposition in MathWorld</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Eigenvalues_and_eigenvectors">eigenvalues &amp; eigenvectors in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Eigenvalues_and_eigenvectors">eigenvalues &amp; eigenvectors in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -12870,7 +12870,7 @@ See also:
 <li><a href="#trimat">trimatu() / trimatl()</a></li>
 <li><a href="#powmat">powmat()</a></li>
 <li><a href="http://mathworld.wolfram.com/MatrixInverse.html">matrix inverse in MathWorld</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Invertible_matrix">invertible matrix in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Invertible_matrix">invertible matrix in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -12966,9 +12966,9 @@ See also:
 <li><a href="#eig_sym">eig_sym()</a></li>
 <li><a href="#is_sympd">.is_sympd()</a></li>
 <li><a href="http://mathworld.wolfram.com/MatrixInverse.html">matrix inverse in MathWorld</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Invertible_matrix">invertible matrix in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Invertible_matrix">invertible matrix in Wikipedia</a></li>
 <li><a href="http://mathworld.wolfram.com/PositiveDefiniteMatrix.html">positive definite matrix in MathWorld</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Positive-definite_matrix">positive definite matrix in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Positive-definite_matrix">positive definite matrix in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -13023,7 +13023,7 @@ mat B = P.t()*L*U;
 See also:
 <ul>
 <li><a href="#chol">chol()</a></li>
-<li><a href="http://en.wikipedia.org/wiki/LU_decomposition">LU decomposition in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/LU_decomposition">LU decomposition in Wikipedia</a></li>
 <li><a href="http://mathworld.wolfram.com/LUDecomposition.html">LU decomposition in MathWorld</a></li>
 </ul>
 </li>
@@ -13085,7 +13085,7 @@ See also:
 <li><a href="#svd">svd()</a></li>
 <li><a href="#rank">rank()</a></li>
 <li><a href="#constants">datum::eps</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Orthonormal_basis">Orthonormal basis in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Orthonormal_basis">Orthonormal basis in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -13144,7 +13144,7 @@ See also:
 <li><a href="#svd">svd()</a></li>
 <li><a href="#rank">rank()</a></li>
 <li><a href="#constants">datum::eps</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Orthonormal_basis">Orthonormal basis in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Orthonormal_basis">Orthonormal basis in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -13287,8 +13287,8 @@ See also:
 <li><a href="#qr_econ">qr_econ()</a></li>
 <li><a href="#chol">chol()</a></li>
 <li><a href="#orth">orth()</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Orthogonal_matrix">orthogonal matrix in Wikipedia</a></li>
-<li><a href="http://en.wikipedia.org/wiki/QR_decomposition">QR decomposition in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Orthogonal_matrix">orthogonal matrix in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/QR_decomposition">QR decomposition in Wikipedia</a></li>
 <li><a href="http://mathworld.wolfram.com/QRDecomposition.html">QR decomposition in MathWorld</a></li>
 </ul>
 </li>
@@ -13330,8 +13330,8 @@ qr_econ(Q, R, X);
 See also:
 <ul>
 <li><a href="#qr">qr()</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Orthogonal_matrix">orthogonal matrix in Wikipedia</a></li>
-<li><a href="http://en.wikipedia.org/wiki/QR_decomposition">QR decomposition in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Orthogonal_matrix">orthogonal matrix in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/QR_decomposition">QR decomposition in Wikipedia</a></li>
 <li><a href="http://octave.sourceforge.net/octave/function/qr.html">QR decomposition in Octave</a></li>
 <li><a href="http://mathworld.wolfram.com/QRDecomposition.html">QR decomposition in MathWorld</a></li>
 </ul>
@@ -13578,7 +13578,7 @@ See also:
 <li><a href="#syl">syl()</a></li>
 <li><a href="#spsolve">spsolve()</a> - solve sparse system of linear equations</li>
 <li><a href="http://mathworld.wolfram.com/LinearSystemofEquations.html">linear system of equations in MathWorld</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Linear_system_of_equations">system of linear equations in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Linear_system_of_equations">system of linear equations in Wikipedia</a></li>
 <li><a href="https://en.wikipedia.org/wiki/Band_matrix">band matrix in Wikipedia</a></li>
 <li><a href="https://en.wikipedia.org/wiki/Definiteness_of_a_matrix">definiteness of a matrix in Wikipedia</a></li>
 <li><a href="http://mathworld.wolfram.com/PositiveDefiniteMatrix.html">positive definite matrix in MathWorld</a></li>
@@ -13658,7 +13658,7 @@ See also:
 <li><a href="#eig_sym">eig_sym()</a></li>
 <li><a href="#princomp">princomp()</a></li>
 <li><a href="#svds">svds()</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Singular_value_decomposition">singular value decomposition in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Singular_value_decomposition">singular value decomposition in Wikipedia</a></li>
 <li><a href="http://mathworld.wolfram.com/SingularValueDecomposition.html">singular value decomposition in MathWorld</a></li>
 </ul>
 </li>
@@ -13738,7 +13738,7 @@ See also:
 <li><a href="#eig_sym">eig_sym()</a></li>
 <li><a href="#princomp">princomp()</a></li>
 <li><a href="#svds">svds()</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Singular_value_decomposition">singular value decomposition in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Singular_value_decomposition">singular value decomposition in Wikipedia</a></li>
 <li><a href="http://mathworld.wolfram.com/SingularValueDecomposition.html">singular value decomposition in MathWorld</a></li>
 </ul>
 </li>
@@ -13782,7 +13782,7 @@ syl(X2, A, B, C);
 See also:
 <ul>
 <li><a href="#solve">solve()</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Sylvester_equation">Sylvester equation in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Sylvester_equation">Sylvester equation in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -13917,7 +13917,7 @@ See also:
 <li><a href="#is_symmetric">.is_symmetric()</a></li>
 <li><a href="https://www.caam.rice.edu/software/ARPACK/UG/node33.html">shift-invert mode in ARPACK</a></li>
 <li><a href="http://mathworld.wolfram.com/EigenDecomposition.html">eigen decomposition in MathWorld</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Eigenvalues_and_eigenvectors">eigenvalues &amp; eigenvectors in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Eigenvalues_and_eigenvectors">eigenvalues &amp; eigenvectors in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -14045,7 +14045,7 @@ See also:
 <li><a href="#svds">svds()</a></li>
 <li><a href="https://www.caam.rice.edu/software/ARPACK/UG/node33.html">shift-invert mode in ARPACK</a></li>
 <li><a href="http://mathworld.wolfram.com/EigenDecomposition.html">eigen decomposition in MathWorld</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Eigenvalues_and_eigenvectors">eigenvalues &amp; eigenvectors in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Eigenvalues_and_eigenvectors">eigenvalues &amp; eigenvectors in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -14194,7 +14194,7 @@ See also:
 <!-- <li><a href="http://crd-legacy.lbl.gov/~xiaoye/SuperLU/">SuperLU home page</a> -->
 <li><a href="https://portal.nersc.gov/project/sparse/superlu/">SuperLU home page</a>
 <li><a href="http://mathworld.wolfram.com/LinearSystemofEquations.html">linear system of equations in MathWorld</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Linear_system_of_equations">system of linear equations in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Linear_system_of_equations">system of linear equations in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -14286,7 +14286,7 @@ See also:
 <li><a href="#eigs_gen">eigs_gen()</a></li>
 <li><a href="#eigs_sym">eigs_sym()</a></li>
 <li><a href="#svd">svd()</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Singular_value_decomposition">singular value decomposition in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Singular_value_decomposition">singular value decomposition in Wikipedia</a></li>
 <li><a href="http://mathworld.wolfram.com/SingularValueDecomposition.html">singular value decomposition in MathWorld</a></li>
 </ul>
 </li>
@@ -14360,8 +14360,8 @@ See also:
 <li><a href="#cor">cor()</a></li>
 <li><a href="#interp1">interp1()</a></li>
 <li><a href="http://mathworld.wolfram.com/Convolution.html">Convolution in MathWorld</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Convolution">Convolution in Wikipedia</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Finite_impulse_response">FIR filter in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Convolution">Convolution in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Finite_impulse_response">FIR filter in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -14414,8 +14414,8 @@ See also:
 <li><a href="#fft2">fft2()</a></li>
 <li><a href="#interp2">interp2()</a></li>
 <li><a href="http://mathworld.wolfram.com/Convolution.html">Convolution in MathWorld</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Convolution">Convolution in Wikipedia</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Kernel_%28image_processing%29">Kernel (image processing) in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Convolution">Convolution in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Kernel_%28image_processing%29">Kernel (image processing) in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -14469,7 +14469,7 @@ See also:
 <li><a href="#conv">conv()</a></li>
 <li><a href="#imag_real">real()</a></li>
 <li><a href="http://mathworld.wolfram.com/FastFourierTransform.html">fast Fourier transform in MathWorld</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Fast_Fourier_transform">fast Fourier transform in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Fast_Fourier_transform">fast Fourier transform in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -14515,7 +14515,7 @@ See also:
 <li><a href="#conv2">conv2()</a></li>
 <li><a href="#imag_real">real()</a></li>
 <li><a href="http://mathworld.wolfram.com/FastFourierTransform.html">fast Fourier transform in MathWorld</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Fast_Fourier_transform">fast Fourier transform in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Fast_Fourier_transform">fast Fourier transform in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -15353,7 +15353,7 @@ See also:
 <li><a href="#eig_sym">eig_sym()</a></li>
 <li><a href="#svd">svd()</a></li>
 <li><a href="#svd_econ">svd_econ()</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Principal_component_analysis">principal components analysis in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Principal_component_analysis">principal components analysis in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -15421,7 +15421,7 @@ See also:
 <li><a href="#normcdf">normcdf()</a></li>
 <li><a href="#randn">randn()</a></li>
 <li><a href="#gmm_diag">gmm_diag&nbsp;/&nbsp;gmm_full</a> - model and evaluate data using Gaussian Mixture Models (GMMs)</li>
-<li><a href="http://en.wikipedia.org/wiki/Normal_distribution">normal distribution in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Normal_distribution">normal distribution in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -15497,7 +15497,7 @@ See also:
 <ul>
 <li><a href="#normpdf">normpdf()</a></li>
 <li><a href="#gmm_diag">gmm_diag&nbsp;/&nbsp;gmm_full</a> - model and evaluate data using Gaussian Mixture Models (GMMs)</li>
-<li><a href="http://en.wikipedia.org/wiki/Normal_distribution">normal distribution in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Normal_distribution">normal distribution in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -15541,7 +15541,7 @@ See also:
 <li><a href="#normpdf">normpdf()</a></li>
 <li><a href="#quantile">quantile()</a></li>
 <li><a href="#randn">randn()</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Normal_distribution">normal distribution in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Normal_distribution">normal distribution in Wikipedia</a></li>
 <li><a href="https://en.wikipedia.org/wiki/Cumulative_distribution_function">cumulative distribution function in Wikipedia</a></li>
 </ul>
 </li>
@@ -16245,9 +16245,9 @@ means.print("means:");
 <li><a href="#gmm_diag">gmm_diag&nbsp;/&nbsp;gmm_full</a> - model and evaluate data using Gaussian Mixture Models (GMMs)</li>
 <li><a href="#stats_fns">statistics functions</a></li>
 <li><a href="#running_stat_vec">running_stat_vec</a></li>
-<li><a href="http://en.wikipedia.org/wiki/K-means_clustering">k-means clustering in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/K-means_clustering">k-means clustering in Wikipedia</a></li>
 <li><a href="http://mathworld.wolfram.com/K-MeansClusteringAlgorithm.html">k-means clustering in MathWorld</a></li>
-<li><a href="http://en.wikipedia.org/wiki/OpenMP">OpenMP in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/OpenMP">OpenMP in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -16930,15 +16930,15 @@ model.save("my_model.gmm");
 <li><a href="#kmeans">kmeans()</a></li>
 <li><a href="https://en.wikipedia.org/wiki/Covariance_matrix">covariance matrix in Wikipedia</a></li>
 <li><a href="http://mathworld.wolfram.com/CovarianceMatrix.html">covariance matrix in MathWorld</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Mahalanobis_distance">Mahalanobis distance in Wikipedia</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Multivariate_normal_distribution">multivariate normal distribution in Wikipedia</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Mixture_model">mixture model in Wikipedia</a></li>
-<li><a href="http://en.wikipedia.org/wiki/K-means_clustering">k-means clustering in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Mahalanobis_distance">Mahalanobis distance in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Multivariate_normal_distribution">multivariate normal distribution in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Mixture_model">mixture model in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/K-means_clustering">k-means clustering in Wikipedia</a></li>
 <li><a href="http://mathworld.wolfram.com/K-MeansClusteringAlgorithm.html">k-means clustering in MathWorld</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Expectation-maximization_algorithm">Expectation-Maximisation algorithm in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Expectation-maximization_algorithm">Expectation-Maximisation algorithm in Wikipedia</a></li>
 <li><a href="http://mathworld.wolfram.com/MaximumLikelihood.html">maximum likelihood in MathWorld</a></li>
-<li><a href="http://en.wikipedia.org/wiki/Vector_quantization">vector quantisation in Wikipedia</a></li>
-<li><a href="http://en.wikipedia.org/wiki/OpenMP">OpenMP in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/Vector_quantization">vector quantisation in Wikipedia</a></li>
+<li><a href="https://en.wikipedia.org/wiki/OpenMP">OpenMP in Wikipedia</a></li>
 </ul>
 </li>
 <br>
@@ -17524,8 +17524,8 @@ See also:
 <ul>
 <li><a href="#is_finite">.is_finite()</a></li>
 <li><a href="#fill">.fill()</a></li>
-<li><a href="http://en.wikipedia.org/wiki/NaN">NaN</a> in Wikipedia</li>
-<li><a href="http://en.wikipedia.org/wiki/Physical_constant">physical constant</a> in Wikipedia</li>
+<li><a href="https://en.wikipedia.org/wiki/NaN">NaN</a> in Wikipedia</li>
+<li><a href="https://en.wikipedia.org/wiki/Physical_constant">physical constant</a> in Wikipedia</li>
 <li><a href="https://en.wikipedia.org/wiki/Turn_(angle)">replacement of 2&pi; with &tau;</a> in Wikipedia</li>
 <li><a href="https://tauday.com/tau-manifesto">The Tau Manifesto</a> by Michael Hartl</li>
 <li><a href="http://cplusplus.com/reference/std/limits/numeric_limits/">std::numeric_limits</a> in cplusplus.com</li>