libceedsetup.c (129d873643707da2348b9f89b9ef383c571b232a) libceedsetup.c (de1229c50b0fa287ef0d97965becc5ec55ca5a76)
1#include <stdio.h>
2#include "../include/libceedsetup.h"
3#include "../include/petscutils.h"
4
5// -----------------------------------------------------------------------------
6// Destroy libCEED operator objects
7// -----------------------------------------------------------------------------
8PetscErrorCode CeedDataDestroy(CeedInt i, CeedData data) {

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

50 CeedOperator op_setup_geo, op_apply;
51 CeedVector x_coord, q_data, x_ceed, y_ceed;
52 CeedInt num_qpts, c_start, c_end, num_elem,
53 q_data_size = bp_data.q_data_size;
54 CeedScalar R = 1, // radius of the sphere
55 l = 1.0/PetscSqrtReal(3.0); // half edge of the inscribed cube
56
57<<<<<<< HEAD
1#include <stdio.h>
2#include "../include/libceedsetup.h"
3#include "../include/petscutils.h"
4
5// -----------------------------------------------------------------------------
6// Destroy libCEED operator objects
7// -----------------------------------------------------------------------------
8PetscErrorCode CeedDataDestroy(CeedInt i, CeedData data) {

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

50 CeedOperator op_setup_geo, op_apply;
51 CeedVector x_coord, q_data, x_ceed, y_ceed;
52 CeedInt num_qpts, c_start, c_end, num_elem,
53 q_data_size = bp_data.q_data_size;
54 CeedScalar R = 1, // radius of the sphere
55 l = 1.0/PetscSqrtReal(3.0); // half edge of the inscribed cube
56
57<<<<<<< HEAD
58<<<<<<< HEAD
58 PetscFunctionBeginUser;
59 // CEED bases
60 P = degree + 1;
61 Q = P + q_extra;
62 CeedBasisCreateTensorH1Lagrange(ceed, topo_dim, num_comp_u, P, Q,
63 bp_data.q_mode,
64 &basis_u);
65 CeedBasisCreateTensorH1Lagrange(ceed, topo_dim, num_comp_x, 2, Q,
66 bp_data.q_mode,
67 &basis_x);
68 CeedBasisGetNumQuadraturePoints(basis_u, &num_qpts);
69
70 // CEED restrictions
71=======
72>>>>>>> 158419b6 (example/petsc: added CreateBasisFromPlex and tested with tensor basis)
73 ierr = DMSetCoordinateDim(dm, topo_dim); CHKERRQ(ierr);
59 PetscFunctionBeginUser;
60 // CEED bases
61 P = degree + 1;
62 Q = P + q_extra;
63 CeedBasisCreateTensorH1Lagrange(ceed, topo_dim, num_comp_u, P, Q,
64 bp_data.q_mode,
65 &basis_u);
66 CeedBasisCreateTensorH1Lagrange(ceed, topo_dim, num_comp_x, 2, Q,
67 bp_data.q_mode,
68 &basis_x);
69 CeedBasisGetNumQuadraturePoints(basis_u, &num_qpts);
70
71 // CEED restrictions
72=======
73>>>>>>> 158419b6 (example/petsc: added CreateBasisFromPlex and tested with tensor basis)
74 ierr = DMSetCoordinateDim(dm, topo_dim); CHKERRQ(ierr);
75=======
76 //ierr = DMSetCoordinateDim(dm, topo_dim); CHKERRQ(ierr);
77>>>>>>> 0fa86f50 (example/petsc: Added CreateDistributedDM in petscutils.c and some cleanup)
74 ierr = DMGetCoordinateDM(dm, &dm_coord); CHKERRQ(ierr);
75
76 // CEED bases
77 ierr = CreateBasisFromPlex(ceed, dm_coord, 0, 0, 0, 0, &basis_x); CHKERRQ(ierr);
78 ierr = CreateBasisFromPlex(ceed, dm, 0, 0, 0, 0, &basis_u); CHKERRQ(ierr);
79
80 // CEED restrictions
81 ierr = DMPlexSetClosurePermutationTensor(dm_coord, PETSC_DETERMINE, NULL);

--- 189 unchanged lines hidden --- 		78 ierr = DMGetCoordinateDM(dm, &dm_coord); CHKERRQ(ierr);
79
80 // CEED bases
81 ierr = CreateBasisFromPlex(ceed, dm_coord, 0, 0, 0, 0, &basis_x); CHKERRQ(ierr);
82 ierr = CreateBasisFromPlex(ceed, dm, 0, 0, 0, 0, &basis_u); CHKERRQ(ierr);
83
84 // CEED restrictions
85 ierr = DMPlexSetClosurePermutationTensor(dm_coord, PETSC_DETERMINE, NULL);

--- 189 unchanged lines hidden ---