Home
last modified time | relevance | path

Searched refs:w (Results 1 – 25 of 318) sorted by relevance

12345678910>>...13

/petsc/include/petsc/private/kernels/
H A Dblockmatmult.h5 #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 Dblockinvert.h112 #define PetscKernel_v_gets_v_plus_Atranspose_times_w(bs, v, A, w) \ argument
119 …StaticAnalyzerForExpressionUnderstandingThatThisIsDangerousAndBugprone(PetscAssertPointer(w, 4)); \
120 …PetscCallBLAS("BLASgemv", BLASgemv_("T", &(_bbs), &(_bbs), &_one, A, &(_bbs), w, &_ione, &_one, v,…
130 #define PetscKernel_v_gets_v_minus_A_times_w(bs, v, A, w) \ argument
137 …StaticAnalyzerForExpressionUnderstandingThatThisIsDangerousAndBugprone(PetscAssertPointer(w, 4)); \
138 …PetscCallBLAS("BLASgemv", BLASgemv_("N", &(_bbs), &(_bbs), &_mone, A, &(_bbs), w, &_ione, &_one, v…
148 #define PetscKernel_v_gets_v_minus_transA_times_w(bs, v, A, w) \ argument
155 …StaticAnalyzerForExpressionUnderstandingThatThisIsDangerousAndBugprone(PetscAssertPointer(w, 4)); \
156 …PetscCallBLAS("BLASgemv", BLASgemv_("T", &(_bbs), &(_bbs), &_mone, A, &(_bbs), w, &_ione, &_one, v…
166 #define PetscKernel_v_gets_v_plus_A_times_w(bs, v, A, w) \ argument
[all …]
/petsc/src/dm/impls/plex/tests/output/
H A Dex18_9_seq_hdf5_simple_int_err.out6 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 Dex7.c9 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/src/dm/dt/tests/output/
H A Dex7_1234.out7 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 Dsfwindow.c72 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/share/petsc/datafiles/meshes/
H A Dsurfacesphere_bin.msh7 …��������������������������+^Do�?۰/7���?��������������w��?�,u#��?�������� �����u#��?֗�w��?�…
8 ����07���?NDo�?�������� ���۰/7��ӿ+^Do�?�������� ����,u#�⿨��w��?��������
9 …w�鿅�u#��?�����������NDo07���?�����������+^Do�۰/7��ӿ��������������w���,u#������������…
/petsc/lib/petsc/bin/maint/generators/
H A Djava.py120 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]*)\)',…
148 f.write(w)
[all …]
/petsc/src/sys/classes/draw/impls/x/
H A Dxtext.c19 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()
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()
126 font->font_w = w; in PetscDrawXiMatchFontSize()
[all …]
H A Dxops.c10 #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 …]
/petsc/src/binding/petsc4py/test/
H A Dtest_nsp.py35 w = v.copy()
36 self.nullsp.remove(w)
37 return (v, w)
40 v, w = self._remove()
42 self.assertTrue(allclose(w.array, [0, 0, 9]))
43 del v, w
46 v, w = self._remove()
48 self.assertTrue(v.equal(w))
49 del v, w
56 v, w = self._remove()
[all …]
/petsc/src/dm/impls/da/
H A Ddasub.c119 PetscCall(PetscMalloc1(dd->w, &indices)); in DMDAGetRay()
120 indices[0] = dd->w * gp; in DMDAGetRay()
121 for (i = 1; i < dd->w; ++i) indices[i] = indices[i - 1] + 1; in DMDAGetRay()
122 PetscCall(AOApplicationToPetsc(ao, dd->w, indices)); in DMDAGetRay()
124 PetscCall(VecSetBlockSize(*newvec, dd->w)); in DMDAGetRay()
125 PetscCall(VecSetSizes(*newvec, dd->w, PETSC_DETERMINE)); in DMDAGetRay()
127 PetscCall(ISCreateGeneral(PETSC_COMM_SELF, dd->w, indices, PETSC_OWN_POINTER, &is)); in DMDAGetRay()
134 PetscCall(PetscMalloc1(dd->w * dd->M, &indices)); in DMDAGetRay()
135 indices[0] = gp * dd->M * dd->w; in DMDAGetRay()
136 for (i = 1; i < dd->M * dd->w; i++) indices[i] = indices[i - 1] + 1; in DMDAGetRay()
[all …]
/petsc/src/dm/impls/patch/tests/output/
H A Dex1_3.out3 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 …]
/petsc/src/dm/tests/output/
H A Dex34_1.out3 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/sys/classes/draw/impls/image/
H A Ddrawimage.c13 #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()
221 static PetscErrorCode PetscDrawStringGetSize_Image(PetscDraw draw, PetscReal *w, PetscReal *h) in PetscDrawStringGetSize_Image() argument
[all …]
/petsc/src/vec/vec/impls/mpi/kokkos/
H A Dmpikok.kokkos.cxx283 static PetscErrorCode VecDuplicateVecs_MPIKokkos_GEMV(Vec w, PetscInt m, Vec *V[]) in VecDuplicateVecs_MPIKokkos_GEMV() argument
288 Vec_MPI *wmpi = (Vec_MPI *)w->data; in VecDuplicateVecs_MPIKokkos_GEMV()
294 w->ops->duplicatevecs = VecDuplicateVecs_Default; in VecDuplicateVecs_MPIKokkos_GEMV()
295 PetscCall(VecDuplicateVecs(w, m, V)); in VecDuplicateVecs_MPIKokkos_GEMV()
298 PetscCall(VecGetLayout(w, &map)); in VecDuplicateVecs_MPIKokkos_GEMV()
299 VecGetLocalSizeAligned(w, 64, &lda); // get in lda the 64-bytes aligned local size in VecDuplicateVecs_MPIKokkos_GEMV()
311 PetscCall(PetscObjectListDuplicate(((PetscObject)w)->olist, &((PetscObject)v)->olist)); in VecDuplicateVecs_MPIKokkos_GEMV()
312 PetscCall(PetscFunctionListDuplicate(((PetscObject)w)->qlist, &((PetscObject)v)->qlist)); in VecDuplicateVecs_MPIKokkos_GEMV()
313 v->ops[0] = w->ops[0]; in VecDuplicateVecs_MPIKokkos_GEMV()
314 v->stash.donotstash = w->stash.donotstash; in VecDuplicateVecs_MPIKokkos_GEMV()
[all …]
/petsc/src/dm/impls/stag/tests/output/
H A Dex50_dim-3_coords-false.out3 Processor [0] M 4 N 3 P 3 m 1 n 1 p 1 w 3 s 1
117 Processor [0] M 4 N 3 P 3 m 1 n 1 p 1 w 1 s 1
159 Processor [0] M 4 N 3 P 3 m 1 n 1 p 1 w 1 s 1
201 Processor [0] M 4 N 3 P 2 m 1 n 1 p 1 w 3 s 1
279 Processor [0] M 4 N 3 P 2 m 1 n 1 p 1 w 1 s 1
309 Processor [0] M 4 N 3 P 2 m 1 n 1 p 1 w 1 s 1
339 Processor [0] M 4 N 3 P 2 m 1 n 1 p 1 w 1 s 1
369 Processor [0] M 4 N 2 P 3 m 1 n 1 p 1 w 3 s 1
447 Processor [0] M 4 N 2 P 3 m 1 n 1 p 1 w 1 s 1
477 Processor [0] M 4 N 2 P 3 m 1 n 1 p 1 w 1 s 1
[all …]
H A Dex50_dim-2_coords-true.out3 Processor [0] M 4 N 3 m 1 n 1 w 3 s 1
45 Processor [0] M 4 N 3 m 1 n 1 w 1 s 1
63 Processor [0] M 4 N 3 m 1 n 1 w 1 s 1
81 Processor [0] M 4 N 2 m 1 n 1 w 3 s 1
111 Processor [0] M 4 N 2 m 1 n 1 w 1 s 1
125 Processor [0] M 3 N 3 m 1 n 1 w 3 s 1
158 Processor [0] M 3 N 3 m 1 n 1 w 1 s 1
173 Processor [0] M 3 N 2 m 1 n 1 w 3 s 1
197 Processor [0] M 3 N 2 m 1 n 1 w 1 s 1
209 Processor [0] M 3 N 2 m 1 n 1 w 1 s 1
[all …]
H A Dex50_dim-2_coords-false.out3 Processor [0] M 4 N 3 m 1 n 1 w 3 s 1
45 Processor [0] M 4 N 3 m 1 n 1 w 1 s 1
63 Processor [0] M 4 N 3 m 1 n 1 w 1 s 1
81 Processor [0] M 4 N 2 m 1 n 1 w 3 s 1
111 Processor [0] M 4 N 2 m 1 n 1 w 1 s 1
125 Processor [0] M 3 N 3 m 1 n 1 w 3 s 1
158 Processor [0] M 3 N 3 m 1 n 1 w 1 s 1
173 Processor [0] M 3 N 2 m 1 n 1 w 3 s 1
197 Processor [0] M 3 N 2 m 1 n 1 w 1 s 1
209 Processor [0] M 3 N 2 m 1 n 1 w 1 s 1
[all …]
H A Dex50_dim-3_coords-true.out3 Processor [0] M 4 N 3 P 3 m 1 n 1 p 1 w 3 s 1
118 Processor [0] M 4 N 3 P 3 m 1 n 1 p 1 w 1 s 1
161 Processor [0] M 4 N 3 P 3 m 1 n 1 p 1 w 1 s 1
204 Processor [0] M 4 N 3 P 2 m 1 n 1 p 1 w 3 s 1
283 Processor [0] M 4 N 3 P 2 m 1 n 1 p 1 w 1 s 1
314 Processor [0] M 4 N 3 P 2 m 1 n 1 p 1 w 1 s 1
345 Processor [0] M 4 N 3 P 2 m 1 n 1 p 1 w 1 s 1
376 Processor [0] M 4 N 2 P 3 m 1 n 1 p 1 w 3 s 1
455 Processor [0] M 4 N 2 P 3 m 1 n 1 p 1 w 1 s 1
486 Processor [0] M 4 N 2 P 3 m 1 n 1 p 1 w 1 s 1
[all …]
/petsc/src/vec/vec/tutorials/
H A Dperformance.c10 Vec v, w, x; in main() local
25 PetscCall(VecDuplicate(v, &w)); in main()
27 PetscCall(VecSetRandom(w, rctx)); in main()
42 PetscCall(VecNorm(w, NORM_1, &norm1)); /* send w to GPU */ in main()
45 PetscCall(VecAXPY(w, 1.0, v)); in main()
48 PetscCall(VecDot(w, v, &val)); in main()
54 PetscCall(VecCopy(v, w)); in main()
59 PetscCall(VecDestroy(&w)); in main()
H A Dex1.c14 Vec x, y, w; /* vectors */ in main() local
39 PetscCall(VecDuplicate(x, &w)); in main()
82 PetscCall(VecCopy(x, w)); in main()
83 PetscCall(VecNorm(w, NORM_2, &norm)); in main()
110 PetscCall(VecWAXPY(w, two, x, y)); in main()
111 PetscCall(VecNorm(w, NORM_2, &norm)); in main()
116 PetscCall(VecPointwiseMult(w, y, x)); in main()
117 PetscCall(VecNorm(w, NORM_2, &norm)); in main()
122 PetscCall(VecPointwiseDivide(w, x, y)); in main()
123 PetscCall(VecNorm(w, NORM_2, &norm)); in main()
[all …]
/petsc/src/sys/classes/draw/utils/
H A Dimage.c11 …mageSavePPM(const char filename[], unsigned char palette[][3], unsigned int w, unsigned int h, con… in PetscDrawImageSavePPM() argument
24 int k, p, n = (int)(w * h); in PetscDrawImageSavePPM()
26 PetscCall(PetscMalloc1(3 * w * h, &rgb)); in PetscDrawImageSavePPM()
38 PetscCall(PetscSNPrintf(header, sizeof(header), "P6\n%d %d\n255\n%c", (int)w, (int)h, '\0')); in PetscDrawImageSavePPM()
42 PetscCall(PetscBinaryWrite(fd, rgb, 3 * w * h, PETSC_CHAR)); in PetscDrawImageSavePPM()
48 …ageSave_PPM(const char filename[], unsigned char palette[][3], unsigned int w, unsigned int h, con… in PetscDrawImageSave_PPM() argument
50 return PetscDrawImageSavePPM(filename, palette, w, h, pixels); in PetscDrawImageSave_PPM()
66 …mageSavePNG(const char filename[], unsigned char palette[][3], unsigned int w, unsigned int h, con… in PetscDrawImageSavePNG() argument
71 unsigned int row, stride = palette ? w : 3 * w; in PetscDrawImageSavePNG()
96 …png_set_IHDR(png_ptr, info_ptr, w, h, /*depth*/ 8, palette ? PNG_COLOR_TYPE_PALETTE : PNG_COLOR_TY… in PetscDrawImageSavePNG()
[all …]
/petsc/src/dm/tutorials/output/
H A Dex19_1.out4 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/ts/tutorials/
H A Dex43.c17 PetscReal w = PetscSqrtReal(1 - xi * xi) * omega; in Exact() local
18 PetscReal C1 = (v0 + a * u0) / w; in Exact()
20 u = PetscExpReal(-a * t) * (C1 * PetscSinReal(w * t) + C2 * PetscCosReal(w * t)); in Exact()
21 …* (C1 * PetscSinReal(w * t) + C2 * PetscCosReal(w * t)) + w * PetscExpReal(-a * t) * (C1 * PetscCo… in Exact()
23 PetscReal w = PetscSqrtReal(xi * xi - 1) * omega; in Exact() local
24 PetscReal C1 = (w * u0 + xi * u0 + v0) / (2 * w); in Exact()
25 PetscReal C2 = (w * u0 - xi * u0 - v0) / (2 * w); in Exact()
26 u = C1 * PetscExpReal((-xi + w) * t) + C2 * PetscExpReal((-xi - w) * t); in Exact()
27 …v = C1 * (-xi + w) * PetscExpReal((-xi + w) * t) + C2 * (-xi - w) * PetscExpReal((-xi -… in Exact()

12345678910>>...13