amgcl icon indicating copy to clipboard operation
amgcl copied to clipboard

Not divisible block size with nullspace vectors

Open njsyw1997 opened this issue 1 year ago • 5 comments

Hi @ddemidov , I am using AMGCL to solve a 2D linear elasticity problem. It has a 2 $\times$ 2 block structure and works fine without nullspace vectors. But when I pass the coordinates to the solver, it returns what(): Matrix size should be divisible by block_size Do you have any idea of this problem? Here are the parameters and sample data.

./amgcl/build/examples/solver \   
-A ./stiffness.mtx \
-f ./rhs.mtx \
-C ./vecpoints.mtx \
-p precond.coarsening.type="smoothed_aggregation" \
-p precond.coarsening.aggr.eps_strong=0.00 \
-p precond.relax.type="ilu0" \
-p solver.maxiter=1000 \
-p solver.type="cg" \
-b 2

data.zip

njsyw1997 avatar Mar 15 '23 14:03 njsyw1997