Note: This discussion is about an older version of the COMSOL Multiphysics® software. The information provided may be out of date.

Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.

Question regarding assembled matrices

Please login with a confirmed email address before reporting spam

Hi everyone

I have question about the assemble command
[K,L,M,N] = assemble(fem);
which actually return the matrices of the linear system being solved.
If I have got it right the M corresponds to the values of dirichlet BC.
What I want to do is to change iteratively the values of M i.e. change my constant BC and resolve using
fem.sol = femstatic('In',{'K' K 'L' L 'M' M 'N' N});
without assembling the matrices again as the problem is rather large 1e6 dof.

My issue now is that the format of M is nx1 where n are the number of nodes with the imposed BC.

After I solve the system fem.sol.u has of course the size of ndof x 1.
How can I identify the n elements of fem.sol.u which are actually constant?

Thank you in advance

Giorgos


0 Replies Last Post 11 ott 2011, 20:05 GMT-4
COMSOL Moderator

Hello Giorgos Kourakos

Your Discussion has gone 30 days without a reply. If you still need help with COMSOL and have an on-subscription license, please visit our Support Center for help.

If you do not hold an on-subscription license, you may find an answer in another Discussion or in the Knowledge Base.

Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.