| /petsc/include/petsc/private/kernels/ |
| H A D | blockmatmult.h | 5 #define PetscKernel_v_gets_A_times_w_1_exp(v, A, w, exp) \ argument
7 v[0] exp A[0] * w[0]; \
10 #define PetscKernel_v_gets_A_times_w_2_exp(v, A, w, exp) \ argument
12 v[0] exp A[0] * w[0] + A[2] * w[1]; \
13 v[1] exp A[1] * w[0] + A[3] * w[1]; \
16 #define PetscKernel_v_gets_A_times_w_3_exp(v, A, w, exp) \ argument
18 v[0] exp A[0] * w[0] + A[3] * w[1] + A[6] * w[2]; \
19 v[1] exp A[1] * w[0] + A[4] * w[1] + A[7] * w[2]; \
20 v[2] exp A[2] * w[0] + A[5] * w[1] + A[8] * w[2]; \
23 #define PetscKernel_v_gets_A_times_w_4_exp(v, A, w, exp) \ argument
[all …]
|
| H A D | blockinvert.h | 42 W - bs by 1 work array
44 …W, allowzeropivot, zeropivotdetected) ((PetscErrorCode)(PetscLINPACKgefa((A), (bs), (pivots), (all… argument
56 W - square bs by bs work array
59 #define PetscKernel_A_gets_A_times_B(bs, A, B, W) \ argument
64 PetscCall(PetscArraycpy((W), (A), (bs) * (bs))); \
67 …StaticAnalyzerForExpressionUnderstandingThatThisIsDangerousAndBugprone(PetscAssertPointer(W, 4)); \
68 …PetscCallBLAS("BLASgemm", BLASgemm_("N", "N", &(_bbs), &(_bbs), &(_bbs), &_one, (W), &(_bbs), (B),…
106 v = v + A^T w v_gets_v_plus_Atranspose_times_w
110 w - array of length bs
112 #define PetscKernel_v_gets_v_plus_Atranspose_times_w(bs, v, A, w) \ argument
[all …]
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | baij.h | 270 - W - bs*bs work array
276 …PetscErrorCode PetscKernel_A_gets_A_times_B_2(PetscScalar *A, const PetscScalar *B, PetscScalar *W) in PetscKernel_A_gets_A_times_B_2() argument
278 PetscCall(PetscArraycpy(W, A, 4)); in PetscKernel_A_gets_A_times_B_2()
279 A[0] = W[0] * B[0] + W[2] * B[1]; in PetscKernel_A_gets_A_times_B_2()
280 A[1] = W[1] * B[0] + W[3] * B[1]; in PetscKernel_A_gets_A_times_B_2()
281 A[2] = W[0] * B[2] + W[2] * B[3]; in PetscKernel_A_gets_A_times_B_2()
282 A[3] = W[1] * B[2] + W[3] * B[3]; in PetscKernel_A_gets_A_times_B_2()
310 - W - bs*bs work array
316 …PetscErrorCode PetscKernel_A_gets_A_times_B_3(PetscScalar *A, const PetscScalar *B, PetscScalar *W) in PetscKernel_A_gets_A_times_B_3() argument
318 PetscCall(PetscArraycpy(W, A, 9)); in PetscKernel_A_gets_A_times_B_3()
[all …]
|
| /petsc/src/dm/impls/plex/tests/output/ |
| H A D | ex18_9_seq_hdf5_simple_int_err.out | 6 vertex 14 w/ coordinates ( 0.000, 50000.000, 83333.333)
7 vertex 13 w/ coordinates ( 0.000, 50000.000, 75000.000)
9 vertex 20 w/ coordinates ( 25000.000, 50000.000, 83333.333)
10 vertex 14 w/ coordinates ( 0.000, 50000.000, 83333.333)
12 vertex 24 w/ coordinates ( 50000.000, 50000.000, 75000.000)
13 vertex 20 w/ coordinates ( 25000.000, 50000.000, 83333.333)
15 vertex 13 w/ coordinates ( 0.000, 50000.000, 75000.000)
16 vertex 24 w/ coordinates ( 50000.000, 50000.000, 75000.000)
19 vertex 30 w/ coordinates ( 75000.000, 50000.000, 83333.333)
20 vertex 24 w/ coordinates ( 50000.000, 50000.000, 75000.000)
[all …]
|
| /petsc/src/sys/classes/draw/tests/ |
| H A D | ex7.c | 9 int i, j, w, h; in main() local
27 PetscCall(PetscDrawGetWindowSize(draw, &w, &h)); in main()
28 PetscCall(PetscDrawSetCoordinates(draw, 0, 0, --w, --h)); in main()
32 PetscCall(PetscDrawLine(draw, w, 0, w, 0, r)); in main()
34 PetscCall(PetscDrawLine(draw, w, h, w, h, r)); in main()
36 PetscCall(PetscDrawLine(draw, 0 + 2, 0, w - 2, 0, k)); in main()
37 PetscCall(PetscDrawLine(draw, 0 + 2, h, w - 2, h, k)); in main()
39 PetscCall(PetscDrawLine(draw, w, 0 + 2, w, h - 2, k)); in main()
41 PetscCall(PetscDrawLine(draw, 0 + 2, h / 2, w - 2, h - 2, b)); in main()
42 PetscCall(PetscDrawLine(draw, 0 + 1, h - 1, w - 1, 0 + 1, b)); in main()
[all …]
|
| /petsc/share/petsc/datafiles/meshes/ |
| H A D | surfacesphere_bin.msh | 7 …��������������������������+^Do�?۰/7���?��������������w��?�,u#��?�������� �����u#��?֗�w��?�…
8 ����07���? NDo�?�������� ���۰/7��ӿ+^Do�?�������� ����,u#���w��?��������
9 …w�鿅�u#��?����������� NDo07���?�����������+^Do�۰/7��ӿ��������������w���,u#������������…
|
| H A D | gmsh-3d-binary-32.msh | 36 ���T���U���V���W���X���Y���Z���[���\���]���^���_���`���^������?�������?^������?������?�������?���…
44 ���n���o���p���q���r���s���t���u���v���w���x���y���z��� ������?@������?�������?������?������?����…
47 …>���?����>�?�>���?W���>�? �>���?W���>�?�������?�������?�������? ���>�? ���>�? ���>�?X���>…
56 …W���6���W���\������7������Z������8������Z���X���9���X���]������:������Z������;���V���Z����…
|
| H A D | gmsh-3d-binary-64.msh | 33 �������T�������U�������V�������W�������X�������Y�������Z�������[�������\�������]�������^�������_���…
39 �������n�������o�������p�������q�������r�������s�������t�������u�������v�������w�������x�������y���…
41 …X���>�?�>���?�������?�������?�������?X���>�? �>���?W���>�? ���>�? ���>�? ���>�? �>���?W�…
50 …W�������9�������W�������_��������������:��������������\��������������;�������X�������\��������…
54 …����h�������u�������b�������a�������g�������v�������a�������b�������i�������w�������c�������b�����…
|
| H A D | mesh-3d-box-innersphere_bin.msh | 9 (g�������w�ɿ ���UDo�?�Q��ɐ�5�/7�Ƴ�
11 (g�������w��?�����������������VUUUUUſ�����������������TUUUUU�?����������VUUUUUſ�������?��������…
17 …w�ҹ�}m�wLq�?F�����k*�����Z���q���ߗʿG����"E ���?=��Yи?/Đ��@˿H���#���,�? ��[���?� }A}�?I�����ϓ…
26 …W������������������������W���T������������������������R���W�������������������…
54 …���U������ ���p������������������V������ ��� ���p���������������W������ ������o����…
|
| H A D | gmsh-3d-binary.msh2 | 45 …���?�������?A������?U���^������?�������?^������?V���������?�������?������?W���`������?�������?`�…
51 …w����f�x���?�f�x���?�������?x����MUUU��?�������?�������?y����������?�MUUU��?�������?z����}�&#�?�}…
61 …W���Z���������������<������(������`���W������������=������(���%���_���$������������>�…
|
| /petsc/src/dm/dt/tests/output/ |
| H A D | ex7_1234.out | 7 w:
10 w(v): 0.440064
14 u wedge w:
17 u wedge w(x): -0.385677
20 star w:
22 star star w:
30 w:
33 w(v): -0.710721
37 u wedge w:
40 u wedge w(x): 0.145303
[all …]
|
| /petsc/src/vec/is/sf/impls/window/ |
| H A D | sfwindow.c | 72 PetscSF_Window *w = (PetscSF_Window *)sf->data; in PetscSFWindowGetDataTypes() local
79 for (link = w->link; link; link = link->next) { in PetscSFWindowGetDataTypes()
119 link->next = w->link; in PetscSFWindowGetDataTypes()
120 w->link = link; in PetscSFWindowGetDataTypes()
181 PetscSF_Window *w = (PetscSF_Window *)sf->data; in PetscSFWindowSetFlavorType_Window() local
184 w->flavor = flavor; in PetscSFWindowSetFlavorType_Window()
214 PetscSF_Window *w = (PetscSF_Window *)sf->data; in PetscSFWindowGetFlavorType_Window() local
217 *flavor = w->flavor; in PetscSFWindowGetFlavorType_Window()
248 PetscSF_Window *w = (PetscSF_Window *)sf->data; in PetscSFWindowSetSyncType_Window() local
251 w->sync = sync; in PetscSFWindowSetSyncType_Window()
[all …]
|
| /petsc/src/sys/classes/draw/impls/x/ |
| H A D | ximpl.h | 41 int x, y, w, h; /* Location and size window */ member
44 #define PetscDrawXiDrawable(w) ((w)->drw ? (w)->drw : (w)->win) argument
46 static inline void PetscDrawXiSetPixVal(PetscDraw_X *W, PetscDrawXiPixVal pix) in PetscDrawXiSetPixVal() argument
48 if (W->gc.cur_pix != pix) { in PetscDrawXiSetPixVal()
49 XSetForeground(W->disp, W->gc.set, pix); in PetscDrawXiSetPixVal()
50 W->gc.cur_pix = pix; in PetscDrawXiSetPixVal()
55 …#define PetscDrawXiValidColor(W, color) PetscCheck((color) >= 0 && (color) < PETSC_DRAW_MAXCOLOR, … argument
57 #define PetscDrawXiValidColor(W, color) \ argument
62 #define PetscDrawXiSetColor(W, color) \ argument
64 PetscDrawXiValidColor(W, color); \
[all …]
|
| H A D | xops.c | 10 #define XTRANS(draw, xwin, x) ((int)(((xwin)->w - 1) * ((draw)->port_xl + (((x - (draw)->coor_xl) *…
13 …oor_xl + (((PetscReal)(i)) * ((draw)->coor_xr - (draw)->coor_xl) / ((xwin)->w - 1) - (draw)->port_…
19 int xa, ya, xb, yb, xmax = XiWin->w - 1, ymax = XiWin->h - 1; in PetscDrawSetViewport_X()
134 int x, y, w, h, c = (c1 + c2 + c3 + c4) / 4; in PetscDrawRectangle_X() local
139 w = XTRANS(draw, XiWin, xr) + 1 - x; in PetscDrawRectangle_X()
140 if (w <= 0) w = 1; in PetscDrawRectangle_X()
144 XFillRectangle(XiWin->disp, PetscDrawXiDrawable(XiWin), XiWin->gc.set, x, y, w, h); in PetscDrawRectangle_X()
151 int xA, yA, w, h; in PetscDrawEllipse_X() local
156 w = XTRANS(draw, XiWin, x + a / 2) + 1 - xA; in PetscDrawEllipse_X()
157 w = PetscAbs(w); in PetscDrawEllipse_X()
[all …]
|
| H A D | xtext.c | 19 PetscErrorCode PetscDrawXiFontFixed(PetscDraw_X *XBWin, int w, int h, PetscDrawXiFont **outfont) in PetscDrawXiFontFixed() argument
26 PetscCall(PetscDrawXiMatchFontSize(font, w, h)); in PetscDrawXiFontFixed()
37 int w, h, descent; member
84 nfonts[j].w = info[i].max_bounds.width; in PetscDrawXiInitFonts()
87 if (nfonts[j].w <= 0 || nfonts[j].h <= 0) continue; in PetscDrawXiInitFonts()
105 nfonts[j].w = info[i].max_bounds.width; in PetscDrawXiInitFonts()
106 /* nfonts[j].w = info[i].max_bounds.lbearing + info[i].max_bounds.rbearing; */ in PetscDrawXiInitFonts()
109 if (nfonts[j].w <= 0 || nfonts[j].h <= 0) continue; in PetscDrawXiInitFonts()
119 static PetscErrorCode PetscDrawXiMatchFontSize(PetscDrawXiFont *font, int w, int h) in PetscDrawXiMatchFontSize() argument
125 if (nfonts[i].w == w && nfonts[i].h == h) { in PetscDrawXiMatchFontSize()
[all …]
|
| /petsc/src/ksp/ksp/impls/cr/pipecr/ |
| H A D | pipecr.c | 25 Vec X, B, Z, P, W, Q, U, M, N; in KSPSolve_PIPECR() local
39 W = ksp->work[4]; in KSPSolve_PIPECR()
46 /* we don't have an R vector, so put the (unpreconditioned) residual in w for now */ in KSPSolve_PIPECR()
48 PetscCall(KSP_MatMult(ksp, Amat, X, W)); /* w <- b - Ax */ in KSPSolve_PIPECR()
49 PetscCall(VecAYPX(W, -1.0, B)); in KSPSolve_PIPECR()
51 PetscCall(VecCopy(B, W)); /* w <- b (x is 0) */ in KSPSolve_PIPECR()
53 PetscCall(KSP_PCApply(ksp, W, U)); /* u <- Bw */ in KSPSolve_PIPECR()
59 PetscCall(KSP_MatMult(ksp, Amat, U, W)); /* w <- Au */ in KSPSolve_PIPECR()
63 PetscCall(KSP_MatMult(ksp, Amat, U, W)); in KSPSolve_PIPECR()
77 PetscCall(KSP_PCApply(ksp, W, M)); /* m <- Bw */ in KSPSolve_PIPECR()
[all …]
|
| /petsc/lib/petsc/bin/maint/generators/ |
| H A D | java.py | 41 outfile = open('src/java/petsc.java','w')
120 w = f.read()
122 w = w.replace("#include <jni.h>",'#include "'+i+'.h"')
123 w = w.replace("#ifndef _Included_","#ifndef _string_that_does_not_exist_")
124 w = w.replace("#define _Included_","#define _another_string_that_does_not_exist_")
125 w = w.replace("#include <jni.h>","#include <JavaVM/jni.h>")
127 … t = re.search('JNIEXPORT void JNICALL Java_'+i+'[_0-9]*'+j+'\n [ ]* \(JNIEnv \*([, A-Za-z]*)\)',w)
142 …w = re.sub('JNIEXPORT void JNICALL Java_'+i+'([_0-9]*)'+j+'\n [ ]* \(JNIEnv \*, ([a-z]*)\)','JNIEX…
144 …w = re.sub('JNIEXPORT void JNICALL Java_'+i+'([_0-9]*)'+j+'\n [ ]* \(JNIEnv \*([, A-Za-z]*)\)','JN…
146 …w = re.sub('JNIEXPORT jint JNICALL Java_'+i+'([_0-9]*)'+i+'Create\n [ ]* \(JNIEnv \*, ([a-z]*)\)',…
[all …]
|
| /petsc/src/ksp/pc/impls/deflation/ |
| H A D | deflation.c | 198 static PetscErrorCode PCDeflationSetSpace_Deflation(PC pc, Mat W, PetscBool transpose) in PCDeflationSetSpace_Deflation() argument
203 /* possibly allows W' = Wt (which is valid but not tested) */ in PCDeflationSetSpace_Deflation()
204 PetscCall(PetscObjectReference((PetscObject)W)); in PCDeflationSetSpace_Deflation()
207 def->Wt = W; in PCDeflationSetSpace_Deflation()
209 PetscCall(MatDestroy(&def->W)); in PCDeflationSetSpace_Deflation()
210 def->W = W; in PCDeflationSetSpace_Deflation()
222 . W - deflation matrix
223 - transpose - indicates that W is an explicit transpose of the deflation matrix
228 Setting W as a multipliplicative `MATCOMPOSITE` enables use of the multilevel
229 deflation. If W = W0*W1*W2*...*Wn, W0 is taken as the first deflation space and
[all …]
|
| /petsc/src/ksp/ksp/impls/bcgs/pipebcgs/ |
| H A D | pipebcgs.c | 19 Vec X, B, S, R, RP, Y, Q, P2, Q2, R2, S2, W, Z, W2, Z2, T, V; in KSPSolve_PIPEBCGS() local
36 W = ksp->work[9]; in KSPSolve_PIPEBCGS()
78 PetscCall(KSP_MatMult(ksp, pc->mat, R2, W)); /* w <- A r2 */ in KSPSolve_PIPEBCGS()
81 PetscCall(VecDotBegin(W, RP, &d2)); /* d2 <- (w,rp) */ in KSPSolve_PIPEBCGS()
82 PetscCall(PetscCommSplitReductionBegin(PetscObjectComm((PetscObject)W))); in KSPSolve_PIPEBCGS()
83 PetscCall(KSP_PCApply(ksp, W, W2)); /* w2 <- K w */ in KSPSolve_PIPEBCGS()
85 PetscCall(VecDotEnd(W, RP, &d2)); in KSPSolve_PIPEBCGS()
95 PetscCall(VecCopy(W, S)); /* s <- w */ in KSPSolve_PIPEBCGS()
100 …PetscCall(VecAXPBYPCZ(S, 1.0, -beta * omega, beta, W, Z)); /* s <- beta * s + w - beta * ome… in KSPSolve_PIPEBCGS()
106 PetscCall(VecWAXPY(Y, -alpha, Z, W)); /* y <- w - alpha z */ in KSPSolve_PIPEBCGS()
[all …]
|
| /petsc/src/sys/classes/draw/impls/image/ |
| H A D | drawimage.c | 13 #define XTRANS(draw, img, x) ((int)(((img)->w - 1) * ((draw)->port_xl + ((((x) - (draw)->coor_xl) *…
16 …coor_xl + (((PetscReal)(i)) * ((draw)->coor_xr - (draw)->coor_xl) / ((img)->w - 1) - (draw)->port_…
25 int xmax = img->w - 1, ymax = img->h - 1; in PetscDrawSetViewport_Image()
111 *width = lw * (draw->coor_xr - draw->coor_xl) / (img->w * (draw->port_xr - draw->port_xl)); in PetscDrawLineGetWidth_Image()
172 int x = XTRANS(draw, img, xl), w = XTRANS(draw, img, xr) + 1 - x; in PetscDrawRectangle_Image() local
175 PetscImageDrawRectangle(img, x, y, w, h, c); in PetscDrawRectangle_Image()
189 int xc = XTRANS(draw, img, x), w = XTRANS(draw, img, x + a / 2) + 0 - xc; in PetscDrawEllipse_Image() local
191 if (PetscAbsReal(a - b) <= 0) w = h = PetscMin(w, h); /* workaround truncation errors */ in PetscDrawEllipse_Image()
192 PetscImageDrawEllipse(img, xc, yc, w, h, c); in PetscDrawEllipse_Image()
214 static PetscErrorCode PetscDrawStringSetSize_Image(PetscDraw draw,PetscReal w,PetscReal h)
[all …]
|
| /petsc/src/dm/dt/interface/ |
| H A D | dtaltv.c | 11 …A k-form w (k is called the "form degree" of w) is an alternating k-linear map acting on tuples (v…
13 …y two vectors in a tuple reverses the sign of the result, e.g. w(v_1, v_2, ..., v_k) = -w(v_2, v_1…
14 …- k-linear: w acts linear in each vector separately, e.g. w(a*v + b*y, v_2, ..., v_k) = a*w(v,v_2,…
36 …If k > 0, a k-form w and a vector v can combine to make a (k-1)-form through the interior product,…
37 defined by (w int v)(v_1,...,v_{k-1}) = w(v,v_1,...,v_{k-1}).
39 …The interior product is implemented for either two inputs (w and v) in PetscDTAltVInterior, for on…
40 matrix to multiply against multiple choices of w) in `PetscDTAltVInteriorMatrix()`,
43 …When there is a linear map L: V -> W from an N dimensional vector space to an M dimensional vector…
44 …it induces the linear pullback map L^* : Alt^k W -> Alt^k V, defined by L^* w(v_1,...,v_k) = w(L v…
45 …The pullback is implemented as `PetscDTAltVPullback()` (acting on a known w) and `PetscDTAltVPullb…
[all …]
|
| /petsc/src/snes/linesearch/impls/basic/ |
| H A D | linesearchbasic.c | 7 Vec X, F, Y, W; in SNESLineSearchApply_Basic() local
12 PetscCall(SNESLineSearchGetVecs(linesearch, &X, &F, &Y, &W, NULL)); in SNESLineSearchApply_Basic()
22 PetscCall(VecWAXPY(W, -lambda, Y, X)); in SNESLineSearchApply_Basic()
23 if (linesearch->ops->viproject) PetscCall((*linesearch->ops->viproject)(snes, W)); in SNESLineSearchApply_Basic()
26 PetscCall(SNESLineSearchPostCheck(linesearch, X, Y, W, &changed_y, &changed_w)); in SNESLineSearchApply_Basic()
28 if (!changed_w) PetscCall(VecWAXPY(W, -lambda, Y, X)); in SNESLineSearchApply_Basic()
29 if (linesearch->ops->viproject) PetscCall((*linesearch->ops->viproject)(snes, W)); in SNESLineSearchApply_Basic()
32 PetscCall((*linesearch->ops->snesfunc)(snes, W, F)); in SNESLineSearchApply_Basic()
38 PetscCall(VecNormBegin(W, NORM_2, &linesearch->xnorm)); in SNESLineSearchApply_Basic()
40 PetscCall(VecNormEnd(W, NORM_2, &linesearch->xnorm)); in SNESLineSearchApply_Basic()
[all …]
|
| /petsc/src/dm/tests/output/ |
| H A D | ex34_1.out | 3 Processor [0] M 11 N 5 P 7 m 3 n 2 p 2 w 2 s 1
5 Processor [1] M 11 N 5 P 7 m 3 n 2 p 2 w 2 s 1
7 Processor [2] M 11 N 5 P 7 m 3 n 2 p 2 w 2 s 1
9 Processor [3] M 11 N 5 P 7 m 3 n 2 p 2 w 2 s 1
11 Processor [4] M 11 N 5 P 7 m 3 n 2 p 2 w 2 s 1
13 Processor [5] M 11 N 5 P 7 m 3 n 2 p 2 w 2 s 1
15 Processor [6] M 11 N 5 P 7 m 3 n 2 p 2 w 2 s 1
17 Processor [7] M 11 N 5 P 7 m 3 n 2 p 2 w 2 s 1
19 Processor [8] M 11 N 5 P 7 m 3 n 2 p 2 w 2 s 1
21 Processor [9] M 11 N 5 P 7 m 3 n 2 p 2 w 2 s 1
[all …]
|
| /petsc/src/dm/tutorials/output/ |
| H A D | ex19_1.out | 4 Processor [0] M 8 m 1 w 1 s 1
9 Processor [0] M 8 m 1 w 2 s 1
14 Processor [0] M 8 m 1 w 2 s 1
19 Processor [0] M 8 N 8 m 1 n 1 w 2 s 1
24 Processor [0] M 8 N 8 m 1 n 1 w 2 s 1
29 Processor [0] M 8 N 8 m 1 n 1 w 2 s 1
34 Processor [0] M 8 N 8 P 8 m 1 n 1 p 1 w 2 s 1
39 Processor [0] M 8 N 8 P 8 m 1 n 1 p 1 w 2 s 1
44 Processor [0] M 8 N 8 P 8 m 1 n 1 p 1 w 2 s 1
|
| /petsc/src/dm/impls/patch/tests/output/ |
| H A D | ex1_3.out | 3 Processor [0] M 4 N 4 m 2 n 2 w 1 s 1
5 Processor [1] M 4 N 4 m 2 n 2 w 1 s 1
7 Processor [2] M 4 N 4 m 2 n 2 w 1 s 1
9 Processor [3] M 4 N 4 m 2 n 2 w 1 s 1
12 Processor [0] M 3 N 3 m 2 n 2 w 1 s 0
14 Processor [1] M 3 N 3 m 2 n 2 w 1 s 0
16 Processor [2] M 3 N 3 m 2 n 2 w 1 s 0
18 Processor [3] M 3 N 3 m 2 n 2 w 1 s 0
46 Processor [0] M 3 N 3 m 2 n 2 w 1 s 0
48 Processor [1] M 3 N 3 m 2 n 2 w 1 s 0
[all …]
|