Re: [deal.II] Quasi Static Compressible Finit Strain for Heterogenous Hyperelastic material

2020-10-26 Thread Jean-Paul Pelteret
Dear Animesh, Let me preface my suggestion by saying that you can introduce some heterogeneity in more than one way. That is, you're not bound to using a material ID (which is what I'm going to suggest), but you could use some geometric arguments to determine which class or material

Re: [deal.II] Linear Operator

2020-10-26 Thread Nikki Holtzer
Matthias, When you helped me implement the operator N_c which is added to my linear_operator , i.e. auto N_c = temp1_coef*M + temp_coef*Kronecker; auto system_matrix = linear_operator(linear_part) + N_c; the first term is not complete. Really I want, auto N_c = temp1_coef*|Psi^m|^2 *M +

Re: [deal.II] Linear Operator

2020-10-26 Thread Matthias Maier
Dear Nikki, A "linear operator" in deal.II is an object that linearly transforms a vector into another one, i.e., it has a vector space as domain of definition and a vector space as range space. In your case M.matrix_scalar_product(u,u) is (u^T M u), this is a quadratic operation, not a linear

[deal.II] Linear Operator

2020-10-26 Thread Nikki Holtzer
Hello all, I am struggling to turn a matrix scalar product into a linear operator. I currently have: temp = M.matrix_scalar_product(u,u); and I need to turn 'temp' into a linear operator. I have tried turning 'M' into a linear operator first and then performing the multiplication and

Re: [deal.II] dealii installation issue on NERSC Cori

2020-10-26 Thread Timo Heister
>Thanks for you suggestion. I tested step-1 and step-2 again and they ran > successfully after I set the number of OpenMP threads to 1. I only plan plan > on running MPI codes without any multi-threading. The BLAS implementation I > am using is a part of the Intel MKL library which is

Re: [deal.II] dealii installation issue on NERSC Cori

2020-10-26 Thread Aaditya Lakshmanan
Hi Timo, Thanks for you suggestion. I tested step-1 and step-2 again and they ran successfully after I set the number of OpenMP threads to 1. I only plan plan on running MPI codes without any multi-threading. The BLAS implementation I am using is a part of the Intel MKL library which is

Re: [deal.II] dealii installation issue on NERSC Cori

2020-10-26 Thread Timo Heister
Can you try setting export OMP_NUM_THREADS=1 before running any of the examples? If you are planning on running MPI-only codes, you don't need to use OpenMP in your BLAS implementation. On Sun, Oct 25, 2020 at 7:25 PM Aaditya Lakshmanan wrote: > > Hi Timo, >I cleaned the previous

[deal.II] Quasi Static Compressible Finit Strain for Heterogenous Hyperelastic material

2020-10-26 Thread Animesh Rastogi IIT Gandhinagar
Hello All, I am trying to solve the wrinkling problem of a thin stiff film attached to a soft substrate in dealii. I am using the code - from the code gallery for this purpose. For the simulations, I would have to consider the film and the substrate with different material properties. I