| /petsc/src/dm/dt/dualspace/impls/lagrange/tutorials/ |
| H A D | ex2.c | 15 DM K; in main() local 28 …esh(PETSC_COMM_SELF, dim, PETSC_FALSE, cells, lo, hi, periodicity, PETSC_TRUE, 0, PETSC_TRUE, &K)); in main() 30 PetscCall(DMSetField(K, 0, NULL, (PetscObject)fe)); in main() 32 PetscCall(DMCreateDS(K)); in main() 35 PetscCall(DMCreateGlobalVector(K, &X)); in main() 36 PetscCall(DMProjectFunction(K, 0.0, initu, NULL, INSERT_ALL_VALUES, X)); in main() 37 PetscCall(DMViewFromOptions(K, NULL, "-dual_dm_view")); in main() 39 PetscCall(DMDestroy(&K)); in main()
|
| H A D | ex1.c | 11 DM K; in main() local 33 …CreateReferenceCell(PETSC_COMM_WORLD, DMPolytopeTypeSimpleShape(dim, (PetscBool)!tensorCell), &K)); in main() 34 PetscCall(PetscDualSpaceSetDM(dsp, K)); in main() 59 PetscCall(DMDestroy(&K)); in main()
|
| /petsc/src/ksp/ksp/impls/hpddm/cuda/ |
| H A D | hpddm.cu | 11 typedef HPDDM::downscaled_type<PetscScalar> K; in KSPSolve_HPDDM_CUDA_Private() typedef 14 typedef HPDDM::upscaled_type<PetscScalar> K; in KSPSolve_HPDDM_CUDA_Private() typedef 21 K *ptr, *host_ptr; in KSPSolve_HPDDM_CUDA_Private() 22 thrust::device_ptr<K> dptr[2]; in KSPSolve_HPDDM_CUDA_Private() 25 PetscCallCUDA(cudaMalloc((void **)&ptr, 2 * N * sizeof(K))); in KSPSolve_HPDDM_CUDA_Private() 30 PetscCallCUDA(cudaMemcpy(host_ptr, ptr, 2 * N * sizeof(K), cudaMemcpyDeviceToHost)); in KSPSolve_HPDDM_CUDA_Private() 32 PetscCallCUDA(cudaMemcpy(ptr + N, host_ptr + N, N * sizeof(K), cudaMemcpyHostToDevice)); in KSPSolve_HPDDM_CUDA_Private() 36 PetscCall(PetscLogGpuToCpu(2 * N * sizeof(K))); in KSPSolve_HPDDM_CUDA_Private() 37 PetscCall(PetscLogCpuToGpu(N * sizeof(K))); in KSPSolve_HPDDM_CUDA_Private()
|
| /petsc/share/petsc/datafiles/meshes/ |
| H A D | gmsh-3d-binary-32.msh | 32 GHIJKLMNOPQRS������?@������?@������?���… 54 …KOSK"L!!LIIP!#L"G… 56 …+_,G,_VHV`,ITYUJU^TKUYXLX… 60 …K}�-]9��-9y��m\k|�mkD|�SK…
|
| H A D | square_periodic_bin.msh | 9 $@�K���R�?$@hz���@$@4t���@$@����… 11 …4t���@hz���@ �K���R�?!… 14 …�����?Ht�牡] @`��v@I�|�X]@���@J�Vً|"@D�9J�@Kb 20 1��@�K�bK]@g�#@Ho�-'�@h�V~�T@t*㑐B 21 @iK����"@��mF�|@j��E�~��?SL�COu@k7 34 …K/J0/`=1… 36 …r:LKs?t@bS… 37 …7cC�Z7�3KL�…
|
| H A D | gmsh-3d-binary-64.msh | 30 GHIJKLMNOPQR… 48 …KOSK"M!!… 50 …,ITYVJVZTKUY… 54 …mcD{�SEK|�SK… 55 W_1}=KE|;D…
|
| H A D | mesh-3d-box-innersphere_bin.msh | 14 $��¿�aB����?1�K��(��?k+&� 17 …���,�? ��[���?�}A}�?I��ϓ����-HD��?9�hI�?�?Jn��ދ�����:�?�?�ZS�οK��_M'ڶ�q� ��[�?��t… 19 �g�?�-�)��?�(${9��?|���}��?�@S>�K�?���_/���?I�C`�v�?s>GPm%ӿ��v���)�?�xI!��?���۵�?��… 44 iK 88 ��~�=K;p��OD8… 91 …�lLN��K;pr�P… 92 …S(�)I=��=Kp�… 93 9B��K1M��O�D… 96 …Ic6�LlN.�KPl.… 97 …2_-D5y�KlM.…
|
| H A D | gmsh-3d-binary.msh2 | 43 …���?^������?I��?��?J������?A������?K`������?`�… 59 …KOSK"M!… 61 …I(V^UJ(U_VK(X^V… 65 …dmbD~�dSEK��dSK… 66 dT]1d=KE�d;D…
|
| H A D | square_bin_physnames.msh | 25 ��?���1�h�?/��1�h�?���m���?0H�5�K�?�Y�� 27 ���?KeQ8 44 …I Y=dJ X;eK DA�… 46 …Ki Lj Mk … 48 …K� O8x� Jr6� …
|
| H A D | surfacesphere_bin.msh | 13 �?֏o ��H�>�pGο����Q�?ޡ�|���I��[%H��.k݉��"=�Ċ�J_�h�ܿŧI���I��+�߿KĈtR���F����ؿ�0~… 19 …?~P}���Iݿ����.��?������?���=�鿴I��r��?�"w����?��S�.߿�o�b��?��K��?�$=�Ċ�?��[%H… 21 …}��̿u���X*�7��?����|�AKm���Կ�B�X��?�g`*3տl���Z�����:��?��䨟K��]��ƿ"�!�@�?�#�… 22 …�V4��?�����SϿ�����꿫���RE�?1`V����'�����ǿ�Y�{P�N�?���:�q���ſ�K ��0�?7d��R�ԿS8���… 26 …¿v�?��w��ѿ��U9���?�������?P}���Iݿ����.��?��"w����?���=�鿴I��r��?���K��?�S�.߿�o�b��? 35 …KDUIJEIVLF…
|
| /petsc/src/dm/dt/fe/impls/composite/ |
| H A D | fecomposite.c | 19 DM K; in PetscFESetUp_Composite() local 31 PetscCall(PetscDualSpaceGetDM(fem->dualSpace, &K)); in PetscFESetUp_Composite() 32 PetscCall(DMGetDimension(K, &dim)); in PetscFESetUp_Composite() 33 PetscCall(DMPlexGetCellType(K, 0, &ct)); in PetscFESetUp_Composite() 42 PetscCall(DMGetWorkArray(K, dim, MPIU_REAL, &subpoint)); in PetscFESetUp_Composite() 49 PetscCall(DMPlexGetTransitiveClosure(K, s, PETSC_TRUE, &closureSize, &closure)); in PetscFESetUp_Composite() 58 PetscCall(DMPlexRestoreTransitiveClosure(K, s, PETSC_TRUE, &closureSize, &closure)); in PetscFESetUp_Composite() 61 PetscCall(DMRestoreWorkArray(K, dim, MPIU_REAL, &subpoint)); in PetscFESetUp_Composite() 100 …n_Composite(PetscFE fem, PetscInt npoints, const PetscReal points[], PetscInt K, PetscTabulation T) in PetscFEComputeTabulation_Composite() argument 110 PetscReal *B = K >= 0 ? T->T[0] : NULL; in PetscFEComputeTabulation_Composite() [all …]
|
| /petsc/src/snes/tutorials/ |
| H A D | ex46.c | 14 PetscReal K; /* The advection coefficient */ member 44 user.K = 1.0; in main() 45 PetscCall(PetscOptionsReal("-K", "The advection coefficient K", __FILE__, user.K, &user.K, NULL)); in main() 125 PetscReal D, K, hx, hy, hxdhy, hydhx; in FormFunctionLocal() local 130 K = user->K; in FormFunctionLocal() 150 …f[j][i] = D * (uxx + uyy) - (K * funcA(x[j][i], user) * PetscSqrtScalar(ux * ux + uy * uy) + funcU… in FormFunctionLocal() 166 PetscScalar D, K, A, v[5], hx, hy, hxdhy, hydhx, ux, uy; in FormJacobianLocal() local 172 K = user->K; in FormJacobianLocal() 211 …v[2] = D * 2.0 * (hydhx + hxdhy) + K * (funcADer(x[j][i], user) * normGradZ - A / normGradZ) *… in FormJacobianLocal() 214 v[3] = -D * hydhx + K * A * hx * hy / (2.0 * normGradZ); in FormJacobianLocal() [all …]
|
| H A D | ex33.c | 37 PetscScalar *K; in FormPermeability() local 45 PetscCall(DMDAVecGetArray(da, Kappa, &K)); in FormPermeability() 49 K[i] = 1.0; in FormPermeability() 52 if (i == (xs + xm) / 2) K[i] = 0.00000001; in FormPermeability() 53 else K[i] = 1.0; in FormPermeability() 56 PetscCall(DMDAVecRestoreArray(da, Kappa, &K)); in FormPermeability() 90 PetscScalar K = 2 * dx / (dx / Kappa[i] + dx / Kappa[i - 1]); in FormFunctionLocal() local 98 f[i].v = u[i].v + K * lambda * (u[i].p - u[i - 1].p) / dx; in FormFunctionLocal()
|
| /petsc/src/mat/tests/ |
| H A D | ex248.c | 7 Mat A, B, C, K, Ad, Bd; in main() local 27 PetscCall(MatCreateKAIJ(A, p, q, NULL, Bv, &K)); in main() 29 PetscCall(MatMultEqual(C, K, 10, &flg)); in main() 36 PetscCall(MatKAIJSetT(K, p, q, Bv)); in main() 38 PetscCall(MatMultEqual(C, K, 10, &flg)); in main() 40 PetscCall(MatDestroy(&K)); in main()
|
| H A D | ex268.c | 5 …s(PetscInt sdim, PetscInt M, PetscInt N, const PetscInt *J, const PetscInt *K, PetscScalar *ptr, P… in GenEntries() argument 14 …; d++) diff += (coords[J[j] * sdim + d] - coords[K[k] * sdim + d]) * (coords[J[j] * sdim + d] - co… in GenEntries() 25 PetscInt m = 100, dim = 3, M, K = 10, begin, n = 0; in main() local 38 PetscCall(PetscOptionsGetInt(NULL, NULL, "-K", &K, NULL)); in main() 48 PetscCall(MatCreateDense(PETSC_COMM_WORLD, m, PETSC_DECIDE, M, K, NULL, &B)); in main() 61 PetscCall(MatCreateDense(PETSC_COMM_WORLD, m, PETSC_DECIDE, M, K, NULL, &X)); in main() 62 PetscCall(MatCreateDense(PETSC_COMM_WORLD, m, PETSC_DECIDE, M, K, NULL, &Xd)); in main()
|
| H A D | ex70.c | 206 PetscInt m, n, k, M = 10, N = 10, K = 5, ldx = 3, ldb = 5, ldr = 4; in main() local 220 PetscCall(PetscOptionsGetInt(NULL, NULL, "-K", &K, NULL)); in main() 286 PetscCall(PetscMalloc1((m + ldx) * K, &dataX)); in main() 287 PetscCall(PetscMalloc1((n + ldb) * K, &dataB)); in main() 288 for (i = 0; i < (m + ldx) * K; i++) dataX[i] = MAGIC_NUMBER; in main() 289 for (i = 0; i < (n + ldb) * K; i++) dataB[i] = MAGIC_NUMBER; in main() 291 PetscCall(MatCreateDense(PETSC_COMM_WORLD, n, PETSC_DECIDE, N, K, dataB, &B)); in main() 292 PetscCall(MatCreateDense(PETSC_COMM_WORLD, m, PETSC_DECIDE, M, K, dataX, &X)); in main() 304 PetscCall(MatCreateDense(PETSC_COMM_WORLD, k, n, K, N, dataBt, &Bt)); in main() 345 …PetscCall(MatDenseGetSubMatrix(B, PETSC_DECIDE, PETSC_DECIDE, PetscMin(1, K - 1), PetscMin(2, K), … in main() [all …]
|
| /petsc/src/binding/petsc4py/demo/legacy/perftest/ |
| H A D | App.f90 | 24 real(kind=8), intent(out) :: K local 26 K = 1.0 28 K = 1.0/(1+P**2) 34 real(kind=8), intent(out) :: K(-1:1) local 38 K(-1) = -Ka 39 K( 0) = Ka + Kb 40 K(+1) = -Kb
|
| /petsc/src/dm/dt/tests/ |
| H A D | ex17.c | 8 DM K; in main() local 14 PetscCall(DMPlexCreateReferenceCell(PETSC_COMM_SELF, DM_POLYTOPE_QUADRILATERAL, &K)); in main() 18 PetscCall(PetscDualSpaceSetDM(h_1, K)); in main() 27 PetscCall(PetscDualSpaceSetDM(h_div, K)); in main() 46 PetscCall(DMDestroy(&K)); in main()
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex69.c | 93 Mat K; in main() local 155 PetscCall(DMCreateMatrix(da, &K)); in main() 156 PetscCall(MatSetOption(K, MAT_NEW_NONZERO_ALLOCATION_ERR, PETSC_FALSE)); in main() 166 PetscCall(MatSetValues(K, n, rows, n, rows, &A[0][0], ADD_VALUES)); in main() 168 PetscCall(MatAssemblyBegin(K, MAT_FINAL_ASSEMBLY)); in main() 169 PetscCall(MatAssemblyEnd(K, MAT_FINAL_ASSEMBLY)); in main() 171 PetscCall(MatCreateVecs(K, &x, &b)); in main() 179 PetscCall(MatZeroRowsColumns(K, 2, rows, 1.0, x, b)); in main() 183 PetscCall(KSPSetOperators(ksp, K, K)); in main() 208 PetscCall(MatDestroy(&K)); in main()
|
| H A D | ex68.c | 56 Mat K; in main() local 88 PetscCall(MatCreateSeqDense(PETSC_COMM_SELF, n, n, &A[0][0], &K)); in main() 92 PetscCall(MatCreateVecs(K, &x, &b)); in main() 97 PetscCall(MatZeroRowsColumns(K, 2, rows, 1.0, x, b)); in main() 98 PetscCall(KSPSetOperators(ksp, K, K)); in main() 122 PetscCall(MatDestroy(&K)); in main()
|
| H A D | ex77f.F90 | 15 PetscInt m, K 24 K = 5 25 PetscCallA(PetscOptionsGetInt(PETSC_NULL_OPTIONS, PETSC_NULL_CHARACTER, '-n', K, flg, ierr)) 33 …PetscCallA(MatCreateDense(PETSC_COMM_WORLD, m, PETSC_DECIDE, PETSC_DECIDE, K, PETSC_NULL_SCALAR_AR… 34 …PetscCallA(MatCreateDense(PETSC_COMM_WORLD, m, PETSC_DECIDE, PETSC_DECIDE, K, PETSC_NULL_SCALAR_AR…
|
| /petsc/include/petsc/private/cpp/ |
| H A D | unordered_map.hpp | 1211 template <typename K, typename T, typename H = std::hash<K>, 1215 typename KE = std::equal_to<K> 1277 template <typename K, typename T, typename H, typename KE> 1279 inline std::pair<typename UnorderedMap<K, T, H, KE>::iterator, bool> UnorderedMap<K, T, H, KE>::emp… in emplace_() 1291 template <typename K, typename T, typename H, typename KE> 1293 inline std::pair<typename UnorderedMap<K, T, H, KE>::iterator, bool> UnorderedMap<K, T, H, KE>::emp… in emplace_() 1302 template <typename K, typename T, typename H, typename KE> 1303 inline typename UnorderedMap<K, T, H, KE>::iterator UnorderedMap<K, T, H, KE>::find(const key_type … in find() 1312 template <typename K, typename T, typename H, typename KE> 1313 inline typename UnorderedMap<K, T, H, KE>::const_iterator UnorderedMap<K, T, H, KE>::find(const key… in find() [all …]
|
| /petsc/src/dm/tests/ |
| H A D | noflux_check.c | 7 PetscErrorCode globalKMat_3d(Mat K, DMDALocalInfo info) in globalKMat_3d() argument 57 PetscCall(MatSetValuesStencil(K, 1, &row, ncols, col, vals, ADD_VALUES)); in globalKMat_3d() 61 PetscCall(MatAssemblyBegin(K, MAT_FINAL_ASSEMBLY)); in globalKMat_3d() 62 PetscCall(MatAssemblyEnd(K, MAT_FINAL_ASSEMBLY)); in globalKMat_3d() 66 PetscErrorCode globalKMat_2d(Mat K, DMDALocalInfo info) in globalKMat_2d() argument 99 PetscCall(MatSetValuesStencil(K, 1, &row, ncols, col, vals, ADD_VALUES)); in globalKMat_2d() 102 PetscCall(MatAssemblyBegin(K, MAT_FINAL_ASSEMBLY)); in globalKMat_2d() 103 PetscCall(MatAssemblyEnd(K, MAT_FINAL_ASSEMBLY)); in globalKMat_2d()
|
| /petsc/src/ksp/ksp/tests/ |
| H A D | lostnullspace.c | 14 Mat A, K; in main() local 92 PetscCall(KSPGetOperators(sub_ksp[1], &K, PETSC_NULLPTR)); in main() 93 PetscCall(MatSetNearNullSpace(K, near_null_space)); in main() 94 PetscCall(MatSetBlockSize(K, 3)); in main() 96 PetscCall(MatGetSize(K, &rows, &cols)); in main() 97 PetscCall(MatGetBlockSize(K, &bs)); in main() 98 PetscCall(MatGetNearNullSpace(K, &nullsp)); in main()
|
| H A D | ex10.c | 18 extern PetscErrorCode paulintegrate20(PetscReal K[60][60]); 88 PetscReal **K, norm; in GetElasticityMatrix() local 98 PetscCall(PetscMalloc1(81, &K)); in GetElasticityMatrix() 99 for (i = 0; i < 81; i++) PetscCall(PetscMalloc1(81, &K[i])); in GetElasticityMatrix() 100 PetscCall(Elastic20Stiff(K)); in GetElasticityMatrix() 120 PetscCall(AddElement(mat, r1, r2, K, h1, h2)); in GetElasticityMatrix() 133 for (i = 0; i < 81; i++) PetscCall(PetscFree(K[i])); in GetElasticityMatrix() 134 PetscCall(PetscFree(K)); in GetElasticityMatrix() 168 PetscErrorCode AddElement(Mat mat, PetscInt r1, PetscInt r2, PetscReal **K, PetscInt h1, PetscInt h… in AddElement() argument 179 if (K[h1 + l1][h2 + l2] != 0.0) { in AddElement() [all …]
|