LibCEED.jl/test/runtests.jl

23 const run_dev_tests = !isrelease() || ("--run-dev-tests" in ARGS)
89             witharray(x -> x .= 1.0, v)
92             @test summarystr(v) == "$n-element CeedVector"
154             @test q1d ≈ CeedScalar[-1.0, 0.0, 1.0]
158             @test q1d ≈ CeedScalar[-sqrt(3/5), 0.0, sqrt(3/5)]
162             d1d = CeedScalar[-0.5 0.5; -0.5 0.5; -0.5 0.5]
163             q1d = CeedScalar[-1.0, 0.0, 1.0]
232             offsets = Vector{CeedInt}(0:n-1)
233             r = create_elem_restriction(c, 1, n, 1, 1, n, offsets)
238                     (:input, r, b, CeedVectorActive()),
239                     (:output, r, b, CeedVectorActive()),
251             diag_vector = create_lvector(r)
254             # TODO: change this test after bug-fix in libCEED
280             offsets = Vector{CeedInt}([0:n-1; n-1:2*n-2])
281             lsize = 2*n - 1
282             r = create_elem_restriction(c, 2, n, 1, lsize, lsize, offsets)
283             @test getcompstride(r) == lsize
284             @test getnumelements(r) == 2
285             @test getelementsize(r) == n
286             @test getlvectorsize(r) == lsize
287             @test getnumcomponents(r) == 1
288             @test length(create_lvector(r)) == lsize
289             @test length(create_evector(r)) == 2*n
290             lv, ev = create_vectors(r)
293             mult = getmultiplicity(r)
299             @test apply(r, rand_lv) == rand_ev
300             @test apply(r, rand_ev; tmode=TRANSPOSE) == rand_lv.*mult
301             @test showstr(r) == string(
319                  Gallery CeedQFunction - Poisson3DApply
391             d1d = CeedScalar[-0.5 0.5; -0.5 0.5; -0.5 0.5]
392             q1d = CeedScalar[-1.0, 0.0, 1.0]
429                 @test v1[i] ≈ -1.0
457             u2[1] = -1.0
474                 offsets[1, i] = i - 1
477                 orients[1, i] = (i - 1)%2 > 0
478                 orients[2, i] = (i - 1)%2 > 0
480             r = create_elem_restriction_oriented(
493                 lv[i] = 10 + i - 1
496             ev = apply(r, lv)
500                     k = j + elemsize*(i - 1)
501                     @test 10 + k÷2 == ev[k]*(-1)^((i - 1)%2)
508                 if (i - 1)%2 > 0
509                     # T = [0  -1]
510                     #     [-1  0]
512                     curlorients[3, i] = -1
513                     curlorients[4, i] = -1
523             r = create_elem_restriction_curl_oriented(
534             ev = apply(r, lv)
538                     k = j + elemsize*(i - 1)
539                     if (i - 1)%2 > 0
540                         @test j == 2 || 10 + i == -ev[k]
541                         @test j == 1 || 10 + i - 1 == -ev[k]