Thank you for the prompt response. I checked and my usage of KSPSolve is correct.
Besides, the assertion line is: if (*mat->rmap->N* != *y->map->N*) SETERRQ2(PETSC_COMM_SELF,PETSC_ERR_ARG_SIZ,"Mat mat,Vec y: global dim %D %D",mat->rmap->N,y->map->N); On Wed, Nov 12, 2014 at 2:29 PM, Barry Smith <[email protected]> wrote: > > The calling sequence for KSPSolve is KSPSolve(KSP ksp,Vec b,Vec x), > perhaps you have the b and the x backwards in your call? > > Barry > > > On Nov 12, 2014, at 2:23 PM, Brian Yang <[email protected]> wrote: > > > > Hi, > > > > I am solving an over-determined linear system Ax=b and b=0. > > > > A size: 296856 x 14430 > > x (non-zero guess): 14430 > > b (zero): 296856 > > > > All of them are using seq mode to create, so one node will solve this. > However I got this error: > > > > [0]PETSC ERROR: Nonconforming object sizes! > > [0]PETSC ERROR: Mat mat,Vec y: global dim 296856 14430! > > > > [0]PETSC ERROR: > ------------------------------------------------------------------------ > > [0]PETSC ERROR: MatMult() line 2171 in src/mat/interface/matrix.c > > [0]PETSC ERROR: KSP_MatMult() line 204 in > src/ksp/ksp/impls/cg//datadb/raafat/src/petsc-3.4.3/include/petsc-private/kspimpl.h > > [0]PETSC ERROR: KSPSolve_CG() line 132 in src/ksp/ksp/impls/cg/cg.c > > [0]PETSC ERROR: KSPSolve() line 441 in src/ksp/ksp/interface/itfunc.c > > > > > > Any restrictions here or I misunderstood something? Thanks! > > > > -- > > Brian Yang > > U of Houston > > > > > > > > -- Brian Yang U of Houston
