15   PetscCall(PetscInitialize(&argc, &argv, NULL, help));  in main()
19   PetscCall(MatCreate(PETSC_COMM_WORLD, &mat));  in main()
20   PetscCall(MatSetSizes(mat, PETSC_DECIDE, PETSC_DECIDE, 4, 4));  in main()
21   PetscCall(MatSetFromOptions(mat));  in main()
22   PetscCall(MatSetUp(mat));  in main()
28   PetscCall(MatSetValues(mat, 1, &i, 1, &j, &v, INSERT_VALUES));  in main()
32   PetscCall(MatSetValues(mat, 1, &i, 1, &j, &v, INSERT_VALUES));  in main()
36   PetscCall(MatSetValues(mat, 1, &i, 1, &j, &v, INSERT_VALUES));  in main()
40   PetscCall(MatSetValues(mat, 1, &i, 1, &j, &v, INSERT_VALUES));  in main()
41   PetscCall(MatAssemblyBegin(mat, MAT_FINAL_ASSEMBLY));  in main()
42   PetscCall(MatAssemblyEnd(mat, MAT_FINAL_ASSEMBLY));  in main()
44   PetscCall(PetscViewerASCIIGetStdout(PETSC_COMM_WORLD, &viewer));  in main()
45   PetscCall(PetscViewerPushFormat(viewer, PETSC_VIEWER_ASCII_DENSE));  in main()
46   PetscCall(PetscViewerASCIIPrintf(viewer, "Original matrix\n"));  in main()
47   PetscCall(MatView(mat, viewer));  in main()
49   PetscCall(MatGetOrdering(mat, MATORDERINGNATURAL, &isrow, &iscol));  in main()
51   PetscCall(MatPermute(mat, isrow, iscol, &B));  in main()
52   PetscCall(PetscViewerASCIIPrintf(viewer, "Original matrix permuted by identity\n"));  in main()
53   PetscCall(MatView(B, viewer));  in main()
54   PetscCall(MatDestroy(&B));  in main()
56   PetscCall(MatReorderForNonzeroDiagonal(mat, 1.e-8, isrow, iscol));  in main()
57   PetscCall(MatPermute(mat, isrow, iscol, &B));  in main()
58 …PetscCall(PetscViewerASCIIPrintf(viewer, "Original matrix permuted by identity + NonzeroDiagonal()…  in main()
59   PetscCall(MatView(B, viewer));  in main()
60   PetscCall(PetscViewerASCIIPrintf(viewer, "Row permutation\n"));  in main()
61   PetscCall(ISView(isrow, viewer));  in main()
62   PetscCall(PetscViewerASCIIPrintf(viewer, "Column permutation\n"));  in main()
63   PetscCall(ISView(iscol, viewer));  in main()
64   PetscCall(MatDestroy(&B));  in main()
66   PetscCall(ISDestroy(&isrow));  in main()
67   PetscCall(ISDestroy(&iscol));  in main()
69   PetscCall(MatGetOrdering(mat, MATORDERINGND, &isrow, &iscol));  in main()
70   PetscCall(MatPermute(mat, isrow, iscol, &B));  in main()
71   PetscCall(PetscViewerASCIIPrintf(viewer, "Original matrix permuted by ND\n"));  in main()
72   PetscCall(MatView(B, viewer));  in main()
73   PetscCall(MatDestroy(&B));  in main()
74   PetscCall(PetscViewerASCIIPrintf(viewer, "ND row permutation\n"));  in main()
75   PetscCall(ISView(isrow, viewer));  in main()
76   PetscCall(PetscViewerASCIIPrintf(viewer, "ND column permutation\n"));  in main()
77   PetscCall(ISView(iscol, viewer));  in main()
79   PetscCall(MatReorderForNonzeroDiagonal(mat, 1.e-8, isrow, iscol));  in main()
80   PetscCall(MatPermute(mat, isrow, iscol, &B));  in main()
81   PetscCall(PetscViewerASCIIPrintf(viewer, "Original matrix permuted by ND + NonzeroDiagonal()\n"));  in main()
82   PetscCall(MatView(B, viewer));  in main()
83   PetscCall(MatDestroy(&B));  in main()
84   PetscCall(PetscViewerASCIIPrintf(viewer, "ND + NonzeroDiagonal() row permutation\n"));  in main()
85   PetscCall(ISView(isrow, viewer));  in main()
86   PetscCall(PetscViewerASCIIPrintf(viewer, "ND + NonzeroDiagonal() column permutation\n"));  in main()
87   PetscCall(ISView(iscol, viewer));  in main()
89   PetscCall(ISDestroy(&isrow));  in main()
90   PetscCall(ISDestroy(&iscol));  in main()
92   PetscCall(MatGetOrdering(mat, MATORDERINGRCM, &isrow, &iscol));  in main()
93   PetscCall(MatPermute(mat, isrow, iscol, &B));  in main()
94   PetscCall(PetscViewerASCIIPrintf(viewer, "Original matrix permuted by RCM\n"));  in main()
95   PetscCall(MatView(B, viewer));  in main()
96   PetscCall(MatDestroy(&B));  in main()
97   PetscCall(PetscViewerASCIIPrintf(viewer, "RCM row permutation\n"));  in main()
98   PetscCall(ISView(isrow, viewer));  in main()
99   PetscCall(PetscViewerASCIIPrintf(viewer, "RCM column permutation\n"));  in main()
100   PetscCall(ISView(iscol, viewer));  in main()
102   PetscCall(MatReorderForNonzeroDiagonal(mat, 1.e-8, isrow, iscol));  in main()
103   PetscCall(MatPermute(mat, isrow, iscol, &B));  in main()
104 …PetscCall(PetscViewerASCIIPrintf(viewer, "Original matrix permuted by RCM + NonzeroDiagonal()\n"));  in main()
105   PetscCall(MatView(B, viewer));  in main()
106   PetscCall(PetscViewerASCIIPrintf(viewer, "RCM + NonzeroDiagonal() row permutation\n"));  in main()
107   PetscCall(ISView(isrow, viewer));  in main()
108   PetscCall(PetscViewerASCIIPrintf(viewer, "RCM + NonzeroDiagonal() column permutation\n"));  in main()
109   PetscCall(ISView(iscol, viewer));  in main()
112     PetscCall(MatSetOption(B, MAT_SYMMETRIC, PETSC_TRUE));  in main()
113     PetscCall(ISCreateStride(PETSC_COMM_SELF, 4, 0, 1, &identity));  in main()
114     PetscCall(MatPermute(B, identity, identity, &C));  in main()
115     PetscCall(MatConvert(C, MATSEQSBAIJ, MAT_INPLACE_MATRIX, &C));  in main()
116     PetscCall(MatDestroy(&C));  in main()
117     PetscCall(ISDestroy(&identity));  in main()
119   PetscCall(MatDestroy(&B));  in main()
123     PetscCall(MatSetValues(mat, 1, &i, 1, &i, &v, INSERT_VALUES));  in main()
125   PetscCall(MatAssemblyBegin(mat, MAT_FINAL_ASSEMBLY));  in main()
126   PetscCall(MatAssemblyEnd(mat, MAT_FINAL_ASSEMBLY));  in main()
127   PetscCall(MatLUFactor(mat, isrow, iscol, NULL));  in main()
130   PetscCall(ISDestroy(&isrow));  in main()
131   PetscCall(ISDestroy(&iscol));  in main()
132   PetscCall(MatDestroy(&mat));  in main()
134   PetscCall(PetscFinalize());  in main()