Advanced Search

mehdibgh

240 Reputation

6 Badges

7 years, 274 days
Contact mehdibgh

MaplePrimes Activity


These are replies submitted by mehdibgh

@tomleslie As I said first, Originaly my...

mehdibgh 240
March 28 2017
0 0

@tomleslie As I said first, Originaly my own matrix have 1000*1000 dimensions, I could not unassign them one by one.

May matrix have the variables as below.

var1 := [seq(seq(U[i, j], j = 0 .. JJ), i = 0 .. II), seq(seq(V[i, j], j = 0 .. JJ), i = 0 .. II), seq(seq(W[i, j], j = 0 .. JJ), i = 0 .. II), seq(seq(Phi[i, j], j = 0 .. JJ), i = 0 .. II), seq(seq(Xi[i, j], j = 0 .. JJ), i = 0 .. II)]

How to unassign them?

is there another simple way?

@Carl Love Where is the problem, this co...

mehdibgh 240
March 28 2017
0 0

@Carl Love Where is the problem, this command dont work at for loop?


 

``

> 

restart

Error, invalid input: with expects its 1st argument, pname, to be of type {`module`, package}, but received shareman

  
> 

A := Matrix([[x, y], [z, t]])

A := Matrix(2, 2, {(1, 1) = x, (1, 2) = y, (2, 1) = z, (2, 2) = t})

 (1)
> 

B := Matrix([[1, 2*f], [3, 4*f]])

B := Matrix(2, 2, {(1, 1) = 1, (1, 2) = 2*f, (2, 1) = 3, (2, 2) = 4*f})

 (2)
> 

r := 1:
> 

for f to 11 do C := A; `~`[assign](`~`[`=`](C, B)); F[r] := x+y+z+t; r := r+1 end do

Error, invalid left hand side in assignment

  

``

 

 

 

``


 

Download assgnSoal.mw

@tomleslie What do you mean by breakpoin...

mehdibgh 240
March 24 2017
0 0

@tomleslie What do you mean by breakpoint? Could you please give an example?

@tomleslie What about pause then resume ...

mehdibgh 240
March 24 2017
0 1

@tomleslie What about pause then resume the running?

 

@Preben Alsholm Is not there a simple wa...

mehdibgh 240
March 24 2017
0 0

@Preben Alsholm Is not there a simple way???

for example a short command (e.g. >stop)?

 

@tomleslie My matrices have variable par...

mehdibgh 240
March 23 2017
0 0

@tomleslie My matrices have variable parameter (x) within itself, it seems this command could not export matrices with variable entries.

Any other way to copy and paste these matrices???

@AmusingYeti corrected, please check aga...

mehdibgh 240
March 23 2017
0 1

@AmusingYeti corrected, please check again

let me explain more. In my worksheet th...

mehdibgh 240
March 23 2017
0 0

let me explain more.

In my worksheet there are 147*147 matrix was computated during long process as an output. Now I want to use that structure in another worksheet. I couldnt transfer it to that new worksheet. It seemsthe copy-paste dont work here to transfer, what do you recommend???

 

@John May Thanks, Are there are ways to ...

mehdibgh 240
March 13 2017
0 0

@John May Thanks, Are there are ways to cope with this problem?

This is because my matrix have variable dimensions and number of parameters that have created within the loop.

@acer ،Thanks your Kind answer. As you ...

mehdibgh 240
March 10 2017
0 5

@acer ،Thanks your Kind answer.

As you see in my previous post, my matrix M is positive definite and symmetric so it seems your second recommended approach is to be implemented by Maple. But as I said the order of matrix M__11 is varying between 200-1000+.

I tested Matlab sqrtm(A) command for such matrices, it returned the result up to 1 min for 1000+ ones and 3 seconds for 200.

But unfortunately couldnt receive any output by Maple. Please help me what exactly have to do to get the results by Maple?

It will be great if you test your note on my askkk.mw  file posted before.

Thanks,

@acer let me explain my problem moreAs y...

mehdibgh 240
March 10 2017
0 0

@acer let me explain my problem more

As you see in order to do that transmission I have to compute the square root of matrix K and M. As I said in my previous post the order of matrices are 200-1000, I couldnt compute the  square root of of them by Maple untill now. So please help if is it possible. Otherwise I have to switch to Matlab.

@mehdibaghaee َAny comment yet? Where t...

mehdibgh 240
March 10 2017
0 0

@mehdibaghaee َAny comment yet?

Where the problem is?

@acer regarding to your comments I still...

mehdibgh 240
March 10 2017
0 0

@acer regarding to your comments I still couldnt calculate the square root of a 111*111 random matrix by Maple with my machine
 

NULL

> 

restart

Error, invalid input: with expects its 1st argument, pname, to be of type {`module`, package}, but received shareman

  
> 

with(LinearAlgebra):
> 

B := RandomMatrix(111, density = .75, generator = 0 .. .5)
> 

``
> 

MatrixFunction(evalf(B), sqrt(v), v):
> 

``
> 

NULL
> 

NULL


 

Download askkk.mw

, Whereas Matlab calculate similar square root of a 111*111 random matrix in 2 seconds.

What to do to calculate it with Maple?

@Kitonum Ok lets do it for my problem ...

mehdibgh 240
March 06 2017
0 4

@Kitonum Ok lets do it for my problem


 

 

> 

restart

Error, invalid input: with expects its 1st argument, pname, to be of type {`module`, package}, but received shareman

  
> 

with(LinearAlgebra):
> 

sys := [-(1/2)*a*b*omega^2*(2*I__0*U[0, 0]+2*I__1*Phi[0, 0])-lambda[1, 0]*K__ux0+b*U[0, 0]*(K__ux0+K[uxa])+2*a*U[0, 0]*(K__uy0+K[uyb])+(1/2)*b*V[0, 0]*(K__vx0+K[vxa])+(1/2)*b*V[0, 0]*(K__wx0+K[wxa])+lambda[6, 0]*K__uxa-lambda[12, 0]*K__uy0+lambda[17, 0]*K__uyb, -(1/2)*a*b*omega^2*(2*I__0*V[0, 0]+2*I__1*Xi[0, 0])-lambda[2, 0]*K__vx0+(1/2)*b*U[0, 0]*(K__vx0+K[vxa])+a*V[0, 0]*(K__vy0+K[vyb])+(1/2)*b*U[0, 0]*(K__wx0+K[wxa])+lambda[7, 0]*K__vxa-lambda[11, 0]*K__vy0+lambda[16, 0]*K__vyb, -a*b*omega^2*I__0*W[0, 0]+b*W[0, 0]*(K__vx0+K[vxa])+a*W[0, 0]*(K__wy0+K[wyb])-lambda[5, 0]*K__wx0+lambda[10, 0]*K__wxa-lambda[15, 0]*K__wy0+lambda[20, 0]*K__wyb, (1/2)*a*b*(2*A[4, 5]*Xi[0, 0]+2*A[5, 5]*Phi[0, 0]+2*Phi[0, 0]*R[3, 3])-(1/2)*a*b*omega^2*(2*I__1*U[0, 0]+2*I__2*Phi[0, 0])+lambda[25, 0]*E[2, 5]+lambda[26, 0]*F[2, 5]+lambda[27, 0]*E[2, 5]+lambda[28, 0]*F[2, 5]+b*Phi[0, 0]*(`K__φx0`+K[`φxa`])+a*Phi[0, 0]*(`K__φy0`+K[`φyb`])+lambda[3, 0]*(-`K__φx0`+F[3, 1])+lambda[5, 0]*A[5, 5]+lambda[8, 0]*`K__φxa`+lambda[10, 0]*A[5, 5]-lambda[14, 0]*`K__φy0`+lambda[15, 0]*A[4, 5]+lambda[19, 0]*`K__φyb`+lambda[20, 0]*A[4, 5]+lambda[21, 0]*E[1, 5]+lambda[22, 0]*F[1, 5]+lambda[23, 0]*E[1, 5]+lambda[24, 0]*F[1, 5], (1/2)*a*b*(2*A[4, 4]*Xi[0, 0]+2*A[4, 5]*Phi[0, 0])-(1/2)*a*b*omega^2*(2*I__1*V[0, 0]+2*I__2*Xi[0, 0])+lambda[25, 0]*E[2, 4]+lambda[26, 0]*F[2, 4]+lambda[27, 0]*E[2, 4]+lambda[28, 0]*F[2, 4]+b*Xi[0, 0]*(`K__ψx0`+K[`ψxa`])+a*Xi[0, 0]*(`K__ψy0`+K[`ψyb`])-lambda[4, 0]*`K__ψx0`+lambda[5, 0]*A[4, 5]+lambda[9, 0]*`K__ψxa`+lambda[10, 0]*A[4, 5]+lambda[13, 0]*`K__ψy0`+lambda[15, 0]*A[4, 4]+lambda[18, 0]*`K__ψyb`+lambda[20, 0]*A[4, 4]+lambda[21, 0]*E[1, 4]+lambda[22, 0]*F[1, 4]+lambda[23, 0]*E[1, 4]+lambda[24, 0]*F[1, 4], 0, E[3, 1]*lambda[1, 0]+E[3, 1]*lambda[6, 0]+E[3, 2]*lambda[11, 0]+E[3, 2]*lambda[16, 0]+E[3, 6]*lambda[2, 0]+E[3, 6]*lambda[7, 0]+E[3, 6]*lambda[12, 0]+E[3, 6]*lambda[17, 0]+F[3, 1]*lambda[8, 0]+F[3, 2]*lambda[13, 0]+F[3, 2]*lambda[18, 0]+F[3, 6]*lambda[4, 0]+F[3, 6]*lambda[9, 0]+F[3, 6]*lambda[14, 0]+F[3, 6]*lambda[19, 0]]:
> 

var := [U[0, 0], V[0, 0], W[0, 0], Phi[0, 0], Xi[0, 0], Theta[0, 0], Omega[0, 0], lambda[1, 0], lambda[2, 0], lambda[3, 0], lambda[4, 0], lambda[5, 0], lambda[6, 0], lambda[7, 0], lambda[8, 0], lambda[9, 0], lambda[10, 0], lambda[21, 0], lambda[22, 0], lambda[23, 0], lambda[24, 0], lambda[11, 0], lambda[12, 0], lambda[13, 0], lambda[14, 0], lambda[15, 0], lambda[16, 0], lambda[17, 0], lambda[18, 0], lambda[19, 0], lambda[20, 0], lambda[25, 0], lambda[26, 0], lambda[27, 0], lambda[28, 0]]:
> 

A, b := GenerateMatrix(sys, var)

Error, recursive assignment

  
> 

``
> 

``

``

``


 

Download soal.mw

Where the problem is?

 

@vv Many Thanks,...

mehdibgh 240
March 05 2017
0 0

@vv Many Thanks,

5 6 7 8 9 Page 7 of 9