Re: [deal.II] Efficient Matrix-based Calculation

2021-10-22 Thread Jean-Paul Pelteret
Dear Masoud, > As you mentioned, we expect this to be slower than the equivalent code in > step-8; however, it is very slower than the same code that I wrote before in > Matlab. A quick question: are you sure that you’re running your problem in release mode, rather than debug mode? I had a

Re: [deal.II] Efficient Matrix-based Calculation

2021-10-20 Thread Masoud Ahmadi
Dear Prof. Bangerth, Thanks for your answer. Indeed, I took the declaration of matrices out of the loop; but still, as I mentioned, the part that is substantially reducing the speed is the product between the matrices. As you mentioned, we expect this to be slower than the equivalent code in

Re: [deal.II] Efficient Matrix-based Calculation

2021-10-19 Thread Wolfgang Bangerth
Masoud, Here is a part of the tutorial code for example 8 to calculate stiffness matrix of an element: Screenshot 2021-10-15 at 12.48.39.png And here is my equivalent code for the same purpose with a matrix-based viewpoint: Screenshot 2021-10-15 at 23.12.13.png it is based on the

Re: [deal.II] Efficient Matrix-based Calculation

2021-10-13 Thread Wolfgang Bangerth
On 10/12/21 10:33 AM, Masoud Ahmadi wrote: I'm trying to use a Matrix-based solution for a general 3D FE elastic problem (just like example-8). In that example the authors used a linear scaler solution for calculating the stiffness matrix of each element; but, due to some reasons I want to

[deal.II] Efficient Matrix-based Calculation

2021-10-12 Thread Masoud Ahmadi
Dear all, I'm trying to use a Matrix-based solution for a general 3D FE elastic problem (just like example-8). In that example the authors used a linear scaler solution for calculating the stiffness matrix of each element; but, due to some reasons I want to use a matrix-based solution so that