diff --git a/.github/workflows/fpm.yml b/.github/workflows/fpm.yml
index 9d67128..978351b 100644
--- a/.github/workflows/fpm.yml
+++ b/.github/workflows/fpm.yml
@@ -8,15 +8,12 @@ jobs:
     strategy:
       fail-fast: false
       matrix:
-        os: [ubuntu-latest, macos-11]
+        os: [ubuntu-latest]
         gcc_v: [10] # Version of GFortran we want to use.
         include:
         - os: ubuntu-latest
           os-arch: linux-x86_64
 
-        - os: macos-11
-          os-arch: macos-x86_64
-
     env:
       FC: gfortran
       GCC_V: ${{ matrix.gcc_v }}
diff --git a/CMakeLists.txt b/CMakeLists.txt
index fc71a2f..67c27a0 100644
--- a/CMakeLists.txt
+++ b/CMakeLists.txt
@@ -3,7 +3,7 @@ cmake_minimum_required(VERSION 3.24)
 project(
     nonlin
     LANGUAGES Fortran
-    VERSION 1.5.2
+    VERSION 1.5.3
 )
 
 # Get helper macros and functions
diff --git a/docs/html/index.html b/docs/html/index.html
index 299e493..95b661d 100644
--- a/docs/html/index.html
+++ b/docs/html/index.html
@@ -27,7 +27,7 @@
  <tbody>
  <tr id="projectrow">
   <td id="projectalign">
-   <div id="projectname">nonlin<span id="projectnumber">&#160;1.5.2</span>
+   <div id="projectname">nonlin<span id="projectnumber">&#160;1.5.3</span>
    </div>
    <div id="projectbrief">A library that provides routines to compute the solutions to systems of nonlinear equations.</div>
   </td>
diff --git a/fpm.toml b/fpm.toml
index e362e02..c392f9e 100644
--- a/fpm.toml
+++ b/fpm.toml
@@ -1,5 +1,5 @@
 name = "nonlin"
-version = "1.5.2"
+version = "1.5.3"
 license = "GPL-3.0"
 author = "Jason Christopherson"
 maintainer = "Jason Christopherson"
diff --git a/src/nonlin_least_squares.f90 b/src/nonlin_least_squares.f90
index 5582551..870db5b 100644
--- a/src/nonlin_least_squares.f90
+++ b/src/nonlin_least_squares.f90
@@ -675,7 +675,7 @@ subroutine lmpar(r, ipvt, diag, qtb, delta, par, x, sdiag, wa1, wa2)
             temp = sqrt(par)
             wa1 = temp * diag
             call lmsolve(r(1:n,1:n), ipvt, wa1, qtb, x, sdiag, wa2)
-            wa2 = diag * x
+            wa2(1:n) = diag * x
             dxnorm = norm2(wa2)
             temp = fp
             fp = dxnorm - delta