Add a test

Author: shashi

test/build_function.jl

@@ -16,6 +16,7 @@ h_str = ModelingToolkit.build_function(h, [a], [b], [c1, c2, c3], [d], [e], [g])
 h_oop = eval(h_str[1])
 h_ip! = eval(h_str[2])
 h_ip_skip! = eval(ModelingToolkit.build_function(h, [a], [b], [c1, c2, c3], [d], [e], [g], skipzeros=true, fillzeros=false)[2])
+h_ip_skip_par! = eval(ModelingToolkit.build_function(h, [a], [b], [c1, c2, c3], [d], [e], [g], skipzeros=true, parallel=ModelingToolkit.MultithreadedForm(), fillzeros=false)[2])
 inputs = ([1], [2], [3, 4, 5], [6], [7], [8])
 
 @test h_oop(inputs...) == h_julia(inputs...)
@@ -30,6 +31,12 @@ h_ip_skip!(out_1, inputs...)
 out_1[3] = 0
 @test out_1 == out_2
 
+fill!(out_1, 10)
+h_ip_skip_par!(out_1, inputs...)
+@test out_1[3] == 10
+out_1[3] = 0
+@test out_1 == out_2
+
 # Multiple input matrix, some unused arguments
 h_skip = [a + b + c1; c2 + c3 + g] # skip d, e
 h_julia_skip(a, b, c, d, e, g) = [a[1] + b[1] + c[1]; c[2] + c[3] + g[1]]