Set RHS to same as A matrix

domfournier · domfournier · commit 57a2f2320464 · 2026-04-10T10:34:11.000-07:00
diff --git a/simpeg/potential_fields/magnetics/simulation.py b/simpeg/potential_fields/magnetics/simulation.py
@@ -1862,7 +1862,7 @@ def _getRHS(self, m):
                 ).diagonal()
             )
 
-        return rhs
+        return rhs.astype(self.solver_dtype)
 
     def _getA(self):
         A = self._Div * self.MfMuiI * self._DivT

Original file line number	Diff line number	Diff line change
`@@ -1862,7 +1862,7 @@ def _getRHS(self, m):`
`1862`	`1862`	`).diagonal()`
`1863`	`1863`	`)`
`1864`	`1864`
`1865`		`- return rhs`
	`1865`	`+ return rhs.astype(self.solver_dtype)`
`1866`	`1866`
`1867`	`1867`	`def _getA(self):`
`1868`	`1868`	`A = self._Div * self.MfMuiI * self._DivT`