Mark When I run with exactly the options in this file on either 1 or 8 process and MatView() the matrix you pass to KSP the matrix has ALL zeros. So of course the GAMG solver is not happy. Is the matrix suppose to have all zeros? When I run the command line options in the make file for runex56 the matrix has nonzero values in it and GAMG runs fine. Barry
On Jun 27, 2016, at 6:56 AM, Mark Adams <[email protected]> wrote:
Barry,
I am having problems with your branch on ksp/ex56. First I was getting an error with the increase of the overlap that you added (I don't want that to be the default). So I set this back to 0 in gamg.c.
Next, I get the same error with the vec out of range because ex56 has BCs in the matrix. I created a branch mark/fix-gamg-asm-aggs from your branch and added the dummy ASM block. But I am getting all kinds of errors. I get errors in get MatSubMatrix and a funny error when it tries to free a prefix string.
I've attached the valgrind output of one process, of 8. There is a lot of noise with string operations but it looks like the first non-string thing is in MatGetSubMatrices. ex56 is elasticity. Maybe there is a block size problem?
<out_val_1>