[deal.II] How to Set Two Constraints on A Single dof Index on Boundaries? (P=0 and \delta P = 0)

2024-07-22 Thread Lex Lee
Hello Deal.ii Friends, I want to set 2 constraints (P = 0 and \delta P = 0) on the scalar pressure variable on a single point (4, 0.1) at the outlet. Here, \delta P is the Newton's iteration update value for P. I used "constraints_hp_nonzero" to interpolate the real boundary values for my phys

Re: [deal.II] Error of Installing Deal.ii 9.5.2 on Intel-based Mac Devices

2024-02-26 Thread Lex Lee
ion is this? What OS version do you have (i.e., is this >>> sonoma?). If you upgraded to sonoma, and are still using the old package, >>> this cannot work, since the compiler has changed with the upgrade. >>> >>> If you use the 9.5.2 package, deal.II comes pre-insta

Re: [deal.II] Error of Installing Deal.ii 9.5.2 on Intel-based Mac Devices

2024-02-05 Thread Lex Lee
> this cannot work, since the compiler has changed with the upgrade. > > If you use the 9.5.2 package, deal.II comes pre-installed. Does that work > for you? > > L. > > > On 3 Feb 2024, at 21:35, Lex Lee wrote: > > > > I got the following error message

Re: [deal.II] Error of Installing Deal.ii 9.5.2 on Intel-based Mac Devices

2024-02-04 Thread Lex Lee
rth wrote: > On 2/3/24 13:35, Lex Lee wrote: > >The C++ compiler > > > > > > > > "/Applications/ > deal.II.app/Contents/Resources/spack/lib/spack/env/clang/clang++" > > > > > >is not able to compile a simple test program. > >

[deal.II] Error of Installing Deal.ii 9.5.2 on Intel-based Mac Devices

2024-02-03 Thread Lex Lee
I got the following error messages -- Detecting CXX compiler ABI info - failed -- Check for working CXX compiler: /Applications/deal.II.app/Contents/Resources/spack/lib/spack/env/clang/clang++ -- Check for working CXX compiler: /Applications/deal.II.app/Contents/Resources/spack/lib/spack/

Re: [deal.II] Errors in get_function_gradients() using PETSc

2023-10-25 Thread Lex Lee
these 3 inputs are not on the locally owned cells. Lex On Wednesday, October 25, 2023 at 12:53:09 PM UTC-7 Wolfgang Bangerth wrote: > > On 10/25/23 13:31, Lex Lee wrote: > > > > Actually it is the output error message from the command line. > > Then how do you know that the

Re: [deal.II] Errors in get_function_gradients() using PETSc

2023-10-25 Thread Lex Lee
I suspect this is some output you get from the debugger? No, not from the debugger. Actually it is the output error message from the command line. Lex On Wednesday, October 25, 2023 at 12:23:32 PM UTC-7 Wolfgang Bangerth wrote: > > > On 10/25/23 10:52, Lex Lee wrote: > > &

Re: [deal.II] Errors in get_function_gradients() using PETSc

2023-10-25 Thread Lex Lee
g Bangerth wrote: > > Lex, > you show the backtrace, but what is the actual error message you observe? > Best > WB > > On 10/24/23 23:08, Lex Lee wrote: > > *** Caution: EXTERNAL Sender *** > > > > Hello deal.ii develpoers, > > > > > > In

Re: [deal.II] Getting cmake errors compiling programs

2023-10-24 Thread Lex Lee
I hope deal.ii group could release version 9.5.1 that works with intel processor on MacOS platform. On Monday, October 9, 2023 at 5:19:03 PM UTC-7 deal.II User Group wrote: > Hello, everyone > > I tried the new image in my laptop with M2 processor and the Ventura Mac > OS (13.6). It worked (I

Re: [deal.II] Difference between 2 Wasy to Initialize System Matrix in MPI World

2023-09-28 Thread Lex Lee
parallel computation. It is inefficient and also likely incorrect if > you run it before distribute_sparsity_pattern(). > Instead, use distribute_local_to_global during assembly (see step-40 for > an example). > > -- > Timo Heister > http://www.math.clemson.edu/~heister/ >

[deal.II] Re: Difference between 2 Wasy to Initialize System Matrix in MPI World

2023-09-28 Thread Lex Lee
Hi Abbas, thanks for the reply. I did debug with option A, but no error reported/no exception is thrown. To me, option A is the same as option B. However, the two generated system matrixs in the code are not the same. On Thursday, September 28, 2023 at 10:27:17 AM UTC-7 abbas.b...@gmail.com wr

[deal.II] Suggestions on PETSc SNES Nonlinear Solver?

2023-09-27 Thread Lex Lee
Hello all, In my FSI problem, there are 4 vector-valued variables + 3 scalar variables. Most of them are strongly coupled with each other. It's hard to decouple the system and design a good preconditoner. Previously, I solved the whole system with the direct solver (UMFPACK) + Newton's itera

[deal.II] Difference between 2 Wasy to Initialize System Matrix in MPI World

2023-09-27 Thread Lex Lee
Hello all, Previously, I initialized my system matrix in this way (option A) DoFTools::make_flux_sparsity_pattern(dof_handler, dsp, cell_coupling, face_coupling); constraints_newton_updat

Re: [deal.II] Is It Possible to Make Step-46 Run in Parallel?

2023-09-14 Thread Lex Lee
"You record constraints in an AffineConstraints object, and that continues to be the work in parallel as well" Thanks a lot, Prof. Bangerth. I may misunderstand the error messages. I will try to debug on other parts and keep this topic updated. Best, Lex On Wednesday, September 13, 2023

[deal.II] Is It Possible to Make Step-46 Run in Parallel?

2023-09-11 Thread Lex Lee
Hello all, I developed my FSI code based on the FENothing framework used in step-46. Step-46 shows how to identify the interface between the fluid and solid domain. Also it shows how to set constraints on variables and how to assemble terms on the interface. I just used the same logic in my FS

[deal.II] Request of Step-77-snes.cc File

2023-08-06 Thread Lex Lee
obian-free Newton Krylov method, and want to play with PETSc(snes) library. Thanks a lot. Best, Lex Lee -- 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

[deal.II] Connect UMFPACK with Trilinos in deal.II

2023-02-08 Thread Lex Lee
Hello all, I am using UMFPACK to directly solve my large system equation Ax=b. The solver process is very slow, as only one core of CPU is involved while running. I am trying to connect UMFPACK with Trilinos, I mean, solve the equation in multithread / parallel mode in deal.II. I know some tut

Re: [deal.II] How to transform points on faces from unit cells to real cells?

2020-06-22 Thread Lex Lee
the support points in real space: >> >> quadrature_point() >> get_quadrature_points() >> >> Best, >> Simon >> >> On Monday, June 22, 2020 at 10:53:10 PM UTC+2, Lex Lee wrote: >>> >>> Hello Doug, >>> >>> >>> Than

Re: [deal.II] How to transform points on faces from unit cells to real cells?

2020-06-22 Thread Lex Lee
> Best, > Simon > > On Monday, June 22, 2020 at 10:53:10 PM UTC+2, Lex Lee wrote: >> >> Hello Doug, >> >> >> Thanks for your help. >> >> However, I need to say, I just had done exactly what you described before >> I posted this question.

Re: [deal.II] How to transform points on faces from unit cells to real cells?

2020-06-22 Thread Lex Lee
t; > Doug > > On Sunday, June 21, 2020 at 8:10:30 PM UTC-4, Lex Lee wrote: > Hello Doug, > > Thanks for your kind help. > > I am trying to understand your suggestion in the original email. > > Are you suggesting me using the member function > "transform

[deal.II] Re: How to transform points on faces from unit cells to real cells?

2020-06-21 Thread Lex Lee
o project a "volume" point, but projecting within > unit cell is easy, which you would then follow with > transform_unit_to_real_cell(). > > Best regards, > > Doug > > On Saturday, June 20, 2020 at 8:16:22 PM UTC-4, Lex Lee wrote: >> >> Hello Deal.II

[deal.II] How to transform points on faces from unit cells to real cells?

2020-06-20 Thread Lex Lee
Hello Deal.II Users, I want to get the physical (geometry) coordinates for support points on cell faces. Also I know that there are such member functions that can me map points between reference and real cells in this link: https://www.dealii.org/current/doxygen/deal.II/classMapping.html . H

[deal.II] Re: How to get normal / tangential vectors at nodes, not at quadrature points?

2020-06-05 Thread Lex Lee
t; > You can then use this quadrature in FEFaceValues to get the normals. Here, > the weights of the dummy quadrature are inf, but if the normals are > everything you need that isn't a problem. > > Best, > Simon > > On Thursday, June 4, 2020 at 8:09:08 PM UTC+2, Lex L

[deal.II] How to get normal / tangential vectors at nodes, not at quadrature points?

2020-06-04 Thread Lex Lee
Hello Deal.II Users, I am working on setting a constraint: (V - v_s ) \cdot n = \phi_f (v_f - v_s) \cdot n on three vector variables at the interface with affine constraint class in Deal.II's library. (Where, V, v_s, and v_f are velocity vectors on two different domains, and they are coupl

[deal.II] How to simplify a 2D complex problem to a '1D-like' problem?

2020-04-24 Thread Lex Lee
Hello Deal.II Users, I am using Deal.II to try to simplify/reduce my complex 2D problem to a 1D-like problem, namely, let the component values at x-direction of all vectors be zero. (for example, V=[v_x, v_y]' = [0, v_y]') . Moreover, given that the governing equations, boundary conditions a

[deal.II] Debug to identify pairs of DoFs that are co-lated and correspond to the components.

2020-04-05 Thread Lex Lee
Hello all, I designed a coupled Laplace problem to play with ConstraintMatrix for coupled components. On face \Gamma_4 and \Gamma_2, I let u=0.5v-0.5, v=2u+1, respectively . [image: Screen Shot 2020-04-05 at 4.14.41 PM.png] Both u and v use the same finite element: FE_Q<*dim*>(1),

[deal.II] How to set up constraints for several different variables at the boundary?

2020-03-27 Thread Lex Lee
Hello all, In order to handle the constraint : (v-v_s) \times n=\phi (v_f -v_s) \times n at the boundary in my research problem, I designed the following test problem to play with ConstraintsMatrix in deal.ii. (n is the normal vector, and v, v_s, v_f are the velocities of three phases. They

Re: [deal.II] How to Make Test Functions Satisfy Certain Constraints

2020-02-17 Thread Lex Lee
rent/doxygen/deal.II/Tutorial.html. > > Best, > Daniel > > Am So., 16. Feb. 2020 um 23:45 Uhr schrieb Lex Lee >: > >> Hello all, >> >> >> To numerically solve the governing equations I am studying now, I notice >> that the chosen three test functions s

[deal.II] How to Make Test Functions Satisfy Certain Constraints

2020-02-16 Thread Lex Lee
. If so, can anyone give me some hints, I mean, study/reading material? Best, Lex Lee -- 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