[deal.II] Re: writing values in the solution vector for constrained dofs

2021-06-23 Thread Bruno Turcksin
Simon, Yes, you should be able to skip the constrained hanging nodes. Like you said the value will be overwritten when you call distribute(). Let us know if you get something unexpected but that should work. Best, Bruno On Wednesday, June 23, 2021 at 5:00:06 AM UTC-4 Simon wrote: > Dear

[deal.II] Deal.II compilation error issued at example steps phase

2021-06-23 Thread Debbahi Saad
Dear Deal.ii community, I'm trying to compile the Deal.II (9.2.0, 9.3.0 or dev branch) with examples enabled, the compilation goes well until it reaches the example phase, for example step-9 in version dev branch breaks the compilation with an error *../../lib/libdeal_II.g.so.10.0.0-pre:

Re: [deal.II] 5x5 Matrix Inversion Problem

2021-06-23 Thread Wells, David
Hi Lucas, C++ doesn't allow us to overload operator[] with multiple arguments so that indexing has to be with single values - so yes, it's a language thing. This matches the C usage of square brackets where they only accept one value at a time. The state attribute is automatically set by

Re: [deal.II] Robust way of setting an affine constraint on a single point

2021-06-23 Thread blais...@gmail.com
Thank you Wolfgang, this is exactly what I needed :) Best Bruno On Tuesday, June 22, 2021 at 6:04:42 p.m. UTC-4 Wolfgang Bangerth wrote: > On 6/22/21 3:12 PM, blais...@gmail.com wrote: > > > > This is most likely a very simple question, but what would be a good way > to > > robustly set an

[deal.II] writing values in the solution vector for constrained dofs

2021-06-23 Thread Simon
Dear all, I am dealing with adaptivity, so hanging nodes are present in my mesh which are handled by a constraints object like in the following: AffineConstraints...constraints; make_hanging_node_constraints...(); constraints.close(); My solution vector is filled via a local approach. Basically

Re: [deal.II] Outputting DataOut Object into Personal Object Rather Than Outputing to File

2021-06-23 Thread Wolfgang Bangerth
On 6/23/21 1:57 AM, Chen R wrote: I am mainly trying to integrate the dealii system with our own geometric modeling program to do FE calculations. Right now created a basic elasticity solver based on the tutorials. The thing is I want to create a triangulation from our geometry and calculate

[deal.II] Marking Boundary Nodes in Triangulation

2021-06-23 Thread Chen R
Hi all, So I am trying to use dealii together with our own geometric modeling environment to do different FE calculations on. I am transferring our geometry into the Triangulation object based on the way that the read_msh function does. The problem is that I want to keep track of boundary

Re: [deal.II] Outputting DataOut Object into Personal Object Rather Than Outputing to File

2021-06-23 Thread Chen R
Hi David, I am mainly trying to integrate the dealii system with our own geometric modeling program to do FE calculations. Right now created a basic elasticity solver based on the tutorials. The thing is I want to create a triangulation from our geometry and calculate on it then present the

Re: [deal.II] dealii 9.3.0 make install fails at "Generating mpi.inst" with Invalid instantiation list: missing 'for'

2021-06-23 Thread vachan potluri
I had noticed that this make error would occur if the file being expanded doesn't have a prefix 'for'. The following snippet is from line 453- of dealii/cmake/scripts/expand_instantiations.cc if (!has_prefix(whole_file, "for")) { std::cerr << "Invalid instantiation list: missing 'for'" <<