[deal.II] Output the solution to .h5 file (hdf5)

Hello everyone! For the step-3 hdf5 to output the solution data, I have an error message returned, like the following: [ 66%] Built target step-3 [100%] Run step-3 with Debug configuration Number of active cells: 1024 Number of degrees of freedom: 1089 DEAL:cg::Starting value 0.121094

Re: [deal.II] Re: How to fix the CMake Error: The library at ~/dealii-9.3.2/build was not configured with DEAL_II_WITH_UMFPACK = ON, thanks.

Dr. Bangerth, Thank you so much! Best, Judy On Monday, February 28, 2022 at 1:08:53 PM UTC-5 Wolfgang Bangerth wrote: > On 2/24/22 16:42, Judy Lee wrote: > > I just removed with /var/lib/apt/lists/lock, and installed two libraries > > BLAS and LAPACK. The CMake Error does no

Re: [deal.II] How to fix the error '/home/~/libdeal_II.g.so.9.3.2', thanks.

Dr. Bangerth, The tutorials are compiled through now, with build. Best, Judy On Friday, February 25, 2022 at 9:00:57 PM UTC-5 Wolfgang Bangerth wrote: > On 2/24/22 12:18, Judy Lee wrote: > > > > I have a quick question in running a tutorial file: error shows no rule > t

[deal.II] Re: How to fix the CMake Error: The library at ~/dealii-9.3.2/build was not configured with DEAL_II_WITH_UMFPACK = ON, thanks.

I just removed with /var/lib/apt/lists/lock, and installed two libraries BLAS and LAPACK. The CMake Error does not change, after the two libraries were installed. On Thursday, February 24, 2022 at 4:50:21 PM UTC-5 Judy Lee wrote: > Hello everyone! > > I have a CMake Error at CMakeLis

[deal.II] How to fix the CMake Error: The library at ~/dealii-9.3.2/build was not configured with DEAL_II_WITH_UMFPACK = ON, thanks.

Hello everyone! I have a CMake Error at CMakeLists.txt:27 (MESSAGE): Error! The deal.II library found at /home/~/dealii-9.3.2/build was not configured with DEAL_II_WITH_UMFPACK = ON One or all of these are OFF in your installation but are required for this tutorial step. --

Re: [deal.II] Evaluating an expression u(x), can that take *q_index* information

d weighted. This is wonderful! I will keep up my work. Happy Holidays! Best, Judy On Wednesday, December 22, 2021 at 12:27:59 AM UTC-5 Wolfgang Bangerth wrote: > On 12/21/21 8:59 AM, Judy Lee wrote: > > > > Can I make an expression u(x) for 1D analytical solution, th

Re: [deal.II] Evaluating an expression u(x), can that take *q_index* information

/21 8:59 AM, Judy Lee wrote: > > > > Can I make an expression u(x) for 1D analytical solution, that can take > > *q_index* for evaluating u(@ *q_index*), while skip obtaining *q_index* > > locations/points value? > > What is 'q_index'? It's usually the index (an

[deal.II] Evaluating an expression u(x), can that take *q_index* information

Hello everyone! Can I make an expression u(x) for 1D analytical solution, that can take *q_index* for evaluating u(@ *q_index*), while skip obtaining *q_index* locations/points value? I am doing the residual: *u(x) - u_h(x)* under each cell, in which *u_h(x)* comes from approximation that has

Re: [deal.II] Store data from a Tensor<1, dim> to a Vector

Dr. Bangerth, Thank you so much! Points are printed, I will keep up my work. Best, Judy On Tuesday, December 14, 2021 at 5:49:37 PM UTC-5 Wolfgang Bangerth wrote: > On 12/14/21 08:32, Judy Lee wrote: > > > > I tried *DoFTools::map_dofs_to_support_points(mappi

Re: [deal.II] Store data from a Tensor<1, dim> to a Vector

Hello everyone! I tried *DoFTools::map_dofs_to_support_points(mapping, dof_handler, support_points)*, to obtain coordinates on *local_dof_indices*, not running thru yet. My .cc file is attached with reference to step-3 in 2D, in which I added Lines 137-138 with reference to Lines 126-127, as

Re: [deal.II] Store data from a Tensor<1, dim> to a Vector

Dr. Bangerth, Thank you so much! I will keep up my work. Best, Judy On Wednesday, December 8, 2021 at 6:36:46 PM UTC-5 Wolfgang Bangerth wrote: > > Judy: > > > (1) Is there an internal command/function that outputs > > higher-order-elemental (Q2, Q3 elements) nodal physical coordinate? > >

Re: [deal.II] Store data from a Tensor<1, dim> to a Vector

ed thru finite element program, like, including "phi_i * x_i" or “grad_phi_i * x_i” for any i of each cell? Thank you very much! Best, Judy On Tuesday, December 7, 2021 at 1:44:38 PM UTC-5 Matthias Maier wrote: > Dear Judy, > > On Tue, Dec 7, 2021, at 11:45 CST, Judy Lee wrote:

Re: [deal.II] Store data from a Tensor<1, dim> to a Vector

o much! I will keep up my work. Best, Judy On Tuesday, December 7, 2021 at 1:44:38 PM UTC-5 Matthias Maier wrote: > Dear Judy, > > On Tue, Dec 7, 2021, at 11:45 CST, Judy Lee wrote: > > > cell_x(i) += fe_values.shape_value(i, q_index) * x(i); > > Here, you are multiplying a sc

Re: [deal.II] Store data from a Tensor<1, dim> to a Vector

to compile, which was made with reference to step-3. The error stops running at Line 132. Thank you so much! Best regards, Judy On Tuesday, December 7, 2021 at 9:54:19 AM UTC-5 Wolfgang Bangerth wrote: > On 12/7/21 03:30, Judy Lee wrote: > > > > But how to use "Point x"

[deal.II] Store data from a Tensor<1, dim> to a Vector

Hello everyone! Can I create "Vector x(n)", to store data values from a Tensor<1, dim> or Point, for physical coordinate after mesh/grid? I can create "Point x", to store data values from a Tensor<1, dim>, it works thru like: Point x; for (const unsigned int i: fe_values.dof_indices())) { x =

Re: [deal.II] Obtain quadrature point value locally

Dr. Bangerth, Yes, your comment is exactly what I need, to resolve my question. I just gave a tiny adjustment on the quadrature parameter, guess what, that returned me exactly my expected data information. This is wonderful, thank you so much! I will keep up my work. Have a great weekend! Judy

Re: [deal.II] Obtain quadrature point value locally

Dr. Bangerth, Thank you so much for the reply that points out a quadrature point may have one or more coordinate. (1) I reset thru step-8 and step-3. Was that sufficient to make the program ready for 1-dimensional(1D) PDE? If not sufficient, I might get more coordinate returned for each

[deal.II] Obtain quadrature point value locally

Hello everyone! How to obtain the value of each quadrature point locally with reference to q_index, (I made this: std::cout << fe_values.quadrature_point <<, in step-3 and step-8, since q_index is looping there), when loops each cell, if my 1D PDE's right-hand-side is a function of x, f(x)?