Greetings

I trust that everyone is well here. I have an inquiry regarding the partitioning of a matrix BA, defined on both regular and irregular domains, into three matrices: Ad, Ab, and Ae, such that NewBA = Ad + Ab + Ae. Here, Ad comprises the entries of BA that reside within the domain, Ab includes the elements of BA located on the boundary, and Ae is derived as BA - Ad - Ab for any values of Mx and My. The specifics of matrix BA are contained in the attached document where NewBA constructed manually for different values of Mx and My for better understanding.
Splitting_a_matrix.mw
The attached file contains the matrix BA constructed in a square format. How may the BA matrix be adapted to create an irregular shape, such as a quarter circle?

The red dots indicate the mesh within the domain Ad, while the meshing along the blue line should occur in Ab.

I await your kind reply. Kindly ensure your well-being

Hi

I hope everyone is fine and doing well. I want to constrcut the set of monomials {p[0],p[1],...,p[m-1]} for any value of m for example, for m=6 the monoial is define as:

p[0]:=1;

p[1]:=x;

p[2]:=y;

p[3]:=x^2;

p[4]:=x*y;

p[5]:=y^2;

and similarly for m=10 the monomials should be given as:

p[0]:=1;

p[1]:=x;

p[2]:=y;

p[3]:=x^2;

p[4]:=x*y;

p[5]:=y^2;

p[6]:=x^3;

p[7]:=x^2*y;

p[8]:=x*y^2;

p[9]:=y^3;

I am waiting for your positive response. Please take care and thanks

Hi Dear,

I hope everyone is fine here. In the attached file, I have generated a square matrix "Q" using two-dimensional polynomials. The dimension of the square matrix "Q" depends on M1 and M2 parameters. In my simulation, sometimes I need this matrix of 1000 by 1000 dimensions. Using the attached method, it took a lot of time to compute two-dimensional polynomials and then to compute the general square matrix "Q." I wanted to write this matrix using proc (procedures). Maybe by using this way, I don't need to compute the polynomials, and it took less time to compute the square matrix "Q." I know how to generate a matrix using proc when its dimension depends on one parameter. However, here, the dimension of matrix "Q" depends on two parameters, M1 and M2. So, I am a little bit confused about how to adjust them in proc. Please see the attached file and share your useful ideas. 

help.mw

Thanks in advance

Hi
I hope you are doing well. I have plotted (in the attached file) the contour plot of the function and its density plot; both have the same behavior but different appearances (error in direction may be rotation needs to apply). I don't know why it happens because this code works well for other solutions. Kindly have a look and fix the issue. I shall be waiting for your positive response. Please take care.
Help.mw

Hi,

I hope everyone is fine here. I want to write the following nested loop in my code

restart; TOL := 10^(-3); v[i] := 10^(-4); u[i] := 10^(-4);
if abs(v[i]-1)<=TOL then    if abs(u[i])<=TOL then break    else t[i+1]:=t[i]+alpha; else; fi; fi;

it is not working. If the first condition (v[i]-1)<=TOL is fulfilled, then we have to check the second condition, u[i]<=TOL. That’s why I don’t use “and” here. If both conditions are fulfilled, then stop the loop and if second the condition not fulfilled (but first fulfill) then t[j+1]=t[j]+alpha. 

I shall be waiting for your positive response.