Симметричная матрица не симметрична

Недавно я написал код для метода конечных элементов. Поскольку мой алгоритм генерирует симметричную матрицу, после сбора данных для каждого элемента результирующая матрица должна быть симметричной.

Однако после того, как я запустил цикл над элементом, результирующая глобальная матрица не является симметричной. Базовая структура кода выглядит так:

A=zeros(dof,dof)

for (each element)
    loc_A = v'*(diagonal matrix)*v
         % (v is 1xN row vector)
         % loc_A is symmetric matrix

    A(K,K) = A(K,K)+loc_A
                 % (K is Nx1 column vector)
end

Поскольку я добавляю симметричные матрицы, результирующая матрица также должна быть симметричной.

Однако полученная матрица не является симметричной (я проверял с помощью issymmetric(A)). Я думаю, что это из-за ошибки округления. Если я добавлю A=(A+A')/2, результирующая матрица будет симметричной (конечно...). Однако я не хочу делать такие вещи. Есть ли какое-либо другое средство, которое делает A симметричным естественным образом без какой-либо постобработки?