Home
last modified time | relevance | path

Searched refs:agmres (Results 1 – 4 of 4) sorted by relevance

/petsc/src/ksp/ksp/impls/gmres/agmres/
H A Dagmres.c27 KSP_AGMRES *agmres = (KSP_AGMRES *)ksp->data; in KSPSetUp_AGMRES() local
28 PetscInt neig = agmres->neig; in KSPSetUp_AGMRES()
29 const PetscInt max_k = agmres->max_k; in KSPSetUp_AGMRES()
36 agmres->q_preallocate = PETSC_TRUE; /* No allocation on the fly */ in KSPSetUp_AGMRES()
39 agmres->max_k = N; /* Set the augmented size to be allocated in KSPSetup_GMRES */ in KSPSetUp_AGMRES()
41 agmres->max_k = max_k; in KSPSetUp_AGMRES()
45 …PetscCall(PetscCalloc4(max_k, &agmres->Rshift, max_k, &agmres->Ishift, hes, &agmres->Rloc, (N + 1)… in KSPSetUp_AGMRES()
46 PetscCall(PetscMalloc3(N + 1, &agmres->tau, lwork, &agmres->work, N + 1, &agmres->nrs)); in KSPSetUp_AGMRES()
47 …PetscCall(PetscCalloc4(N + 1, &agmres->Scale, N + 1, &agmres->sgn, N + 1, &agmres->tloc, N + 1, &a… in KSPSetUp_AGMRES()
50 PetscCall(VecGetLocalSize(agmres->vecs[0], &nloc)); in KSPSetUp_AGMRES()
[all …]
H A Dagmresimpl.h42 #define VEC_TMP agmres->vecs[0]
43 #define VEC_TMP_MATOP agmres->vecs[1]
44 #define VEC_V(i) agmres->vecs[VEC_OFFSET + i]
46 #define MAXKSPSIZE (agmres->DeflPrecond ? agmres->max_k : (agmres->max_k + agmres->max_neig))
47 #define KSPSIZE (agmres->DeflPrecond ? agmres->max_k : (agmres->max_k + agmres->r))
48 #define H(a, b) (agmres->hh_origin + (b) * (MAXKSPSIZE + 2) + (a))
49 #define HS(a, b) (agmres->hes_origin + (b) * (MAXKSPSIZE + 1) + (a))
50 #define RLOC(a, b) (agmres->Rloc + (b) * (MAXKSPSIZE + 1) + (a))
H A Dagmresdeflation.c94 KSP_AGMRES *agmres = (KSP_AGMRES *)ksp->data; in KSPAGMRESSchurForm() local
95 PetscInt max_k = agmres->max_k; in KSPAGMRESSchurForm()
97 PetscInt neig = agmres->neig; in KSPAGMRESSchurForm()
98 PetscScalar *wr = agmres->wr; in KSPAGMRESSchurForm()
99 PetscScalar *wi = agmres->wi; in KSPAGMRESSchurForm()
100 PetscScalar *beta = agmres->beta; in KSPAGMRESSchurForm()
101 PetscScalar *Q = agmres->Q; in KSPAGMRESSchurForm()
102 PetscScalar *Z = agmres->Z; in KSPAGMRESSchurForm()
103 PetscScalar *work = agmres->work; in KSPAGMRESSchurForm()
104 PetscBLASInt *select = agmres->select; in KSPAGMRESSchurForm()
[all …]
H A Dagmresorthog.c22 KSP_AGMRES *agmres = (KSP_AGMRES *)ksp->data; in KSPAGMRESRoddecInitNeighboor() local
26 PetscCall(PetscObjectGetComm((PetscObject)agmres->vecs[0], &comm)); in KSPAGMRESRoddecInitNeighboor()
33 agmres->Ileft = rank - 1; in KSPAGMRESRoddecInitNeighboor()
34 agmres->Iright = rank + 1; in KSPAGMRESRoddecInitNeighboor()
37 agmres->Ileft = rank - 1; in KSPAGMRESRoddecInitNeighboor()
38 agmres->Iright = First; in KSPAGMRESRoddecInitNeighboor()
41 agmres->Ileft = Last; in KSPAGMRESRoddecInitNeighboor()
42 agmres->Iright = rank + 1; in KSPAGMRESRoddecInitNeighboor()
46 agmres->rank = rank; in KSPAGMRESRoddecInitNeighboor()
47 agmres->size = size; in KSPAGMRESRoddecInitNeighboor()
[all …]