« first day    last day (16 days later) » 

00:07
Be my guest =)
We suppose that $\Omega \subset \mathbb{R}^2$ is an open, bounded region of a smooth boundary $\Gamma=\partial{\Omega}$, $f \in L^2(\Omega), k \in L^{\infty}(\Omega)$ and there is an $a >0$ such that $k(x) \geq a, \forall x \in \Omega$.
We want to find a solution of the elliptic problem

$\left\{\begin{matrix}
k(x)\Delta u(x)=f(x) & \text{ in } \Omega\\
u=0 & \text{in } \Gamma
\end{matrix}\right.$

We consider that the region $\Omega$ is seperated into two subregions $\Omega=\Omega_1 \cup \Omega_2$ and the function $k(x)$ is constant in each region
I have thought the following:
Let h_1 be the partition for the interval [-1,0] and h_2 the partition for the interval [0,1].
At the y-interval [-1,0] k is equal to k_2.
So, since -1+jh_1=0 => j=1/h_1=N_1 at the matrix A for j=0,...,N_1 , k is equal to k_2.
At the y-interval [0,1] ,k is equal to k_1.
So, since jh_2=1 => j=1/h_2=N_2 at the matrix A for j=N_1+1,...,N_2, k is equal to k_1.
Is this right? @ekkilop
I tried to implement it in matlab, by starting to put B in A as follows:
But I get the following warning message:
>> name(2,1, 4, 5)
??? Error using ==> plus
Matrix dimensions must agree.

Error in ==> name at 35
A=diag(v)+diag(w,-1)+diag(u,1)
Do you maybe have an idea where my mistake it? @ekkilop
Tell me, in your code, are the input parameters N and k scalars?
Yes, they are... @ekkilop
00:23
So the problem is in the function name?
@ekkilop Yes
In the code, name is located in the same m-file as the function B, but the function B never calls name. How is name supposed to be used?
@ekkilop Thay are not located in the same m-file... I just wrote both in the same pastebin...
ah, I see
well give me a minute and I'll have a look
Ok :) @ekkilop
00:28
So, the vector v is supposed to be one element longer than the vectors u and w but currently it is two elements longer so somewhere there is a dimensional error
Yes, I thought so too but I can't find the error... Do you maybe have an idea where it is? @ekkilop
It is hard for me to tell exactly what the error is since I don't really understand what the code does. However, from my experience, these types of problems usually arise when you accidentally represent one point in your grid twice.
You have an interface in your domain, Sigma. My guess is that you represent the points along the interface twice; both in your discretisation of Omega_1 and Omega_2. I suspect this is what introduces the inconsistent dimensions in your problem
I hope that helps you. It's getting pretty late where I am so I'm tucking in. Let me know how things proceed =) Good night!
@ekkilop Ok.. I will look at my code and I will tell you tomorrow if I found my mistake... Good night!!! :-)
 
16 hours later…
16:41
@ekkilop Hello
 
4 hours later…
21:10
Hello @MichaelMedvinsky
hi
it is hard to me to read this chat, particularly your problem...do I need to read it or your question is general?
I will tell you my question.. You don't have to read it all.. @MichaelMedvinsky
your question is the the one I see on the last image?
or you have another one to ask?
21:26
@MichaelMedvinsky Yes, that's my question.
why do you use a chat for this?
I have always asked about this subject at the main but I didn't get an answer and I saw that it is your field... :)
Do you have an idea?
do you have a link to your question?
I haven't ask this one...others from this subject... @MichaelMedvinsky
and you never got an answer?
or it wasn't sufficient?
21:32
No answer... For example here: math.stackexchange.com/questions/1491211/…
commented on that one, also please create a question for the question here - I will try to write answer, but not immediately - give me some time, shouldn't be longer then week. also send me links to other unanswered questions.
btw, did you try to post a bounty on your questions?
may I ask about your background?
21:48
Ok, I will think about it again and if I still face difficulties I will ask it at the main. Because it's the first question and I need it for the other ones.
I study applied mathematics.
In Greece @MichaelMedvinsky
You? @MichaelMedvinsky
not too far from you right now - israel
Nice @MichaelMedvinsky
Are you an undergraduate ? @MichaelMedvinsky
I'm teaching undergraduate students:)
22:01
Are you a professor? @MichaelMedvinsky
something like this
Great!!!
Do you like your job? @MichaelMedvinsky
@MichaelMedvinsky Nice...
How many years do you teach by now? @MichaelMedvinsky
Hi @rych !!!

« first day    last day (16 days later) »