amgcl
amgcl copied to clipboard
Not divisible block size with nullspace vectors
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