add prod tests in multiprecision (#108)

Author: tmigot

src/nlp/multiple-precision.jl

@@ -47,10 +47,12 @@ function multiple_precision_nlp(
     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)
     @test obj ∈ exclude || typeof(obj(nlp, x)) == T
     @test grad ∈ exclude || typeof(grad(nlp, x)) == S
     @test hess ∈ exclude || eltype(hess(nlp, x)) == T
     @test hess_op ∈ exclude || eltype(hess_op(nlp, x)) == T
+    @test hprod ∈ exclude || typeof(hprod(nlp, x, v)) == S
     if hess_coord ∉ exclude && hess_op ∉ exclude
       rows, cols = hess_structure(nlp)
       vals = hess_coord(nlp, x)
@@ -63,6 +65,8 @@ function multiple_precision_nlp(
       @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
+      @test jprod ∈ exclude || typeof(jprod(nlp, x, v)) == S
+      @test jtprod ∈ exclude || typeof(jtprod(nlp, x, y)) == S
       if linear_api && nlp.meta.nnln > 0
         @test cons ∈ exclude || typeof(cons_nln(nlp, x)) == S
         @test jac ∈ exclude || eltype(jac_nln(nlp, x)) == T
@@ -100,6 +104,7 @@ function multiple_precision_nlp(
       @test hess ∈ exclude || eltype(hess(nlp, x, y)) == T
       @test hess ∈ exclude || eltype(hess(nlp, x, y, obj_weight = one(T))) == T
       @test hess_op ∈ exclude || eltype(hess_op(nlp, x, y)) == T
+      @test hprod ∈ exclude || typeof(hprod(nlp, x, y, v)) == S
       if hess_coord ∉ exclude && hess_op ∉ exclude
         rows, cols = hess_structure(nlp)
         vals = hess_coord(nlp, x, y)

src/nls/multiple-precision.jl

@@ -47,9 +47,14 @@ function multiple_precision_nls(
     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)
     @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
+    @test jprod_residual ∈ exclude || typeof(jprod_residual(nls, x, v)) == S
+    @test jtprod_residual ∈ exclude || typeof(jtprod_residual(nls, x, w)) == S
     if jac_coord_residual ∉ exclude && jac_op_residual ∉ exclude
       rows, cols = jac_structure_residual(nls)
       vals = jac_coord_residual(nls, x)
@@ -64,6 +69,11 @@ function multiple_precision_nls(
         @test eltype(hess_op_residual(nls, x, i)) == T
       end
     end
+    if hprod_residual ∉ exclude
+      for i = 1:(nls.nls_meta.nequ)
+        @test typeof(hprod_residual(nls, x, i, v)) == S
+      end
+    end
     @test obj ∈ exclude || typeof(obj(nls, x)) == T
     @test grad ∈ exclude || typeof(grad(nls, x)) == S
     if nls.meta.ncon > 0
@@ -87,6 +97,8 @@ function multiple_precision_nls(
         Av = fill!(S(undef, nls.meta.ncon), 0)
         Atv = fill!(S(undef, nls.meta.nvar), 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)