i am using the following code in my 3D FEM analysis
digits(16)
k_bc=vpa(k_bc);
U=k_bcf_bc.
residual=norm(k_bc*U-f_bc).
Now to decrease the residual i have taken the variable to symbolic domain.Now when i use vpa i loose my sparcity of stiffness matrix causing a slow calculation rate but give lower value of residual and when i use sparse matrix the residual value is high because of the large condition number of k_bc of the order 10e14(i think it is a result of lagrange shape function beig used in Galerkin’s FEM).how to use both sparcity and vpa simultaneously to have lowervalue to residual and still have higher speedi am using the following code in my 3D FEM analysis
digits(16)
k_bc=vpa(k_bc);
U=k_bcf_bc.
residual=norm(k_bc*U-f_bc).
Now to decrease the residual i have taken the variable to symbolic domain.Now when i use vpa i loose my sparcity of stiffness matrix causing a slow calculation rate but give lower value of residual and when i use sparse matrix the residual value is high because of the large condition number of k_bc of the order 10e14(i think it is a result of lagrange shape function beig used in Galerkin’s FEM).how to use both sparcity and vpa simultaneously to have lowervalue to residual and still have higher speed i am using the following code in my 3D FEM analysis
digits(16)
k_bc=vpa(k_bc);
U=k_bcf_bc.
residual=norm(k_bc*U-f_bc).
Now to decrease the residual i have taken the variable to symbolic domain.Now when i use vpa i loose my sparcity of stiffness matrix causing a slow calculation rate but give lower value of residual and when i use sparse matrix the residual value is high because of the large condition number of k_bc of the order 10e14(i think it is a result of lagrange shape function beig used in Galerkin’s FEM).how to use both sparcity and vpa simultaneously to have lowervalue to residual and still have higher speed vpa, sparse MATLAB Answers — New Questions

​