mat-ceed.c (f674649b871b4939edeaefdaebac9e77c0983ca8) mat-ceed.c (9bc663991d6482bcb1d60b1f116148f11db83fa1)
1/// @file
2/// MatCEED implementation
3
4#include <ceed.h>
5#include <ceed/backend.h>
6#include <mat-ceed-impl.h>
7#include <mat-ceed.h>
8#include <petsc-ceed-utils.h>

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

1224 PetscCall(VecGetArrayReadAndMemType(X, &x, &(*ctx)->mem_type));
1225 PetscCall(VecRestoreArrayReadAndMemType(X, &x));
1226 PetscCall(DMRestoreLocalVector(dm_x, &X));
1227 }
1228
1229 // libCEED objects
1230 PetscCheck(CeedOperatorGetCeed(op_mult, &(*ctx)->ceed) == CEED_ERROR_SUCCESS, PETSC_COMM_SELF, PETSC_ERR_LIB,
1231 "retrieving Ceed context object failed");
1/// @file
2/// MatCEED implementation
3
4#include <ceed.h>
5#include <ceed/backend.h>
6#include <mat-ceed-impl.h>
7#include <mat-ceed.h>
8#include <petsc-ceed-utils.h>

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

1224 PetscCall(VecGetArrayReadAndMemType(X, &x, &(*ctx)->mem_type));
1225 PetscCall(VecRestoreArrayReadAndMemType(X, &x));
1226 PetscCall(DMRestoreLocalVector(dm_x, &X));
1227 }
1228
1229 // libCEED objects
1230 PetscCheck(CeedOperatorGetCeed(op_mult, &(*ctx)->ceed) == CEED_ERROR_SUCCESS, PETSC_COMM_SELF, PETSC_ERR_LIB,
1231 "retrieving Ceed context object failed");
1232 PetscCallCeed((*ctx)->ceed, CeedReference((*ctx)->ceed));
1233 PetscCallCeed((*ctx)->ceed, CeedOperatorGetActiveVectorLengths(op_mult, &x_loc_len, &y_loc_len));
1234 PetscCallCeed((*ctx)->ceed, CeedOperatorReferenceCopy(op_mult, &(*ctx)->op_mult));
1235 if (op_mult_transpose) PetscCallCeed((*ctx)->ceed, CeedOperatorReferenceCopy(op_mult_transpose, &(*ctx)->op_mult_transpose));
1236 PetscCallCeed((*ctx)->ceed, CeedVectorCreate((*ctx)->ceed, x_loc_len, &(*ctx)->x_loc));
1237 PetscCallCeed((*ctx)->ceed, CeedVectorCreate((*ctx)->ceed, y_loc_len, &(*ctx)->y_loc));
1238
1239 // Flop counting
1240 {

--- 282 unchanged lines hidden --- 		1232 PetscCallCeed((*ctx)->ceed, CeedOperatorGetActiveVectorLengths(op_mult, &x_loc_len, &y_loc_len));
1233 PetscCallCeed((*ctx)->ceed, CeedOperatorReferenceCopy(op_mult, &(*ctx)->op_mult));
1234 if (op_mult_transpose) PetscCallCeed((*ctx)->ceed, CeedOperatorReferenceCopy(op_mult_transpose, &(*ctx)->op_mult_transpose));
1235 PetscCallCeed((*ctx)->ceed, CeedVectorCreate((*ctx)->ceed, x_loc_len, &(*ctx)->x_loc));
1236 PetscCallCeed((*ctx)->ceed, CeedVectorCreate((*ctx)->ceed, y_loc_len, &(*ctx)->y_loc));
1237
1238 // Flop counting
1239 {

--- 282 unchanged lines hidden ---