Do not fill-in with integer in mutli-precision tests (#111)

Author: tmigot

src/nlp/multiple-precision.jl

@@ -46,8 +46,8 @@ function multiple_precision_nlp(
   for T in precisions
     nlp = nlp_from_T(T)
     S = typeof(nlp.meta.x0)
-    x = fill!(S(undef, nlp.meta.nvar), 1)
-    v = fill!(S(undef, nlp.meta.nvar), 2)
+    x = fill!(S(undef, nlp.meta.nvar), T(1))
+    v = fill!(S(undef, nlp.meta.nvar), T(2))
     @test obj ∈ exclude || typeof(obj(nlp, x)) == T
     @test grad ∈ exclude || typeof(grad(nlp, x)) == S
     @test hess ∈ exclude || eltype(hess(nlp, x)) == T
@@ -57,11 +57,11 @@ function multiple_precision_nlp(
       rows, cols = hess_structure(nlp)
       vals = hess_coord(nlp, x)
       @test typeof(vals) == S
-      Hv = fill!(S(undef, nlp.meta.nvar), 1)
+      Hv = fill!(S(undef, nlp.meta.nvar), T(1))
       @test eltype(hess_op!(nlp, rows, cols, vals, Hv)) == T
     end
     if nlp.meta.ncon > 0
-      y = fill!(S(undef, nlp.meta.ncon), 1)
+      y = fill!(S(undef, nlp.meta.ncon), T(1))
       @test cons ∈ exclude || typeof(cons(nlp, x)) == S
       @test jac ∈ exclude || eltype(jac(nlp, x)) == T
       @test jac_op ∈ exclude || eltype(jac_op(nlp, x)) == T
@@ -81,23 +81,23 @@ function multiple_precision_nlp(
         rows, cols = jac_structure(nlp)
         vals = jac_coord(nlp, x)
         @test typeof(vals) == S
-        Av = fill!(S(undef, nlp.meta.ncon), 0)
-        Atv = fill!(S(undef, nlp.meta.nvar), 0)
+        Av = fill!(S(undef, nlp.meta.ncon), T(0))
+        Atv = fill!(S(undef, nlp.meta.nvar), T(0))
         @test eltype(jac_op!(nlp, rows, cols, vals, Av, Atv)) == T
         if linear_api && nlp.meta.nnln > 0
           rows, cols = jac_nln_structure(nlp)
           vals = jac_nln_coord(nlp, x)
           @test typeof(vals) == S
-          Av = fill!(S(undef, nlp.meta.nnln), 0)
-          Atv = fill!(S(undef, nlp.meta.nvar), 0)
+          Av = fill!(S(undef, nlp.meta.nnln), T(0))
+          Atv = fill!(S(undef, nlp.meta.nvar), T(0))
           @test eltype(jac_nln_op!(nlp, rows, cols, vals, Av, Atv)) == T
         end
         if linear_api && nlp.meta.nlin > 0
           rows, cols = jac_lin_structure(nlp)
           vals = jac_lin_coord(nlp, x)
           @test typeof(vals) == S
-          Av = fill!(S(undef, nlp.meta.nlin), 0)
-          Atv = fill!(S(undef, nlp.meta.nvar), 0)
+          Av = fill!(S(undef, nlp.meta.nlin), T(0))
+          Atv = fill!(S(undef, nlp.meta.nvar), T(0))
           @test eltype(jac_lin_op!(nlp, rows, cols, vals, Av, Atv)) == T
         end
       end
@@ -109,7 +109,7 @@ function multiple_precision_nlp(
         rows, cols = hess_structure(nlp)
         vals = hess_coord(nlp, x, y)
         @test typeof(vals) == S
-        Hv = fill!(S(undef, nlp.meta.nvar), 0)
+        Hv = fill!(S(undef, nlp.meta.nvar), T(0))
         @test eltype(hess_op!(nlp, rows, cols, vals, Hv)) == T
       end
       @test jth_hess ∈ exclude || eltype(jth_hess(nlp, x, 1)) == T

src/nls/multiple-precision.jl

@@ -46,10 +46,10 @@ function multiple_precision_nls(
   for T in precisions
     nls = nls_from_T(T)
     S = typeof(nls.meta.x0)
-    x = fill!(S(undef, nls.meta.nvar), 1)
-    v = fill!(S(undef, nls.meta.nvar), 2)
-    y = fill!(S(undef, nls.meta.ncon), 2)
-    w = fill!(S(undef, nls.nls_meta.nequ), 2)
+    x = fill!(S(undef, nls.meta.nvar), T(1))
+    v = fill!(S(undef, nls.meta.nvar), T(2))
+    y = fill!(S(undef, nls.meta.ncon), T(2))
+    w = fill!(S(undef, nls.nls_meta.nequ), T(2))
     @test residual ∈ exclude || typeof(residual(nls, x)) == S
     @test jac_residual ∈ exclude || eltype(jac_residual(nls, x)) == T
     @test jac_op_residual ∈ exclude || eltype(jac_op_residual(nls, x)) == T
@@ -59,8 +59,8 @@ function multiple_precision_nls(
       rows, cols = jac_structure_residual(nls)
       vals = jac_coord_residual(nls, x)
       @test typeof(vals) == S
-      Av = fill!(S(undef, nls.nls_meta.nequ), 0)
-      Atv = fill!(S(undef, nls.meta.nvar), 0)
+      Av = fill!(S(undef, nls.nls_meta.nequ), T(0))
+      Atv = fill!(S(undef, nls.meta.nvar), T(0))
       @test eltype(jac_op_residual!(nls, rows, cols, vals, Av, Atv)) == T
     end
     @test hess_residual ∈ exclude || eltype(hess_residual(nls, x, ones(T, nls.nls_meta.nequ))) == T
@@ -94,25 +94,25 @@ function multiple_precision_nls(
         rows, cols = jac_structure(nls)
         vals = jac_coord(nls, x)
         @test typeof(vals) == S
-        Av = fill!(S(undef, nls.meta.ncon), 0)
-        Atv = fill!(S(undef, nls.meta.nvar), 0)
+        Av = fill!(S(undef, nls.meta.ncon), T(0))
+        Atv = fill!(S(undef, nls.meta.nvar), T(0))
         @test eltype(jac_op!(nls, rows, cols, vals, Av, Atv)) == T
         @test jprod ∈ exclude || typeof(jprod(nls, x, v)) == S
         @test jtprod ∈ exclude || typeof(jtprod(nls, x, y)) == S
         if linear_api && nls.meta.nnln > 0
           rows, cols = jac_nln_structure(nls)
           vals = jac_nln_coord(nls, x)
           @test typeof(vals) == S
-          Av = fill!(S(undef, nls.meta.nnln), 0)
-          Atv = fill!(S(undef, nls.meta.nvar), 0)
+          Av = fill!(S(undef, nls.meta.nnln), T(0))
+          Atv = fill!(S(undef, nls.meta.nvar), T(0))
           @test eltype(jac_nln_op!(nls, rows, cols, vals, Av, Atv)) == T
         end
         if linear_api && nls.meta.nlin > 0
           rows, cols = jac_lin_structure(nls)
           vals = jac_lin_coord(nls, x)
           @test typeof(vals) == S
-          Av = fill!(S(undef, nls.meta.nlin), 0)
-          Atv = fill!(S(undef, nls.meta.nvar), 0)
+          Av = fill!(S(undef, nls.meta.nlin), T(0))
+          Atv = fill!(S(undef, nls.meta.nvar), T(0))
           @test eltype(jac_lin_op!(nls, rows, cols, vals, Av, Atv)) == T
         end
       end