| /petsc/src/mat/impls/baij/seq/ |
| H A D | baij.c | 223 PetscScalar s[7], xw[7] = {0}; /* avoid some compilers thinking xw is uninitialized */ in MatSOR_SeqBAIJ() local
265 xw[0] = x[vi[j]]; in MatSOR_SeqBAIJ()
266 PetscKernel_v_gets_v_minus_A_times_w_1(s, (v + j), xw); in MatSOR_SeqBAIJ()
269 PetscKernel_v_gets_A_times_w_1(xw, idiag, s); in MatSOR_SeqBAIJ()
270 x[i2] = xw[0]; in MatSOR_SeqBAIJ()
290 xw[0] = x[idx]; in MatSOR_SeqBAIJ()
291 xw[1] = x[1 + idx]; in MatSOR_SeqBAIJ()
292 PetscKernel_v_gets_v_minus_A_times_w_2(s, (v + it), xw); in MatSOR_SeqBAIJ()
296 PetscKernel_v_gets_A_times_w_2(xw, idiag, s); in MatSOR_SeqBAIJ()
297 x[i2] = xw[0]; in MatSOR_SeqBAIJ()
[all …]
|
| /petsc/src/dm/dt/interface/ |
| H A D | dt.c | 1915 PetscReal *x, *w, *xw, *ww; in PetscDTGaussTensorQuadrature() local
1942 PetscCall(PetscMalloc2(npoints, &xw, npoints, &ww)); in PetscDTGaussTensorQuadrature()
1943 PetscCall(PetscDTGaussQuadrature(npoints, a, b, xw, ww)); in PetscDTGaussTensorQuadrature()
1946 x[(i * npoints + j) * dim + 0] = xw[i]; in PetscDTGaussTensorQuadrature()
1947 x[(i * npoints + j) * dim + 1] = xw[j]; in PetscDTGaussTensorQuadrature()
1951 PetscCall(PetscFree2(xw, ww)); in PetscDTGaussTensorQuadrature()
1955 PetscCall(PetscMalloc2(npoints, &xw, npoints, &ww)); in PetscDTGaussTensorQuadrature()
1956 PetscCall(PetscDTGaussQuadrature(npoints, a, b, xw, ww)); in PetscDTGaussTensorQuadrature()
1960 x[((i * npoints + j) * npoints + k) * dim + 0] = xw[i]; in PetscDTGaussTensorQuadrature()
1961 x[((i * npoints + j) * npoints + k) * dim + 1] = xw[j]; in PetscDTGaussTensorQuadrature()
[all …]
|
| /petsc/src/ts/tutorials/hamiltonian/ |
| H A D | ex3.c | 613 PetscReal xv0[3], xJ[9], xinvJ[9], xdetJ, xqr[3], xden, xw = 0.; in InitializeWeights() local
630 xw += xden * (xwq[q] * xdetJ); in InitializeWeights()
632 xwtot += xw; in InitializeWeights()
639 …ering[c], (double)PetscRealPart(xcoords[0]), (double)PetscRealPart(xcoords[0 + dim]), (double)xw)); in InitializeWeights()
669 weight[p] = totalWeight * vw * xw; in InitializeWeights()
671 …" weight exceeded 10: weight=%g, xw=%g, vw=%g, totalWeight=%g", p, weight[p], xw, vw, totalWeight); in InitializeWeights()
673 …scCall(PetscPrintf(comm, "particle %" PetscInt_FMT ": %g, vw: %g xw: %g\n", p, weight[p], vw, xw)); in InitializeWeights()
|
| H A D | ex2.c | 1310 PetscReal xv0[3], xJ[9], xinvJ[9], xdetJ, xqr[3], xden, xw = 0.; in InitializeWeights() local
1327 xw += xden * (xwq[q] * xdetJ); in InitializeWeights()
1329 xwtot += xw; in InitializeWeights()
1336 …ering[c], (double)PetscRealPart(xcoords[0]), (double)PetscRealPart(xcoords[0 + dim]), (double)xw)); in InitializeWeights()
1354 weight[p] = totalWeight * vw * xw; in InitializeWeights()
1356 …ERR_ARG_WRONG, "Particle %" PetscInt_FMT " weight exceeded 1: %g, %g, %g", p, xw, vw, totalWeight); in InitializeWeights()
1358 …scCall(PetscPrintf(comm, "particle %" PetscInt_FMT ": %g, vw: %g xw: %g\n", p, weight[p], vw, xw)); in InitializeWeights()
|
| H A D | ex4.c | 1924 PetscReal xv0[3], xJ[9], xinvJ[9], xdetJ, xqr[3], xden, xw = 0.; in InitializeWeights() local
1941 xw += xden * (xwq[q] * xdetJ); in InitializeWeights()
1943 xwtot += xw; in InitializeWeights()
1950 …ering[c], (double)PetscRealPart(xcoords[0]), (double)PetscRealPart(xcoords[0 + dim]), (double)xw)); in InitializeWeights()
1968 weight[p] = totalWeight * vw * xw; in InitializeWeights()
1970 …ERR_ARG_WRONG, "Particle %" PetscInt_FMT " weight exceeded 1: %g, %g, %g", p, xw, vw, totalWeight); in InitializeWeights()
1972 …scCall(PetscPrintf(comm, "particle %" PetscInt_FMT ": %g, vw: %g xw: %g\n", p, weight[p], vw, xw)); in InitializeWeights()
|