Thank you both for your response. From: Barry Smith <bsm...@petsc.dev> Date: Thursday, 5 September 2024 at 23:19 To: Chockalingam, Karthikeyan (STFC,DL,HC) <karthikeyan.chockalin...@stfc.ac.uk> Cc: petsc-users@mcs.anl.gov <petsc-users@mcs.anl.gov> Subject: Re: [petsc-users] MPIAIJ MatMult and non conforming object sizes
You can use MatCreateVecs() to create vectors sized properly for both x and y in y = Ax when A is not square. On Sep 5, 2024, at 5:20 PM, Karthikeyan Chockalingam - STFC UKRI via petsc-users <petsc-users@mcs.anl.gov> wrote: Hello, I am unsure why the program crashes even while running the code serially petscErr = MatCreate(mesh.comm().get(), &par_G); petscErr = MatSetType(par_G, MATMPIAIJ); petscErr = MatSetSizes(par_G, local_num_rows, local_num_cols, total_num_rows, total_num_cols); PetscInt d_nz = 2; PetscInt o_nz = 2; petscErr = MatMPIAIJSetPreallocation(par_G, d_nz, NULL, o_nz, NULL); petscErr = VecCreateMPI(mesh.comm().get(),local_num_cols,total_num_cols,&par_xcoord); petscErr = VecDuplicate(par_xcoord,&par_xvec); petscErr = MatMult(par_G, par_xcoord, par_xvec); [0]PETSC ERROR: --------------------- Error Message -------------------------------------------------------------- [0]PETSC ERROR: Nonconforming object sizes [0]PETSC ERROR: Mat mat,Vec y: global dim 40 25 [0]PETSC ERROR: WARNING! There are unused option(s) set! Could be the program crashed before usage or a spelling mistake, etc! [0]PETSC ERROR: Option left: name:--keep-cout (no value) source: command line [0]PETSC ERROR: See https://urldefense.us/v3/__https://petsc.org/release/faq/__;!!G_uCfscf7eWS!dkFKYMJODzVEEvczpl2dC8oOLyzskTQ_Q3FJm8ewb0TiuxAwrMXwH1eMCstakS6pUT5C04m47ZecPpQFRWqJ1jRXwEU1LXKGK3E8$ <https://urldefense.us/v3/__https:/petsc.org/release/faq/__;!!G_uCfscf7eWS!eNf_b1At1lAvdYeAlplIXxQN9WwVYKe1gkVUZRvNkdF0V6JCbVOG_mYTuCQrp2rs8FMod_waFjX0LLNuAvMcHmdf7Hrzm4s3JaQT$> for trouble shooting. [0]PETSC ERROR: Petsc Release Version 3.20.3, unknown [0]PETSC ERROR: ./example-dbg on a arch-moose named HC20210312 by karthikeyan.chockalingam Thu Sep 5 22:16:54 2024 [0]PETSC ERROR: Configure options --with-64-bit-indices --with-cxx-dialect=C++17 --with-debugging=no --with-fortran-bindings=0 --with-mpi=1 --with-openmp=1 --with-shared-libraries=1 --with-sowing=0 --download-fblaslapack=1 --download-hypre=1 --download-metis=1 --download-mumps=1 --download-ptscotch=1 --download-parmetis=1 --download-scalapack=1 --download-slepc=1 --download-strumpack=1 --download-superlu_dist=1 --download-bison=1 --download-hdf5=1 --with-hdf5-fortran-bindings=0 --download-hdf5-configure-arguments="--with-zlib" --with-make-np=4 [0]PETSC ERROR: #1 MatMult() at /Users/karthikeyan.chockalingam/moose/petsc/src/mat/interface/matrix.c:2587 Vec Object: 1 MPI process Many thanks for the help! Karthik. -- Karthik Chockalingam, Ph.D. Senior Research Software Engineer High Performance Systems Engineering Group Hartree Centre | Science and Technology Facilities Council karthikeyan.chockalin...@stfc.ac.uk<mailto:karthikeyan.chockalin...@stfc.ac.uk> <image001.png>
