> Very Many Thanks for your efforts on this Barry. The PLAPACK > website looks like it hasn't been updated since 2007. Maybe PLAPACK > is in need of some maintenance? You said "nonsquare", is plapack working for > you for square matrices ?
Yes, it works for square matrices. See ~petsc/src/mat/examples/tests/ex103.c and ex107.c Hong > > >>> >>> >>> >>> [0]PETSC ERROR: --------------------- Error Message >>> ---------------------------- >>> -------- >>> [0]PETSC ERROR: Error in external library! >>> [1]PETSC ERROR: [0]PETSC ERROR: --------------------- Error Message >>> ------------ >>> ------------------------ >>> [1]PETSC ERROR: Error in external library! >>> Due to aparent bugs in PLAPACK,this is not currently supported! >> ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ >> I could not get plapack to properly handle nonsquare matrices properly >> for these operations. >> I tried hard to debug but plapack is a mess of complexity, got very >> frustrated. If you take out the >> generation of this error message, so the code runs, then you can try to >> debug plapack >> (I am pretty sure the problem is in plapack, not in the PETSc interface). >> >> Barry >> >