How allocate a matrix to grids?

9 views (last 30 days)
Hamideh Rouhani
Hamideh Rouhani on 19 Nov 2023
Commented: Matt J on 21 Nov 2023
I want to use staggered grid to write MATLAB code about viscoelastic fluid. You suppose I have a 10×14 node at x and y direction. I have initial value for conformation tensor as A=[1 , 0 ;0  1] I should allocate this value for all grids as matrix.  Then It should calculate eigen values of A matrix and also its normalized vectors. Then put eigen values in matrix as diagonal matrix as "Lamda". Then Z= Rt *Lambda*R Which RT is transpose of normalazed matrix of A. Now I do not know how should I define matrix (A) as the input , that the new one will replace with the previous in the loop. I defined as this: A (i,j)=eye(2) Or A11(i,j)=1, A12(i,j)0 , A21(i,j)=0 , A22(i,j)=1
But it gives error about dimension.
If I define A=eye(2), how MATLAB understand put this value in each node and in the next step replace it with new one which the value can be different from each grid to other one.

Sign in to answer this question.

Answers (1)

Matt J
Matt J on 19 Nov 2023
Edited: Matt J on 19 Nov 2023
Ran in:
You could have a 4 dimensional array
A=repmat(eye(2) ,1,1,10,14)
A =
A(:,:,1,1) = 1 0 0 1 A(:,:,2,1) = 1 0 0 1 A(:,:,3,1) = 1 0 0 1 A(:,:,4,1) = 1 0 0 1 A(:,:,5,1) = 1 0 0 1 A(:,:,6,1) = 1 0 0 1 A(:,:,7,1) = 1 0 0 1 A(:,:,8,1) = 1 0 0 1 A(:,:,9,1) = 1 0 0 1 A(:,:,10,1) = 1 0 0 1 A(:,:,1,2) = 1 0 0 1 A(:,:,2,2) = 1 0 0 1 A(:,:,3,2) = 1 0 0 1 A(:,:,4,2) = 1 0 0 1 A(:,:,5,2) = 1 0 0 1 A(:,:,6,2) = 1 0 0 1 A(:,:,7,2) = 1 0 0 1 A(:,:,8,2) = 1 0 0 1 A(:,:,9,2) = 1 0 0 1 A(:,:,10,2) = 1 0 0 1 A(:,:,1,3) = 1 0 0 1 A(:,:,2,3) = 1 0 0 1 A(:,:,3,3) = 1 0 0 1 A(:,:,4,3) = 1 0 0 1 A(:,:,5,3) = 1 0 0 1 A(:,:,6,3) = 1 0 0 1 A(:,:,7,3) = 1 0 0 1 A(:,:,8,3) = 1 0 0 1 A(:,:,9,3) = 1 0 0 1 A(:,:,10,3) = 1 0 0 1 A(:,:,1,4) = 1 0 0 1 A(:,:,2,4) = 1 0 0 1 A(:,:,3,4) = 1 0 0 1 A(:,:,4,4) = 1 0 0 1 A(:,:,5,4) = 1 0 0 1 A(:,:,6,4) = 1 0 0 1 A(:,:,7,4) = 1 0 0 1 A(:,:,8,4) = 1 0 0 1 A(:,:,9,4) = 1 0 0 1 A(:,:,10,4) = 1 0 0 1 A(:,:,1,5) = 1 0 0 1 A(:,:,2,5) = 1 0 0 1 A(:,:,3,5) = 1 0 0 1 A(:,:,4,5) = 1 0 0 1 A(:,:,5,5) = 1 0 0 1 A(:,:,6,5) = 1 0 0 1 A(:,:,7,5) = 1 0 0 1 A(:,:,8,5) = 1 0 0 1 A(:,:,9,5) = 1 0 0 1 A(:,:,10,5) = 1 0 0 1 A(:,:,1,6) = 1 0 0 1 A(:,:,2,6) = 1 0 0 1 A(:,:,3,6) = 1 0 0 1 A(:,:,4,6) = 1 0 0 1 A(:,:,5,6) = 1 0 0 1 A(:,:,6,6) = 1 0 0 1 A(:,:,7,6) = 1 0 0 1 A(:,:,8,6) = 1 0 0 1 A(:,:,9,6) = 1 0 0 1 A(:,:,10,6) = 1 0 0 1 A(:,:,1,7) = 1 0 0 1 A(:,:,2,7) = 1 0 0 1 A(:,:,3,7) = 1 0 0 1 A(:,:,4,7) = 1 0 0 1 A(:,:,5,7) = 1 0 0 1 A(:,:,6,7) = 1 0 0 1 A(:,:,7,7) = 1 0 0 1 A(:,:,8,7) = 1 0 0 1 A(:,:,9,7) = 1 0 0 1 A(:,:,10,7) = 1 0 0 1 A(:,:,1,8) = 1 0 0 1 A(:,:,2,8) = 1 0 0 1 A(:,:,3,8) = 1 0 0 1 A(:,:,4,8) = 1 0 0 1 A(:,:,5,8) = 1 0 0 1 A(:,:,6,8) = 1 0 0 1 A(:,:,7,8) = 1 0 0 1 A(:,:,8,8) = 1 0 0 1 A(:,:,9,8) = 1 0 0 1 A(:,:,10,8) = 1 0 0 1 A(:,:,1,9) = 1 0 0 1 A(:,:,2,9) = 1 0 0 1 A(:,:,3,9) = 1 0 0 1 A(:,:,4,9) = 1 0 0 1 A(:,:,5,9) = 1 0 0 1 A(:,:,6,9) = 1 0 0 1 A(:,:,7,9) = 1 0 0 1 A(:,:,8,9) = 1 0 0 1 A(:,:,9,9) = 1 0 0 1 A(:,:,10,9) = 1 0 0 1 A(:,:,1,10) = 1 0 0 1 A(:,:,2,10) = 1 0 0 1 A(:,:,3,10) = 1 0 0 1 A(:,:,4,10) = 1 0 0 1 A(:,:,5,10) = 1 0 0 1 A(:,:,6,10) = 1 0 0 1 A(:,:,7,10) = 1 0 0 1 A(:,:,8,10) = 1 0 0 1 A(:,:,9,10) = 1 0 0 1 A(:,:,10,10) = 1 0 0 1 A(:,:,1,11) = 1 0 0 1 A(:,:,2,11) = 1 0 0 1 A(:,:,3,11) = 1 0 0 1 A(:,:,4,11) = 1 0 0 1 A(:,:,5,11) = 1 0 0 1 A(:,:,6,11) = 1 0 0 1 A(:,:,7,11) = 1 0 0 1 A(:,:,8,11) = 1 0 0 1 A(:,:,9,11) = 1 0 0 1 A(:,:,10,11) = 1 0 0 1 A(:,:,1,12) = 1 0 0 1 A(:,:,2,12) = 1 0 0 1 A(:,:,3,12) = 1 0 0 1 A(:,:,4,12) = 1 0 0 1 A(:,:,5,12) = 1 0 0 1 A(:,:,6,12) = 1 0 0 1 A(:,:,7,12) = 1 0 0 1 A(:,:,8,12) = 1 0 0 1 A(:,:,9,12) = 1 0 0 1 A(:,:,10,12) = 1 0 0 1 A(:,:,1,13) = 1 0 0 1 A(:,:,2,13) = 1 0 0 1 A(:,:,3,13) = 1 0 0 1 A(:,:,4,13) = 1 0 0 1 A(:,:,5,13) = 1 0 0 1 A(:,:,6,13) = 1 0 0 1 A(:,:,7,13) = 1 0 0 1 A(:,:,8,13) = 1 0 0 1 A(:,:,9,13) = 1 0 0 1 A(:,:,10,13) = 1 0 0 1 A(:,:,1,14) = 1 0 0 1 A(:,:,2,14) = 1 0 0 1 A(:,:,3,14) = 1 0 0 1 A(:,:,4,14) = 1 0 0 1 A(:,:,5,14) = 1 0 0 1 A(:,:,6,14) = 1 0 0 1 A(:,:,7,14) = 1 0 0 1 A(:,:,8,14) = 1 0 0 1 A(:,:,9,14) = 1 0 0 1 A(:,:,10,14) = 1 0 0 1
  4 Comments
Show 2 older comments
Matt J
Matt J on 21 Nov 2023
You are quite welcome, but please Accept-click the answer if it resolved your question.

Sign in to comment.

Categories

Find more on Programming in Help Center and File Exchange

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!