Solve symmetric positive definite system from CEMEF-UTN
[u0, u, P0] = SPDsolve(K0, K, P, cDOF, cDOFasc, nn)
Complete stiffness matrix of the structure [DOFn*nn x DOFn*nn]
Active stiffness matrix of the structure (reduced to active DOF) [DOFn*nn-cDOF x DOFn*nn-cDOF]
Nodal external loads vector (reduced to active DOF) [DOFn*nn-cDOF x 1]
Constrained degrees of freedom vector (descending order) [cDOF x 1]
Constrained degrees of freedom vector (ascending order) [cDOF x 1]
Number of nodes [1 x 1]
Complete nodal displacements vector [DOFn*nn x 1]
Active nodal displacements vector (reduced to active DOF) [DOFn*nn-cDOF x 1]
Complete nodal loads vector [DOFn*nn x 1]
Grupo de Simulacion y Mecanica Computacional,