zero.b

来自「C 开发的有限元软件」· B 代码 · 共 16 行

function ZeroRowCol ( )
    global K, nodes, dof_numbers, dof_positions

    for n in nodes do
	b = nd * (n.number - 1) + 1
	for d in dof_numbers do
	    if n.constraint.dx (d) then
		a = b + dof_positions (d) - 1
		K (:, a) = zeros (rols (K))
		K (a, :) = zeros (cols (K))
		K (a, a) = 1
	    end
	end
    end
end