Changeset 747 for codes/icosagcm/devel/Python/src/unstructured.pyx

Timestamp:

10/05/18 13:42:49 (6 years ago)

Author:

dubos

Message:

devel/unstructured : more fixes to mixed precision

File:

• codes/icosagcm/devel/Python/src/unstructured.pyx (modified) (2 diffs)

codes/icosagcm/devel/Python/src/unstructured.pyx

@@ -239 +239 @@
         # h*s = h => uniform buoyancy s=1 => shallow-water
         dh, du_slow, du_fast, hs, buf = self.fmass(), self.fu(), self.fu(), h.copy(), self.geopot
+        assert type(tau) is np_num, 'tau must be of type unstructured.np_num.'
         ker.dynamico_caldyn_unstructured(tau, self.ms, h, hs, u, self.geopot, buf,
                   self.s, self.ps, self.pk, self.hflux, self.qv,
@@ -296 +297 @@
               left,right,down,up,trisk_deg,trisk,
               Ai, Av, fv, le_de, Riv2, wee):
+
+    print 'Types of Ai, Av, fv, le_de, Riv2, wee : ', Ai.dtype,Av.dtype,fv.dtype,le_de.dtype,Riv2.dtype,wee.dtype
+    for var,varname in zip((Ai,Av,fv,le_de,Riv2,wee), ('Ai','Av','fv','le_de','Riv2','wee')):
+        assert var.dtype == np.float64, '%s must be double precision'%varname
+
     setvars( ('llm','nqdyn','edge_num','primal_num','dual_num',
               'max_trisk_deg','max_primal_deg','max_dual_deg'),