Lines Matching full:1
69 @test axes(v) == (1:n,)
70 @test ndims(v) == 1
71 @test ndims(CeedVector) == 1
81 for p ∈ [1, 2, Inf]
140 ncomp = 1
155 @test w1d ≈ CeedScalar[1/3, 4/3, 1/3]
164 w1d = CeedScalar[1/3, 4/3, 1/3]
167 d2d[1, :, :] = kron(b1d, d1d)
171 b2 = create_tensor_h1_basis(c, dim2, 1, p, q, b1d, d1d, q1d, w1d)
180 1,
184 reshape(d1d, 1, q, p),
185 reshape(q1d, 1, q),
188 @test getqref(b3) == reshape(q1d, 1, q)
204 for dim = 1:3
211 V = zeros(CeedScalar, dim*(dim + 1)÷2)
230 b = create_tensor_h1_lagrange_basis(c, 3, 1, 3, 3, GAUSS_LOBATTO)
232 offsets = Vector{CeedInt}(0:n-1)
233 r = create_elem_restriction(c, 1, n, 1, 1, n, offsets)
261 1 elements with 27 quadrature points each
263 1 input field:
266 Size: 1
269 1 output field:
272 Size: 1
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)
287 @test getnumcomponents(r) == 1
298 rand_ev = [rand_lv[1:n]; rand_lv[n:end]]
302 "CeedElemRestriction from (19, 1) to 2 elements ",
306 strides = CeedInt[1, n, n]
307 rs = create_elem_restriction_strided(c, 1, n, 1, n, strides)
309 "CeedElemRestriction from (10, 1) to 1 elements ",
310 "with 10 nodes each and strides [1, $n, $n]",
325 Input field 1:
329 1 output field:
335 id = create_identity_qfunction(c, 1, EVAL_INTERP, EVAL_INTERP)
365 @test in_sz == [dim, 1]
366 @test out_sz == [1]
368 v2 = rand(CeedScalar, 1)
371 cv3 = CeedVector(c, 1)
372 apply!(qf, 1, [cv1, cv2], [cv3])
374 @test @witharray_read(v3 = cv3, v3[1] == v2[1]*sum(v1))
383 ncomp = 1
393 w1d = CeedScalar[1/3, 4/3, 1/3]
396 d2d[1, :, :] = kron(b1d, d1d)
400 b2 = create_tensor_h1_basis(c, dim2, 1, p, q, b1d, d1d, q1d, w1d)
406 1,
410 reshape(d1d, 1, q, p),
411 reshape(q1d, 1, q),
417 ncomp = 1
428 for i = 1:q1
453 for i = 1:q2
457 u2[1] = -1.0
462 for i = 1:q2
473 for i = 1:nelem
474 offsets[1, i] = i - 1
476 # flip the dofs on element 1, 3, ...
477 orients[1, i] = (i - 1)%2 > 0
478 orients[2, i] = (i - 1)%2 > 0
484 1,
485 1,
486 nelem + 1,
491 lv = Vector{CeedScalar}(undef, nelem + 1)
492 for i = 1:nelem+1
493 lv[i] = 10 + i - 1
498 for i = 1:nelem
499 for j = 1:elemsize
500 k = j + elemsize*(i - 1)
501 @test 10 + k÷2 == ev[k]*(-1)^((i - 1)%2)
506 for i = 1:nelem
507 curlorients[1, i] = curlorients[6, i] = 0
508 if (i - 1)%2 > 0
509 # T = [0 -1]
510 # [-1 0]
512 curlorients[3, i] = -1
513 curlorients[4, i] = -1
517 curlorients[2, i] = 1
520 curlorients[5, i] = 1
527 1,
528 1,
529 nelem + 1,
536 for i = 1:nelem
537 for j = 1:elemsize
538 k = j + elemsize*(i - 1)
539 if (i - 1)%2 > 0
541 @test j == 1 || 10 + i - 1 == -ev[k]