Re: [deal.II] type/value mismatch error in geometry/geometries/point_xy

Hi Vimukthi, I did end up by resolving that issue by making sure that gcc, g++, and gfortran were all on a single version. Ended up using 9 for them. Seems to be working now. Michael On Wed, Feb 21, 2024 at 8:06 AM Wolfgang Bangerth wrote: > On 2/21/24 02:46, vimukthi Nanayakkara wrote: &g

Re: [deal.II] type/value mismatch error in geometry/geometries/point_xy

Hi Wolfgang, I appreciate you taking a look at it. Before changing the files I'll try with other compilers, currently, I'm using gcc11.4, the latest available on Ubuntu 22.04, so maybe others will work. Thanks again, Michael On Sat, Feb 10, 2024 at 1:21 PM Wolfgang Bangerth wrote: > > M

[deal.II] Re: type/value mismatch error in geometry/geometries/point_xy

.5.1/include/deal.II/bundled/boost/geometry/geometries/point_xy.hpp:78:27: error: type/value mismatch at argument 1 in template parameter list for ‘template class std::set’ 78 | { this->template set<1>(v); } On Friday, February 9, 2024 at 9:58:04 AM UTC-5 Michael Pilipchuk wrote:

[deal.II] type/value mismatch error in geometry/geometries/point_xy

Hello all, I installed deal.II-v9.5.1 using candi and that installed cleanly. Then I tried to install the open-source PRISMS-Plasticity package that uses deall.II, and it compiled with no errors up until using "make release". The chain of errors

RE: [deal.II] How to get the shape_value at the quadrature point in the reference cell?

rad(i, center) << std::endl;  } I wonder if there are other ways to access the shape_grad on the ref cell directly. Do I need to use fe_values.inverse_jacobian(q_index) to convert from shape_grad on the real cell to the reference cell?  Best,Michael From: Wolfgang BangerthSent: Sunday, Septem

RE: [deal.II] Re: How to get the shape_value at the quadrature point in the reference cell?

Hi Dear Thang,Thank you for you code snippet! It is simple and works well and is what I want. The link helps me understand the shape function in deal.ii. Have a good weekend!Michael  From: Thang W PhamSent: Saturday, September 11, 2021 12:20 AMTo: deal.II User GroupSubject: [deal.II] Re: How

RE: [deal.II] How to watch the variable of BlockVector type?

Yeah! .block(0) is essential. Without this, the data obtained is garbage! Thanks for this great hint! From: Wolfgang BangerthSent: Wednesday, August 18, 2021 8:37 PMTo: dealii@googlegroups.comSubject: Re: [deal.II] How to watch the variable of BlockVector type? On 8/18/21 7:18 PM, Michael Lee

RE: [deal.II] How to watch the variable of BlockVector type?

Sorry, it seems to work using *((double(*)[20])residual.begin()).   From: Michael LeeSent: Wednesday, August 18, 2021 7:18 PMTo: deal.II User GroupSubject: [deal.II] How to watch the variable of BlockVector type? I want to check if I have system_rhs assembled correctly. Is there an easy way

[deal.II] How to watch the variable of BlockVector type?

does not work for BlockVector. I'm using VS code. Best, Michael -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "dea

RE: [deal.II] Computing the solution gradient at the quadrature point on a face

into it later. For now I hope I can have a better understanding this topic and implement it in deal.II by hand then I will be comfortable with AD. Thank you!Michael  From: Jean-Paul PelteretSent: Sunday, August 8, 2021 12:16 PMTo: dealii@googlegroups.comSubject: Re: [deal.II] Computing the solution

Re: [deal.II] Debug multithread code

you compiling in debug mode? > >> On Thu, Aug 5, 2021, 18:34 Michael Li wrote: >> I put MultithreadInfo::set_thread_limit(1) in main function, but still >> cannot do step debug in the function >> >> assemble_neumann_contribution_one_cell. Am I doing someth

RE: [deal.II] Debug multithread code

multithread code Take a look at set_thread_limit:https://www.dealii.org/current/doxygen/deal.II/classMultithreadInfo.htmlAs you can see, you can also control it using an environment variable, so you don't have to recompile. On Thu, Aug 5, 2021, 13:36 Michael Lee <lianxi...@gmail.com> wrote:Hi, I fin

[deal.II] Debug multithread code

series one? Regards, Michael -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubs

RE: [deal.II] Computing the solution gradient at the quadrature point on a face

totally right, it converges only if a very small load step is used. The residual grows up fast and the result is meaningless if the load increment is large. Thanks for your great help! Michael    From: Jean-Paul PelteretSent: Wednesday, August 4, 2021 12:19 AMTo: dealii@googlegroups.comSubject: Re

RE: [deal.II] Computing the solution gradient at the quadrature point on a face

} Can you tell me what’s wrong with the code above? Thank you!Michael   From: Jean-Paul PelteretSent: Sunday, August 1, 2021 1:46 PMTo: dealii@googlegroups.comSubject: Re: [deal.II] Computing the solution gradient at the quadrature point on a face Hi Michael, The problem here is that “scratch.solution

[deal.II] Re: Computing the solution gradient at the quadrature point on a face

*get_function_gradients*(solution, >> solution_grads); >> >> for (unsigned int f_q_point = 0; f_q_point < n_q_points_f; ++f_q_point) >> { >>*const Tensor<2,dim,NumberType> _u = >> solution_grads[f_q_point];* >> *...* >>

[deal.II] Computing the solution gradient at the quadrature point on a face

< n_q_points_f; ++f_q_point) { *const Tensor<2,dim,NumberType> _u = solution_grads[f_q_point];* *...* } Can anybody give a hint on that? Thanks, Michael -- The deal.II project is located at http://www.dealii.org/ For mailing l

RE: [deal.II] Matrix Multiplication

nt::push_forward(S, F);} I just undefine the variable to circumvent the Sacado error.  #define ENABLE_SACADO_FORMULATION   Thanks,Michael   From: Jean-Paul PelteretSent: Monday, July 26, 2021 12:46 AMTo: dealii@googlegroups.comSubject: Re: [deal.II] Matrix Multiplication Hi Michael, I revisited t

RE: [deal.II] Matrix Multiplication

. Vienant Kirchhoff material. Is that correct? I compiled cook_membrane.cc a couple of days ago, but for no reason, it fails giving the following errors: /home/michael/dealii-9.2.0/examples/code-gallery/Quasi_static_Finite_strain_Compressible_Elasticity/cook_membrane.cc:1815:5:   required from here/home

[deal.II] Matrix Multiplication

elegant or easier way to do the matrix construction as I need to do other similar operations. Best, Michael -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you

RE: [deal.II] Questions on Step-18

of 180°. So the code I suggest will look like as follows.   const Point<3> curl = …;   const double mag_curl = std::sqrt(curl * curl);   const double  angle = 0.5 * mag_curl; Thank you!Michael  From: Jean-Paul PelteretSent: Monday, July 5, 2021 11:50 PMTo: dealii@googlegroups.comSubje

Re: [deal.II] Questions on Step-18

* curl); const double angle = std::atan(tan_angle); Any comments will be of great help. Best, Michael On Mon, Jun 28, 2021 at 2:21 AM Andrew McBride wrote: > My view here is that the problem is quasi-static and hence time serves to > order events. Hence a displacement inc

RE: [deal.II] Questions on Step-18

the rotation angle (curl of displacement). Do correct me if I’m wrong. Now I’m excited to learn how to make my first patch. - Michael From: Jean-Paul PelteretSent: Saturday, June 26, 2021 2:57 PMTo: dealii@googlegroups.comSubject: Re: [deal.II] Questions on Step-18 Hi Michael, What’s written

Re: [deal.II] Questions on Step-18

, displacement_increment_grads); Should we use the *velocity *([image: \vec{v} = \frac{d \vec{u}}{dt}]), since [image: \vec{\omega} = \frac{1}{2}\nabla \times \vec{v} = \frac{1}{2}\nabla \times \frac{d\vec{u}}{dt}] ? Can you check this as well? Best, Michael On Sat, Jun 26, 2021 at 1:48 AM Jean-Paul

RE: [deal.II] Questions on Step-18

the code roll. Thank you again for your wonderful help! Best,Michael  From: Jean-Paul PelteretSent: Thursday, June 24, 2021 11:43 PMTo: dealii@googlegroups.comSubject: Re: [deal.II] Questions on Step-18 HI Michael. To be honest, I have a hard time to understand the three-filed formulation in Step 44

RE: [deal.II] Questions on Step-18

deformation model. I believe that will make things simpler and help me understand more complex nonlinear models. Meanwhile, I try to compare the results with a reference which only consider the geometrical nonlinearity. -Michael   From: Andrew McBrideSent: Thursday, June 24, 2021 2:23 PMTo: deal.II User

RE: [deal.II] Different get_function_gradients

Hi Jean-Paul,Thank you so much for your explanation. That totally makes sense. These different functions exist for implementing  different weak forms easily. Though I can get everything using (1), but I would prefer using the others when I try to extract components of the gradient.  Best,Michael

[deal.II] Different get_function_gradients

After reading the documentation on Handling vector valued problems (The deal.II Library: Handling vector valued problems (dealii.org) ), I have the following question. What are the differences between the three

[deal.II] Questions on Step-18

grad[i][1] * grad[j][1] + * * grad[i][2] * grad[j][2]*) / 2; return strain; } Thanks, Michael -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this me

RE: [deal.II] Re: Integrated material and spatial traction forces on boundary not equal

4:07 PM, Michael Lee wrote:> > Is there a way to let deal.II know Trillinos is installed? Do I have to > reinstall deal.ii?> Or, can I compile the code without using Trillinos with some modification? No, you have to re-install deal.II so that at compile time, it knows that Trilinos exi

Re: [deal.II] Re: Integrated material and spatial traction forces on boundary not equal

to reinstall deal.ii? Or, can I compile the code without using Trillinos with some modification? Thanks! Michael On Friday, June 11, 2021 at 4:04:20 AM UTC-6 alexanderc...@gmail.com wrote: > Hi Michael, > > I think that is the first point in my code that something from the > automatic dif

RE: [deal.II] How to define boundary conditions on edges and vertices in 3D?

ges with different boundary_ids for other use. There are still many scenarios that we want to fix some points in the domain. I wonder if we can make use of the AffineConstraints class to deal with Dirichlet boundary conditions though it is not possible for Neumman BCs. Best,Michael From: Wolfg

RE: [deal.II] Re: Integrated material and spatial traction forces on boundary not equal

ite element FE_Q(2). I remember the deflection of beam is a cubic function of coordinate. You may try that see if it improves. Best,Michael From: Alex CumberworthSent: Wednesday, June 9, 2021 9:54 AMTo: deal.II User GroupSubject: [deal.II] Re: Integrated material and spatial traction forces on bound

[deal.II] Re: Integrated material and spatial traction forces on boundary not equal

Hi Alex, I'm learning deal.ii and trying do the similar verification. If it is possible for you to share the code with me? Thank you! Michael On Tuesday, May 11, 2021 at 4:46:55 AM UTC-6 alexanderc...@gmail.com wrote: > Hello, > > As a test to validate my code, I am solving the

Re: [deal.II] Test 4 failed, maybe different installation of MPI? then how to solve this?

Hi Daniel, Thank you, I used mpicc to compile. Replace NULL by nullptr won't help because compiler can not recognize nullptr neither. On Saturday, October 28, 2017 at 1:21:25 PM UTC-5, Daniel Arndt wrote: > > Michael, > > are actually compiling with a MPI compiler? > Does

Re: [deal.II] Test 4 failed, maybe different installation of MPI? then how to solve this?

> > On Thu, Oct 26, 2017 at 6:38 PM, Michael <micha...@gmail.com > > wrote: > > Hi, > > > > Test 4 failed after make test. I make my research and find the following > > answer: > > > > > https://urldefense.proofpoint.

[deal.II] Test 4 failed after installation, no clue

Hi, Test 4 failed after make test. I make my research and find the following answer: https://groups.google.com/d/msg/dealii/mV-juFpnybA/DSBP4DArCQAJ And I get the following suggested commands with returns blow:

Re: [deal.II] Re: Installation error, please help, thx!

UTC-5, Wolfgang Bangerth wrote: > > On 10/25/2017 06:43 PM, Michael wrote: > > > /home/superman/software/dealii-8.5.0/include/deal.II/lac/chunk_sparse_matrix.templates.h: > > In instantiation of ‘somenumber > dealii::ChunkSparseMatrix::matrix_norm_square(const

[deal.II] Installation error, please help, thx!

Hi, I've installed dealii before using the default step from README. This time, I am trying to install dealii with trilinos, but getting the following kind of errors during installation: " /usr/include/c++/6/tuple:648:21: error: no type named ‘type’ in ‘struct std::enable_if’

[deal.II] Installation error with c++ lib

Hi, I've installed dealii before using the default step by README. This time, I am trying to install dealii with trilinos, but getting the following kind of errors during installation: " /usr/include/c++/6/tuple:648:21: error: no type named ‘type’ in ‘struct std::enable_if’

Re: [deal.II] Re: Announcing the deal.II Code Gallery

UTC-5, Jean-Paul Pelteret wrote: > > Thanks Wolfgang, a slight permutation of that seemed to work! I'll submit > a PR in a moment. > > Michael, can you tell me if you've built any of the code gallery examples? > I think that this might be the issue. If you have, can you go into

Re: [deal.II] Re: Announcing the deal.II Code Gallery

Thanks! I am glad it wasn't just me!! Mike On Tuesday, March 7, 2017 at 10:26:43 AM UTC-5, Jean-Paul Pelteret wrote: > > Hi Michael, > > I've just tried to build the documentation with the code gallery and have > run into similar problems. I'm going to fiddle around to see if I

Re: [deal.II] Re: Announcing the deal.II Code Gallery

6, 2017 at 10:58:28 PM UTC-5, Wolfgang Bangerth wrote: > > On 03/06/2017 08:14 PM, Michael Harmon wrote: > > The doxygen.log file seems to get deleted... but here's the copy of it I > took > > right before the error gets thrown and it is deleted. > > Hm, yes, that

Re: [deal.II] Re: Announcing the deal.II Code Gallery

, March 6, 2017 at 6:03:47 PM UTC-5, Wolfgang Bangerth wrote: > > On 03/06/2017 03:59 PM, Michael Harmon wrote: > > down and I get the error attached as a picture... any ideas on what's > > going wrong? > > Not sure -- can you look into doxygen.log see if there are errors? M

[dealii-developers] Re: Announcing the deal.II Code Gallery

y. >> >> Best >> W. > > > On Tuesday, February 28, 2017 at 4:24:31 PM UTC+1, Michael Harmon wrote: >> >> Hi, I wanted to submit a example code to the code gallery, but was >> wondering is there a tool to see how the page will look before when all

[deal.II] Re: Announcing the deal.II Code Gallery

irectory and build the deal.II documentation as described in >> the >> readme. It will pick up the code gallery and create joint documentation >> for >> the tutorial and the code gallery. >> >> Best >> W. > > > On Tuesday, February 28, 2017

[dealii-developers] Re: Announcing the deal.II Code Gallery

Hi, I wanted to submit a example code to the code gallery, but was wondering is there a tool to see how the page will look before when all the doxygen comments are taken into account.. like the following: https://dealii.org/developer/doxygen/deal.II/code_gallery_cdr.html On Friday, February

[deal.II] Re: Announcing the deal.II Code Gallery

Hi, I wanted to submit a example code to the code gallery, but was wondering is there a tool to see how the page will look before when all the doxygen comments are taken into account.. like the following: https://dealii.org/developer/doxygen/deal.II/code_gallery_cdr.html On Friday, February

[deal.II] Re: possible memory leak using TrilinosWrappers direct solver?

Is your matrix updating at every time step or just the right hand side vector? If its just the right hand side vector is changing, the development version of dealii has the ability to factorize the matrix once and then just do solve(x,b) at every time step. I'm not sure, but it Might help. On

[deal.II] cout stiffness matrix in step-3

Hello! I'm trying to cout the stiffness matrix $A_{ij}$ in step-3 for a reasonably small domain 2d mesh..is $A_ij$ a function in the FEValues class..? Just trying to learn the basics of dealii.. Thanks in advance, Michael -- The deal.II project is located at http://www.dealii.org

[deal.II] deal.ii workflow question?

W What is a good way to trace/compile parts of deal.ii code without doing cmake/make/makerun? In particular I want to cout matrices to see what their entries are 9:26 PM (7 minutes ago) -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see

[deal.II] deal.ii workflow?

What is a good way to trace/compile parts of deal.ii code without doing cmake/make/makerun? Just an enthusiastic student wanting to really write some deal.ii code! Thanks! Mike -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see

Re: [deal.II] Re: question on doxygen

Just on my own program. On Wednesday, November 16, 2016 at 9:44:59 AM UTC-5, Wolfgang Bangerth wrote: > > On 11/16/2016 07:41 AM, Michael Harmon wrote: > > > > I had to change the http from: > > > > MATHJAX_REPATH = http://cdn.mathjax.org/math

[deal.II] Re: question on doxygen

Thanks Timo! I had to change the http from: MATHJAX_REPATH = http://cdn.mathjax.org/mathjax/latest/ to https: MATHJAX_REPATH = http*s*://cdn.mathjax.org/mathjax/latest/ On Tuesday, November 15, 2016 at 12:03:55 PM UTC-5, Michael Harmon wrote: > > Hi, > > I used deal for my phd

[deal.II] Re: Publications based on deal.II

in photoelectrochemical (PEC) solar cells}, Author = {Michael Harmon}, School = {University of Texas at Austin}, Year = {2016} } On Sunday, January 17, 2016 at 6:09:44 PM UTC-5, Wolfgang Bangerth wrote: > > > All, > > as you m

[deal.II] question on doxygen

Hi, I used deal for my phd work, loved using it and wanted to put my code up on github. I spent a lot of time documenting it with doxygen, but latex equations only work on my local machine and not on github: https://mdh266.github.io/PECS/index.html Anyone else have this issue? Thanks!

[deal.II] Re: Using the solution from one problem as a boundary condition in another problem with matching mesh on the boundary

deal ii. > > On Tuesday, August 23, 2016 at 5:55:32 PM UTC+3, Michael Harmon wrote: >> >> Hey, >> >> I had a similar problem: PDES in separate domains that are coupled >> through an interface as a boundary condition. You can go about it using >> one

[deal.II] Re: Using the solution from one problem as a boundary condition in another problem with matching mesh on the boundary

Hey, I had a similar problem: PDES in separate domains that are coupled through an interface as a boundary condition. You can go about it using one triangulation; I attempted to do this at first, but ended up using multiple meshes. The fact you have matching meshes on the boundary is good.