About operations and tmatrix class

[jgenee]

Hello everyone,

I am a relatively new user of Mfront, so the question might be relatively easy to answer. Though, I was wondering if operations between nonsquare matrix (NxM with N and M different) are handled within Mfront (multiplication, etc).

For example, I defined two matrices A and B as :

@localvariable tfel::math::tmatrix<9,18,real> A;
@localvariable tfel::math::tmatrix<18,9,real> B;

I need to calculate something like B = invert(transpose(A)*A) * transpose(A).
I tried this the hard way but it does not seem to work, as there is no match for " operator * ". I guess the problem will be the same for the invert function..

Am I missing something ?

Thanks in advance !

[thelfer]

@jgenee, Thanks for your interest in TFEL/MFront. You did not miss anything. It just seems that I never needed to implement the matrix product. In pratice, I have implemented the product of tensors. Maybe you could give me more insight on what your are trying to do just to figure out why you need it. The standard implementation is easy to add.

[SCrenSC]

Thank you! The following lines are working.
auto success = bool{};
tie(success, B) = tfel::math::extensions::my_complex_operation(...);
if(!success){
return false; // not sure
}

[thelfer]

Great !

[thelfer]

Let me know if you have other troubles

[SCrenSC]

Some new troubles indeed.
The compilation was successful, however when I tried the behaviours in Amitex, I got:
The behaviour 'mymatbehavior' has thrown a generic standard exception : initialization of behaviour 'mymatbehavior' failed.
It seems it is not recognized in Amitex.

[thelfer]

You may want try:

auto success = bool{};
tie(success, B) = tfel::math::extensions::my_complex_operation(...);
if(!success){
  std::cout << "my_complex_operation failed\n";  
  return false; // not sure
}
std::cout << "my_complex_operation was successful\n";

[jgenee]

Thanks for your answer.

I am implementing a non local behavior law using MFront and Amitex. To that end, I defined a 9x18 matrix (A in the previous message) that contain information about dislocations types. Then, the matrix is used in the calculation of the 18x18 matrix B that relates two state variables in the behavior law (one being an non-local value calculated outside of the behavior law by Amitex). I defined the calculation of B in a second mfront file that I imported into the main one containing the law.

Basically, A is initialized as a 9x18 matrix. And I need to compute B = invert( transpose(A) * A ) * transpose(A).
As you told me, operation * is not defined for matrices, same for invert and transpose functions (not sure about invert though, as it is here used on a square matrix).

Hope I am clear enough !