I am writing a FEM code where I have created the stiffness matrix using sparse command. The problem is implementing boundary condition, no matter what I do, this part is taking much more time than other parts of the code. Is there a faster way?

I have attached the code. The generation of stiffness matrix is taking about 2 seconds, but implementing the boundary condition takes more than 10 seconds. I want to know how to speed up this section.