ex26.c (be7a8307fd44d09b1b5ae337caf3795a085b8258) ex26.c (ed5e4e85792117e4bf875eb7efb50a9b8564ac26)
1static char help[] = "Test FEM layout with DM and ExodusII storage\n\n";
2
3/*
4 In order to see the vectors which are being tested, use
5
6 -ua_vec_view -s_vec_view
7*/
8

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

37 ierr = PetscOptionsBegin(PETSC_COMM_WORLD, NULL, "FEM Layout Options", "ex26");CHKERRQ(ierr);
38 ierr = PetscOptionsString("-i", "Filename to read", "ex26", ifilename, ifilename, sizeof(ifilename), NULL);CHKERRQ(ierr);
39 ierr = PetscOptionsString("-o", "Filename to write", "ex26", ofilename, ofilename, sizeof(ofilename), NULL);CHKERRQ(ierr);
40 ierr = PetscOptionsBoundedInt("-order", "FEM polynomial order", "ex26", order, &order, NULL,1);CHKERRQ(ierr);
41 ierr = PetscOptionsEnd();CHKERRQ(ierr);
42 if ((order > 2) || (order < 1)) SETERRQ1(PETSC_COMM_WORLD, PETSC_ERR_ARG_OUTOFRANGE, "Unsupported polynomial order %D not in [1, 2]", order);
43
44 /* Read the mesh from a file in any supported format */
1static char help[] = "Test FEM layout with DM and ExodusII storage\n\n";
2
3/*
4 In order to see the vectors which are being tested, use
5
6 -ua_vec_view -s_vec_view
7*/
8

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

37 ierr = PetscOptionsBegin(PETSC_COMM_WORLD, NULL, "FEM Layout Options", "ex26");CHKERRQ(ierr);
38 ierr = PetscOptionsString("-i", "Filename to read", "ex26", ifilename, ifilename, sizeof(ifilename), NULL);CHKERRQ(ierr);
39 ierr = PetscOptionsString("-o", "Filename to write", "ex26", ofilename, ofilename, sizeof(ofilename), NULL);CHKERRQ(ierr);
40 ierr = PetscOptionsBoundedInt("-order", "FEM polynomial order", "ex26", order, &order, NULL,1);CHKERRQ(ierr);
41 ierr = PetscOptionsEnd();CHKERRQ(ierr);
42 if ((order > 2) || (order < 1)) SETERRQ1(PETSC_COMM_WORLD, PETSC_ERR_ARG_OUTOFRANGE, "Unsupported polynomial order %D not in [1, 2]", order);
43
44 /* Read the mesh from a file in any supported format */
45 ierr = DMPlexCreateFromFile(PETSC_COMM_WORLD, ifilename, "ex26_plex", PETSC_TRUE, &dm);CHKERRQ(ierr);
45 ierr = DMPlexCreateFromFile(PETSC_COMM_WORLD, ifilename, NULL, PETSC_TRUE, &dm);CHKERRQ(ierr);
46 ierr = DMSetFromOptions(dm);CHKERRQ(ierr);
47 ierr = DMViewFromOptions(dm, NULL, "-dm_view");CHKERRQ(ierr);
48 ierr = DMGetDimension(dm, &sdim);CHKERRQ(ierr);
49
50 /* Create the exodus result file */
51 {
52 PetscInt numstep = 3, step;
53 char *nodalVarName[4];

--- 525 unchanged lines hidden --- 		46 ierr = DMSetFromOptions(dm);CHKERRQ(ierr);
47 ierr = DMViewFromOptions(dm, NULL, "-dm_view");CHKERRQ(ierr);
48 ierr = DMGetDimension(dm, &sdim);CHKERRQ(ierr);
49
50 /* Create the exodus result file */
51 {
52 PetscInt numstep = 3, step;
53 char *nodalVarName[4];

--- 525 unchanged lines hidden ---