Currently, Zhiwen is solving,
\partial_i \Gamma_{ij} \partial_j u = \alpha
This is a scalar equation, not a vector or matrix equation.
On Tue, Dec 16, 2008 at 2:12 PM, R. Edwin Garcia <[email protected]> wrote:
> Hi Dan,
>
> While the equation that Zhiwen sent you might not make sense, her intent is
> to solve for the rotational of a rank two tensor object. Something like:
>
> \nabla \times \beta = \alpha
>
> where \beta and \alpha are three by three matrices.
>
> You have a good point, though. Zhiwen, could you check that what you are
> trying to solve actually is what you intended?
>
> Cheers,
> EG
>
> On Dec 16, 2008, at 1:55 PM, Daniel Wheeler wrote:
>
>>
>> On Tue, Dec 16, 2008 at 11:27 AM, Zhiwen Liang <[email protected]>
>> wrote:
>>
>>> So every element in the mesh should have diffusion coefficient like this:
>>> D=[0 0 1
>>> 0 0 0
>>> -1 0 0]
>>> Does this also lead to 0*x=b?
>>
>> I think so. The equation doesn't make sense,
>>
>> \partial_y \partial_z u - \partial_z \partial_y u
>>
>> is trivially zero. Right or am I losing my mind? So you are trying to
>> solve "0 = \alpha".
>>
>>>>> Hope all is going well.
>>>>> Recently, I am trying to solve for an equation like this:
>>>>> \frac{\partial^2 u}{\partial y \partial z}-\frac{\partial^2 u}{\partial
>>>>> z
>>>>> \partial y}=\alpha
>>
>> The above is essentially "0 = \alpha". Do you agree?
>>
>> --
>> Daniel Wheeler
>>
>
>
--
Daniel Wheeler
