Home
last modified time | relevance | path

Searched refs:PetscConj (Results 1 – 25 of 68) sorted by relevance

123

/petsc/src/vec/vec/impls/seq/
H A Ddvec2.c106 sum0 += x2 * PetscConj(yy0[2]); in VecMDot_Seq()
107 sum1 += x2 * PetscConj(yy1[2]); in VecMDot_Seq()
108 sum2 += x2 * PetscConj(yy2[2]); /* fall through */ in VecMDot_Seq()
111 sum0 += x1 * PetscConj(yy0[1]); in VecMDot_Seq()
112 sum1 += x1 * PetscConj(yy1[1]); in VecMDot_Seq()
113 sum2 += x1 * PetscConj(yy2[1]); /* fall through */ in VecMDot_Seq()
116 sum0 += x0 * PetscConj(yy0[0]); in VecMDot_Seq()
117 sum1 += x0 * PetscConj(yy1[0]); in VecMDot_Seq()
118 sum2 += x0 * PetscConj(yy2[0]); /* fall through */ in VecMDot_Seq()
134 …sum0 += x0 * PetscConj(yy0[0]) + x1 * PetscConj(yy0[1]) + x2 * PetscConj(yy0[2]) + x3 * PetscConj(… in VecMDot_Seq()
[all …]
/petsc/src/ksp/ksp/impls/cg/pipecg2/
H A Dpipecg2.c23 sumwu += PW[j] * PetscConj(PU[j]); in VecMergedDot_Private()
24 sumru += PR[j] * PetscConj(PU[j]); in VecMergedDot_Private()
25 sumuu += PU[j] * PetscConj(PU[j]); in VecMergedDot_Private()
30 sumwu += PW[j] * PetscConj(PU[j]); in VecMergedDot_Private()
31 sumru += PR[j] * PetscConj(PU[j]); in VecMergedDot_Private()
32 sumuu += PR[j] * PetscConj(PR[j]); in VecMergedDot_Private()
37 sumwu += PW[j] * PetscConj(PU[j]); in VecMergedDot_Private()
38 sumru += PR[j] * PetscConj(PU[j]); in VecMergedDot_Private()
68 sumwm += PW[j] * PetscConj(PM[j]); in VecMergedDot2_Private()
69 sumnm += PN[j] * PetscConj(PM[j]); in VecMergedDot2_Private()
[all …]
/petsc/src/vec/vec/impls/seq/ftn-kernels/
H A Dfmdot.F9023 sum1 = sum1 + x(i)*PetscConj(y1(i))
24 sum2 = sum2 + x(i)*PetscConj(y2(i))
25 sum3 = sum3 + x(i)*PetscConj(y3(i))
26 sum4 = sum4 + x(i)*PetscConj(y4(i))
45 sum1 = sum1 + x(i)*PetscConj(y1(i))
46 sum2 = sum2 + x(i)*PetscConj(y2(i))
47 sum3 = sum3 + x(i)*PetscConj(y3(i))
65 sum1 = sum1 + x(i)*PetscConj(y1(i))
66 sum2 = sum2 + x(i)*PetscConj(y2(i))
83 sum1 = sum1 + x(i)*PetscConj(y1(i))
/petsc/src/mat/tests/
H A Dex264.c50 if (i < n - 1) PetscCall(MatSetValue(R, i, i + 1, PetscConj(b), INSERT_VALUES)); in main()
51 if (i < n - 2) PetscCall(MatSetValue(R, i, i + 2, PetscConj(a), INSERT_VALUES)); in main()
70 PetscCall(MatScale(block[2], PetscConj(b))); in main()
73 PetscCall(MatScale(block[3], PetscConj(d))); in main()
/petsc/src/vec/vec/impls/seq/kokkos/
H A Dveckok.kokkos.cxx319 for (size_type j = 0; j < value_count; ++j) sum[j] += PetscConj(yv[j](i)) * xval; in operator ()()
410 …lsum0 += xv(i) * PetscConj(y0(i)); lsum1 += xv(i) * PetscConj(y1(i)); lsum2 += xv(i) * PetscConj(y… in VecMultiDot_Verbose()
411 …lsum4 += xv(i) * PetscConj(y4(i)); lsum5 += xv(i) * PetscConj(y5(i)); lsum6 += xv(i) * PetscConj(y… in VecMultiDot_Verbose()
441 …lsum0 += xv(i) * PetscConj(y0(i)); lsum1 += xv(i) * PetscConj(y1(i)); lsum2 += xv(i) * PetscConj(y… in VecMultiDot_Verbose()
442 …lsum4 += xv(i) * PetscConj(y4(i)); lsum5 += xv(i) * PetscConj(y5(i)); lsum6 += xv(i) * PetscConj(y… in VecMultiDot_Verbose()
449 …lsum0 += xv(i) * PetscConj(y0(i)); lsum1 += xv(i) * PetscConj(y1(i)); lsum2 += xv(i) * PetscConj(y… in VecMultiDot_Verbose()
450 lsum4 += xv(i) * PetscConj(y4(i)); lsum5 += xv(i) * PetscConj(y5(i)); in VecMultiDot_Verbose()
457 …lsum0 += xv(i) * PetscConj(y0(i)); lsum1 += xv(i) * PetscConj(y1(i)); lsum2 += xv(i) * PetscConj(y… in VecMultiDot_Verbose()
458 lsum4 += xv(i) * PetscConj(y4(i)); in VecMultiDot_Verbose()
465 …lsum0 += xv(i) * PetscConj(y0(i)); lsum1 += xv(i) * PetscConj(y1(i)); lsum2 += xv(i) * PetscConj(y… in VecMultiDot_Verbose()
[all …]
/petsc/src/mat/impls/baij/seq/
H A Dbaij2.c2863 z[rval] += PetscConj(*v) * x1; in MatMultHermitianTransposeAdd_SeqBAIJ()
2879 z[rval++] += PetscConj(v[0]) * x1 + PetscConj(v[1]) * x2; in MatMultHermitianTransposeAdd_SeqBAIJ()
2880 z[rval++] += PetscConj(v[2]) * x1 + PetscConj(v[3]) * x2; in MatMultHermitianTransposeAdd_SeqBAIJ()
2897 z[rval++] += PetscConj(v[0]) * x1 + PetscConj(v[1]) * x2 + PetscConj(v[2]) * x3; in MatMultHermitianTransposeAdd_SeqBAIJ()
2898 z[rval++] += PetscConj(v[3]) * x1 + PetscConj(v[4]) * x2 + PetscConj(v[5]) * x3; in MatMultHermitianTransposeAdd_SeqBAIJ()
2899 z[rval++] += PetscConj(v[6]) * x1 + PetscConj(v[7]) * x2 + PetscConj(v[8]) * x3; in MatMultHermitianTransposeAdd_SeqBAIJ()
2917 …z[rval++] += PetscConj(v[0]) * x1 + PetscConj(v[1]) * x2 + PetscConj(v[2]) * x3 + PetscConj(v[3]) … in MatMultHermitianTransposeAdd_SeqBAIJ()
2918 …z[rval++] += PetscConj(v[4]) * x1 + PetscConj(v[5]) * x2 + PetscConj(v[6]) * x3 + PetscConj(v[7]) … in MatMultHermitianTransposeAdd_SeqBAIJ()
2919 …z[rval++] += PetscConj(v[8]) * x1 + PetscConj(v[9]) * x2 + PetscConj(v[10]) * x3 + PetscConj(v[11]… in MatMultHermitianTransposeAdd_SeqBAIJ()
2920 …z[rval++] += PetscConj(v[12]) * x1 + PetscConj(v[13]) * x2 + PetscConj(v[14]) * x3 + PetscConj(v[1… in MatMultHermitianTransposeAdd_SeqBAIJ()
[all …]
/petsc/src/ksp/ksp/impls/gmres/fgmres/
H A Dfgmres.c321 *hh = PetscConj(*cc) * tt + *ss * *(hh + 1); in KSPFGMRESUpdateHessenberg()
340 tt = PetscSqrtScalar(PetscConj(*hh) * *hh + PetscConj(*(hh + 1)) * *(hh + 1)); in KSPFGMRESUpdateHessenberg()
351 *RS(it) = PetscConj(*cc) * *RS(it); in KSPFGMRESUpdateHessenberg()
352 *hh = PetscConj(*cc) * *hh + *ss * *(hh + 1); in KSPFGMRESUpdateHessenberg()
/petsc/src/ksp/ksp/impls/bicg/
H A Dbicg.c73 b = PetscConj(b); in KSPSolve_BiCG()
85 ma = PetscConj(ma); in KSPSolve_BiCG()
/petsc/src/ksp/ksp/impls/gmres/lgmres/
H A Dlgmres.c456 *hh = PetscConj(*cc) * tt + *ss * *(hh + 1); in KSPLGMRESUpdateHessenberg()
475 tt = PetscSqrtScalar(PetscConj(*hh) * *hh + PetscConj(*(hh + 1)) * *(hh + 1)); in KSPLGMRESUpdateHessenberg()
485 *GRS(it) = PetscConj(*cc) * *GRS(it); in KSPLGMRESUpdateHessenberg()
486 *hh = PetscConj(*cc) * *hh + *ss * *(hh + 1); in KSPLGMRESUpdateHessenberg()
/petsc/src/mat/impls/cdiagonal/
H A Dcdiagonal.c71 PetscCall(VecAXPBY(v3, PetscConj(ctx->diag), 1.0, v1)); in MatMultHermitianTransposeAdd_ConstantDiagonal()
73 PetscCall(VecAXPBYPCZ(v3, PetscConj(ctx->diag), 1.0, 0.0, v1, v2)); in MatMultHermitianTransposeAdd_ConstantDiagonal()
160 PetscCall(VecAXPBY(y, PetscConj(ctx->diag), 0.0, x)); in MatMultHermitianTranspose_ConstantDiagonal()
205 ctx->diag = PetscConj(ctx->diag); in MatConjugate_ConstantDiagonal()
/petsc/include/petsc/finclude/
H A Dpetscsysbase.h106 #define PetscConj(a) conjg(a) macro
117 #define PetscConj(a) a
/petsc/src/ksp/ksp/tests/
H A Dex87.c69 if (i < n - 1) PetscCall(MatSetValue(R, i, i + 1, PetscConj(b), INSERT_VALUES)); in main()
70 if (i < n - 2) PetscCall(MatSetValue(R, i, i + 2, PetscConj(a), INSERT_VALUES)); in main()
H A Dex49.c46 PetscCall(MatSetValue(A, J, Ii, PetscConj(v), INSERT_VALUES)); in main()
/petsc/include/petsc/private/
H A Ddmpleximpl.h499 y[0 * ldx] = PetscConj(A[0]) * z[0] + PetscConj(A[3]) * z[1] + PetscConj(A[6]) * z[2]; in DMPlex_MultTranspose3D_Internal()
500 y[1 * ldx] = PetscConj(A[1]) * z[0] + PetscConj(A[4]) * z[1] + PetscConj(A[7]) * z[2]; in DMPlex_MultTranspose3D_Internal()
501 y[2 * ldx] = PetscConj(A[2]) * z[0] + PetscConj(A[5]) * z[1] + PetscConj(A[8]) * z[2]; in DMPlex_MultTranspose3D_Internal()
794 for (d = 0; d < dim; ++d) sum += PetscRealPart(PetscConj(x[d] - y[d]) * (x[d] - y[d])); in DMPlex_DistD_Internal()
/petsc/src/ksp/pc/impls/bjacobi/bjkokkos/
H A Dbjkokkos.kokkos.cxx164 …e(team, Nblk), [=](const int idx, PetscScalar &lsum) { lsum += R[idx] * PetscConj(R[idx]); }, dpi); in BJSolve_TFQMR()
192 …nge(team, Nblk), [=](const int idx, PetscScalar &dot) { dot += R[idx] * PetscConj(RP[idx]); }, rho… in BJSolve_TFQMR()
206 …ange(team, Nblk), [=](const int idx, PetscScalar &dot) { dot += V[idx] * PetscConj(RP[idx]); }, s); in BJSolve_TFQMR()
227 …e(team, Nblk), [=](const int idx, PetscScalar &lsum) { lsum += R[idx] * PetscConj(R[idx]); }, dpi); in BJSolve_TFQMR()
284 …nge(team, Nblk), [=](const int idx, PetscScalar &dot) { dot += R[idx] * PetscConj(RP[idx]); }, rho… in BJSolve_TFQMR()
396 …(team, Nblk), [=](const int idx, PetscScalar &lsum) { lsum += Rr[idx] * PetscConj(Rr[idx]); }, dpi… in BJSolve_BICG()
416 …ge(team, Nblk), [=](const int idx, PetscScalar &dot) { dot += Zr[idx] * PetscConj(Rl[idx]); }, bet… in BJSolve_BICG()
436 t2 = PetscConj(t1); in BJSolve_BICG()
448 …(team, Nblk), [=](const int idx, PetscScalar &lsum) { lsum += Zr[idx] * PetscConj(Pl[idx]); }, dpi… in BJSolve_BICG()
457 t2 = PetscConj(t1); in BJSolve_BICG()
[all …]
/petsc/src/ksp/ksp/impls/gmres/
H A Dborthog2.c86 for (j = 0; j <= it; j++) hnrm += PetscRealPart(lhh[j] * PetscConj(lhh[j])); in KSPGMRESClassicalGramSchmidtOrthogonalization()
H A Dgmres.c362 *hh = PetscConj(*cc) * tt + *ss * *(hh + 1); in KSPGMRESUpdateHessenberg()
374 tt = PetscSqrtScalar(PetscConj(*hh) * *hh + PetscConj(*(hh + 1)) * *(hh + 1)); in KSPGMRESUpdateHessenberg()
383 *GRS(it) = PetscConj(*cc) * *GRS(it); in KSPGMRESUpdateHessenberg()
384 *hh = PetscConj(*cc) * *hh + *ss * *(hh + 1); in KSPGMRESUpdateHessenberg()
/petsc/src/ksp/ksp/impls/gmres/pgmres/
H A Dpgmres.c265 hh[j] = PetscConj(cc[j]) * hhj + ss[j] * hh[j + 1]; in KSPPGMRESUpdateHessenberg()
291 hh[it] = PetscConj(cc[it]) * hh[it] + ss[it] * hh[it + 1]; in KSPPGMRESUpdateHessenberg()
293 rs[it] = PetscConj(cc[it]) * rs[it]; in KSPPGMRESUpdateHessenberg()
/petsc/src/mat/impls/aij/seq/cholmod/
H A Daijcholmod.c26 if (values) ca[k] = PetscConj(aa[j]); in MatWrapCholmod_seqaij()
/petsc/src/ksp/ksp/impls/gmres/pipefgmres/
H A Dpipefgmres.c407 hh[j] = PetscConj(cc[j]) * hhj + ss[j] * hh[j + 1]; in KSPPIPEFGMRESUpdateHessenberg()
433 hh[it] = PetscConj(cc[it]) * hh[it] + ss[it] * hh[it + 1]; in KSPPIPEFGMRESUpdateHessenberg()
435 rs[it] = PetscConj(cc[it]) * rs[it]; in KSPPIPEFGMRESUpdateHessenberg()
/petsc/src/dm/dt/tests/
H A Dex5.c100 err += PetscRealPart(PetscConj(diff) * diff); in main()
/petsc/src/snes/impls/vi/
H A Dvi.c118 …]) - zerotolerance) || PetscRealPart(f[i]) >= 0.0)) rnorm += PetscRealPart(PetscConj(f[i]) * f[i]); in SNESMonitorVI()
339 …) - zerotolerance) || PetscRealPart(f[i]) >= 0.0)) *fnorm += PetscRealPart(PetscConj(f[i]) * f[i]); in SNESVIComputeInactiveSetFnorm()
384 …tscRealPart(xu[i]) - zerotolerance) || PetscRealPart(f[i]) >= 0.0)) *fty += f[i] * PetscConj(y[i]); in SNESVIComputeInactiveSetFtY()
/petsc/src/sys/tests/
H A Dex25.c80 b = PetscConj(a); in main()
/petsc/src/ksp/ksp/impls/gmres/dgmres/
H A Ddgmres.c379 *hh = PetscConj(*cc) * tt + *ss * *(hh + 1); in KSPDGMRESUpdateHessenberg()
391 tt = PetscSqrtScalar(PetscConj(*hh) * *hh + PetscConj(*(hh + 1)) * *(hh + 1)); in KSPDGMRESUpdateHessenberg()
399 *GRS(it) = PetscConj(*cc) * *GRS(it); in KSPDGMRESUpdateHessenberg()
400 *hh = PetscConj(*cc) * *hh + *ss * *(hh + 1); in KSPDGMRESUpdateHessenberg()
/petsc/src/ksp/ksp/tutorials/
H A Dex13.c108 …for (i = 0; i < N; i++) enorm += PetscRealPart(PetscConj(solution[i] - userx[i]) * (solution[i] - … in main()

123