Fix tests

Author: ngoldbaum

src/csrc/scalar.c

@@ -440,21 +440,34 @@ QuadPrecision_conjugate(QuadPrecisionObject *self, PyObject *args)
         return Py_NewRef(self);
     }
 
+    // TODO: Use PyArray_Conjugate when NumPy >= 2.5 is required. NumPy 2.4
+    // can dispatch Quad through a ufunc instead of copying its unchanged value.
+    if (PyArray_FailUnlessWriteable(out, "output array") < 0) {
+        return NULL;
+    }
     PyArray_Descr *dtype = (PyArray_Descr *)new_quaddtype_instance(self->backend);
     if (dtype == NULL) {
         return NULL;
     }
+    if (!PyArray_CanCastTypeTo(dtype, PyArray_DESCR(out), NPY_SAME_KIND_CASTING)) {
+        PyErr_Format(PyExc_TypeError,
+                     "Cannot cast conjugate result from %R to %R with casting='same_kind'",
+                     dtype, PyArray_DESCR(out));
+        Py_DECREF(dtype);
+        return NULL;
+    }
     PyArrayObject *array = (PyArrayObject *)PyArray_SimpleNewFromDescr(0, NULL, dtype);
     if (array == NULL) {
         return NULL;
     }
     quad_value_store(PyArray_BYTES(array), &self->value, self->backend);
-    PyObject *result = PyArray_Conjugate(array, out);
+    int copy_result = PyArray_CopyInto(out, array);
     Py_DECREF(array);
-    if (result == NULL || !PyArray_Check(result)) {
-        return result;
+    if (copy_result < 0) {
+        return NULL;
     }
-    return PyArray_Return((PyArrayObject *)result);
+    Py_INCREF(out);
+    return PyArray_Return(out);
 }
 
 // Method implementations for float compatibility