fixed a bug in amg coarsening

jiyess · Aug 30, 2021 · 912c396 · 912c396
1 parent 11f05d1
commit 912c396
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/solver/coarsenAMGc.m b/solver/coarsenAMGc.m
@@ -27,7 +27,7 @@
 if ~exist('theta','var'), theta = 0.025; end
 N = size(A,1);
 isC = false(N,1);       % C: coarse node
-N0 = min(sqrt(N),25);   % number of the coarest nodes
+N0 = min(floor(sqrt(N)),25);   % number of the coarest nodes
 
 %% Generate strong connectness matrix
 Dinv = spdiags(1./sqrt(diag(A)),0,N,N);