CLMatCont5p4 What is happening in this part of init_EP_EP code?

[Anonymous]

///////////////
jac = cjac(eds.func,eds.Jacobian,x,num2cell(p),eds.ActiveParams);
jac = jac(:,1:cds.ndim-1);
jac(:,end+1) = 0;
A1 = sparse(eds.BiAlt_M1_I,eds.BiAlt_M1_J,jac(eds.BiAlt_M1_V));
A2 = sparse(eds.BiAlt_M2_I,eds.BiAlt_M2_J,jac(eds.BiAlt_M2_V));
A3 = sparse(eds.BiAlt_M3_I,eds.BiAlt_M3_J,jac(eds.BiAlt_M3_V));
A = A1-A2+A3;
[Q,R,E] = qr(full(A));
eds.bigW = Q(:,end);
eds.bigV = E(:,end);
eds.bigD = 0;
////////////////////////

What does this do? I had to define a large system and now I want to perform continuation on it. I get an error that says---

////////
Index in position 2 exceeds array bounds
(must not exceed 3).

Error in init_EP_EP (line 79)
jac = jac(:,1:cds.ndim-1);

Error in continue_my_func (line 6)
[x0,v0]=init_EP_EP(@my_func,X0,p,ap);
/////////

I think I might be able to solve it if I understood what that part of the code does.
Also, I have defined the jacobian for my system, what does cjac do in this case?

[hilmeijer]

cjac will compute the first order partial derivatives with respect to the state variables as well as to the Active Parameters. So if the input is ok, then this will result in a n(n+1) (rowscolumns) array.

---

From: Rohini Vaideswaran rohini-v@users.sourceforge.net
Sent: Sunday, July 7, 2019 4:29 PM
To: [matcont:discussion]
Subject: [matcont:discussion] CLMatCont5p4 What is happening in this part of init_EP_EP code?

///////////////
jac = cjac(eds.func,eds.Jacobian,x,num2cell(p),eds.ActiveParams);
jac = jac(:,1:cds.ndim-1);
jac(:,end+1) = 0;
A1 = sparse(eds.BiAlt_M1_I,eds.BiAlt_M1_J,jac(eds.BiAlt_M1_V));
A2 = sparse(eds.BiAlt_M2_I,eds.BiAlt_M2_J,jac(eds.BiAlt_M2_V));
A3 = sparse(eds.BiAlt_M3_I,eds.BiAlt_M3_J,jac(eds.BiAlt_M3_V));
A = A1-A2+A3;
[Q,R,E] = qr(full(A));
eds.bigW = Q(:,end);
eds.bigV = E(:,end);
eds.bigD = 0;
////////////////////////

What does this do? I had to define a large system and now I want to perform continuation on it. I get an error that says---

////////
Index in position 2 exceeds array bounds
(must not exceed 3).

Error in init_EP_EP (line 79)
jac = jac(:,1:cds.ndim-1);

Error in continue_boom_ns (line 6)
[x0,v0]=init_EP_EP(@boom_ns,X0,p,ap);
/////////

I think I might be able to solve it if I understood what that part of the code does.

---

CLMatCont5p4 What is happening in this part of init_EP_EP code?

---

Sent from sourceforge.net because you indicated interest in https://sourceforge.net/p/matcont/discussion/762214/

To unsubscribe from further messages, please visit https://sourceforge.net/auth/subscriptions/