ceed-elemrestriction.c (63fe681b35f0e91d70b9739fcc4dce8f46b867ae) ceed-elemrestriction.c (db002c03923317a1c3814dcd861330002c00a8ea)
1// Copyright (c) 2017-2022, Lawrence Livermore National Security, LLC and other CEED contributors.
2// All Rights Reserved. See the top-level LICENSE and NOTICE files for details.
3//
4// SPDX-License-Identifier: BSD-2-Clause
5//
6// This file is part of CEED: http://github.com/ceed
7
8#include <ceed-impl.h>

--- 346 unchanged lines hidden (view full) ---

355 return CEED_ERROR_SUCCESS;
356 }
357
358 CeedCheck(elem_size > 0, ceed, CEED_ERROR_DIMENSION, "Element size must be at least 1");
359 CeedCheck(num_comp > 0, ceed, CEED_ERROR_DIMENSION, "ElemRestriction must have at least 1 component");
360 CeedCheck(num_comp == 1 || comp_stride > 0, ceed, CEED_ERROR_DIMENSION, "ElemRestriction component stride must be at least 1");
361
362 CeedCall(CeedCalloc(1, rstr));
1// Copyright (c) 2017-2022, Lawrence Livermore National Security, LLC and other CEED contributors.
2// All Rights Reserved. See the top-level LICENSE and NOTICE files for details.
3//
4// SPDX-License-Identifier: BSD-2-Clause
5//
6// This file is part of CEED: http://github.com/ceed
7
8#include <ceed-impl.h>

--- 346 unchanged lines hidden (view full) ---

355 return CEED_ERROR_SUCCESS;
356 }
357
358 CeedCheck(elem_size > 0, ceed, CEED_ERROR_DIMENSION, "Element size must be at least 1");
359 CeedCheck(num_comp > 0, ceed, CEED_ERROR_DIMENSION, "ElemRestriction must have at least 1 component");
360 CeedCheck(num_comp == 1 || comp_stride > 0, ceed, CEED_ERROR_DIMENSION, "ElemRestriction component stride must be at least 1");
361
362 CeedCall(CeedCalloc(1, rstr));
363 (*rstr)->ceed = ceed;
364 CeedCall(CeedReference(ceed));
363 CeedCall(CeedReferenceCopy(ceed, &(*rstr)->ceed));
365 (*rstr)->ref_count = 1;
366 (*rstr)->num_elem = num_elem;
367 (*rstr)->elem_size = elem_size;
368 (*rstr)->num_comp = num_comp;
369 (*rstr)->comp_stride = comp_stride;
370 (*rstr)->l_size = l_size;
371 (*rstr)->num_blk = num_elem;
372 (*rstr)->blk_size = 1;

--- 38 unchanged lines hidden (view full) ---

411 return CEED_ERROR_SUCCESS;
412 }
413
414 CeedCheck(elem_size > 0, ceed, CEED_ERROR_DIMENSION, "Element size must be at least 1");
415 CeedCheck(num_comp > 0, ceed, CEED_ERROR_DIMENSION, "ElemRestriction must have at least 1 component");
416 CeedCheck(num_comp == 1 || comp_stride > 0, ceed, CEED_ERROR_DIMENSION, "ElemRestriction component stride must be at least 1");
417
418 CeedCall(CeedCalloc(1, rstr));
364 (*rstr)->ref_count = 1;
365 (*rstr)->num_elem = num_elem;
366 (*rstr)->elem_size = elem_size;
367 (*rstr)->num_comp = num_comp;
368 (*rstr)->comp_stride = comp_stride;
369 (*rstr)->l_size = l_size;
370 (*rstr)->num_blk = num_elem;
371 (*rstr)->blk_size = 1;

--- 38 unchanged lines hidden (view full) ---

410 return CEED_ERROR_SUCCESS;
411 }
412
413 CeedCheck(elem_size > 0, ceed, CEED_ERROR_DIMENSION, "Element size must be at least 1");
414 CeedCheck(num_comp > 0, ceed, CEED_ERROR_DIMENSION, "ElemRestriction must have at least 1 component");
415 CeedCheck(num_comp == 1 || comp_stride > 0, ceed, CEED_ERROR_DIMENSION, "ElemRestriction component stride must be at least 1");
416
417 CeedCall(CeedCalloc(1, rstr));
419 (*rstr)->ceed = ceed;
420 CeedCall(CeedReference(ceed));
418 CeedCall(CeedReferenceCopy(ceed, &(*rstr)->ceed));
421 (*rstr)->ref_count = 1;
422 (*rstr)->num_elem = num_elem;
423 (*rstr)->elem_size = elem_size;
424 (*rstr)->num_comp = num_comp;
425 (*rstr)->comp_stride = comp_stride;
426 (*rstr)->l_size = l_size;
427 (*rstr)->num_blk = num_elem;
428 (*rstr)->blk_size = 1;

--- 30 unchanged lines hidden (view full) ---

459 CeedCall(CeedElemRestrictionCreateStrided(delegate, num_elem, elem_size, num_comp, l_size, strides, rstr));
460 return CEED_ERROR_SUCCESS;
461 }
462
463 CeedCheck(elem_size > 0, ceed, CEED_ERROR_DIMENSION, "Element size must be at least 1");
464 CeedCheck(num_comp > 0, ceed, CEED_ERROR_DIMENSION, "ElemRestriction must have at least 1 component");
465
466 CeedCall(CeedCalloc(1, rstr));
419 (*rstr)->ref_count = 1;
420 (*rstr)->num_elem = num_elem;
421 (*rstr)->elem_size = elem_size;
422 (*rstr)->num_comp = num_comp;
423 (*rstr)->comp_stride = comp_stride;
424 (*rstr)->l_size = l_size;
425 (*rstr)->num_blk = num_elem;
426 (*rstr)->blk_size = 1;

--- 30 unchanged lines hidden (view full) ---

457 CeedCall(CeedElemRestrictionCreateStrided(delegate, num_elem, elem_size, num_comp, l_size, strides, rstr));
458 return CEED_ERROR_SUCCESS;
459 }
460
461 CeedCheck(elem_size > 0, ceed, CEED_ERROR_DIMENSION, "Element size must be at least 1");
462 CeedCheck(num_comp > 0, ceed, CEED_ERROR_DIMENSION, "ElemRestriction must have at least 1 component");
463
464 CeedCall(CeedCalloc(1, rstr));
467 (*rstr)->ceed = ceed;
468 CeedCall(CeedReference(ceed));
465 CeedCall(CeedReferenceCopy(ceed, &(*rstr)->ceed));
469 (*rstr)->ref_count = 1;
470 (*rstr)->num_elem = num_elem;
471 (*rstr)->elem_size = elem_size;
472 (*rstr)->num_comp = num_comp;
473 (*rstr)->l_size = l_size;
474 (*rstr)->num_blk = num_elem;
475 (*rstr)->blk_size = 1;
476 (*rstr)->is_oriented = false;

--- 43 unchanged lines hidden (view full) ---

520 return CEED_ERROR_SUCCESS;
521 }
522
523 CeedCheck(elem_size > 0, ceed, CEED_ERROR_DIMENSION, "Element size must be at least 1");
524 CeedCheck(blk_size > 0, ceed, CEED_ERROR_DIMENSION, "Block size must be at least 1");
525 CeedCheck(num_comp > 0, ceed, CEED_ERROR_DIMENSION, "ElemRestriction must have at least 1 component");
526 CeedCheck(num_comp == 1 || comp_stride > 0, ceed, CEED_ERROR_DIMENSION, "ElemRestriction component stride must be at least 1");
527
466 (*rstr)->ref_count = 1;
467 (*rstr)->num_elem = num_elem;
468 (*rstr)->elem_size = elem_size;
469 (*rstr)->num_comp = num_comp;
470 (*rstr)->l_size = l_size;
471 (*rstr)->num_blk = num_elem;
472 (*rstr)->blk_size = 1;
473 (*rstr)->is_oriented = false;

--- 43 unchanged lines hidden (view full) ---

517 return CEED_ERROR_SUCCESS;
518 }
519
520 CeedCheck(elem_size > 0, ceed, CEED_ERROR_DIMENSION, "Element size must be at least 1");
521 CeedCheck(blk_size > 0, ceed, CEED_ERROR_DIMENSION, "Block size must be at least 1");
522 CeedCheck(num_comp > 0, ceed, CEED_ERROR_DIMENSION, "ElemRestriction must have at least 1 component");
523 CeedCheck(num_comp == 1 || comp_stride > 0, ceed, CEED_ERROR_DIMENSION, "ElemRestriction component stride must be at least 1");
524
528 CeedCall(CeedCalloc(1, rstr));
529
530 CeedCall(CeedCalloc(num_blk * blk_size * elem_size, &blk_offsets));
531 CeedCall(CeedPermutePadOffsets(offsets, blk_offsets, num_blk, num_elem, blk_size, elem_size));
532
525 CeedCall(CeedCalloc(num_blk * blk_size * elem_size, &blk_offsets));
526 CeedCall(CeedPermutePadOffsets(offsets, blk_offsets, num_blk, num_elem, blk_size, elem_size));
527
533 (*rstr)->ceed = ceed;
534 CeedCall(CeedReference(ceed));
528 CeedCall(CeedCalloc(1, rstr));
529 CeedCall(CeedReferenceCopy(ceed, &(*rstr)->ceed));
535 (*rstr)->ref_count = 1;
536 (*rstr)->num_elem = num_elem;
537 (*rstr)->elem_size = elem_size;
538 (*rstr)->num_comp = num_comp;
539 (*rstr)->comp_stride = comp_stride;
540 (*rstr)->l_size = l_size;
541 (*rstr)->num_blk = num_blk;
542 (*rstr)->blk_size = blk_size;

--- 37 unchanged lines hidden (view full) ---

580 return CEED_ERROR_SUCCESS;
581 }
582
583 CeedCheck(elem_size > 0, ceed, CEED_ERROR_DIMENSION, "Element size must be at least 1");
584 CeedCheck(blk_size > 0, ceed, CEED_ERROR_DIMENSION, "Block size must be at least 1");
585 CeedCheck(num_comp > 0, ceed, CEED_ERROR_DIMENSION, "ElemRestriction must have at least 1 component");
586
587 CeedCall(CeedCalloc(1, rstr));
530 (*rstr)->ref_count = 1;
531 (*rstr)->num_elem = num_elem;
532 (*rstr)->elem_size = elem_size;
533 (*rstr)->num_comp = num_comp;
534 (*rstr)->comp_stride = comp_stride;
535 (*rstr)->l_size = l_size;
536 (*rstr)->num_blk = num_blk;
537 (*rstr)->blk_size = blk_size;

--- 37 unchanged lines hidden (view full) ---

575 return CEED_ERROR_SUCCESS;
576 }
577
578 CeedCheck(elem_size > 0, ceed, CEED_ERROR_DIMENSION, "Element size must be at least 1");
579 CeedCheck(blk_size > 0, ceed, CEED_ERROR_DIMENSION, "Block size must be at least 1");
580 CeedCheck(num_comp > 0, ceed, CEED_ERROR_DIMENSION, "ElemRestriction must have at least 1 component");
581
582 CeedCall(CeedCalloc(1, rstr));
588
589 (*rstr)->ceed = ceed;
590 CeedCall(CeedReference(ceed));
583 CeedCall(CeedReferenceCopy(ceed, &(*rstr)->ceed));
591 (*rstr)->ref_count = 1;
592 (*rstr)->num_elem = num_elem;
593 (*rstr)->elem_size = elem_size;
594 (*rstr)->num_comp = num_comp;
595 (*rstr)->l_size = l_size;
596 (*rstr)->num_blk = num_blk;
597 (*rstr)->blk_size = blk_size;
598 (*rstr)->is_oriented = false;

--- 352 unchanged lines hidden --- 		584 (*rstr)->ref_count = 1;
585 (*rstr)->num_elem = num_elem;
586 (*rstr)->elem_size = elem_size;
587 (*rstr)->num_comp = num_comp;
588 (*rstr)->l_size = l_size;
589 (*rstr)->num_blk = num_blk;
590 (*rstr)->blk_size = blk_size;
591 (*rstr)->is_oriented = false;

--- 352 unchanged lines hidden ---