[deal.II] Assemble Righthand Side for vector-valued problem

tor, *^* */Users/jaekwangjk/Programs/dealii/include/deal.II/base/symmetric_tensor.h:3086:1: note: * candidate template ignored: could not match 'SymmetricTensor<2, dim, type-parameter-0-1>' against 'const double' operator * (const SymmetricTensor<2,dim,Number> , *^* 1

Re: [deal.II] Assemble Righthand Side for vector-valued problem

017 12:35 PM, Jaekwang Kim wrote: > >*/_local_rhs(i) += fe_values.shape_value(i,q) *_/* > > > > */_rhs_values[q] *_/* > > > > */_fe_values.JxW(q);_/* > > > > rhs_values[q] is

[deal.II] Assemble system with variational right hand side

Hi. I have a question on how we can impose external force term, or right hand side, in Stokes flow. My governing equation is same with that of step-22, stokes flow -nabla^2 u + grad p = f -grad u = 0 but I want to impose 'f', To do this, I made function that gives me f value on each

[deal.II] Re: Assemble system with variational right hand side

Thanks for the reply, I just wanted to make sure that whether I implemented it right, but I still cannot sure whether I am doing correctly. Let me be with more specific questions. First, Yes I am considering 2-dimensional flow. if I writes a code that looks as previous, template

Re: [deal.II] Assemble Righthand Side for vector-valued problem

Thank you for your advice. I would make the problem simpler. Do you get the right solution if you > have a constant viscosity? If you iterate, do you get the solution after > one iteration? > Yes, I checked this. After one iteration, I get the solution after one iteration when constant

Re: [deal.II] Assemble Righthand Side for vector-valued problem

Dr. Bangerth Thank you, I just fixed up what was wrong... in most cases, it is usually easy problems. The manufactured solution of velocities has not satisfies the continuity equation. (i.e. manufactured solution does not satisfy, div.u=0)... Jaekwang Kim -- The deal.II project

[deal.II] Relation between Solution Error Behavior and Polynomial Approximation Degree

o node point for this) is there any suggestion that I need to fix more? or my first prediction that the slope will show '-p' or error will just behave O(h^p) was wrong? I am always thank you for all guys! Jaekwang Kim -- The deal.II project is located at http://www.dealii.org/ For mailin

Re: [deal.II] VectorTools::integrate_difference function

Thank you for all the repliers! All of you gave me a good advice that I may consider deeply when I write my code from now I have figured out what was wrong. Seeing at Step 20 code, I realized that I have to select component of solution before going to compare the difference. const

[deal.II] Re: Can we say that the higher order method, the more accurate?

Thank you for all repliers!! I tried fix code as you mentioned please see redline... template void StokesProblem::assemble_system () { system_matrix=0; system_rhs=0; *QGaussLobatto<**2> quadrature_formula(degree+2);* FEValues fe_values (fe, quadrature_formula,

[deal.II] Re: Can we say that the higher order method, the more accurate?

2016년 9월 18일 일요일 오전 10시 57분 16초 UTC-5, JAEKWANG KIM 님의 말: > > > Hello, I am a starter of dealii and am learning a lot these days with the > help of video lectures and tutorial examples. > > I modified step-22 code (stokes flow code) into my own problem, the flow > aro

Re: [deal.II] VectorTools::integrate_difference function

yes...I tried that before, but I get an error message as follow.. As you mentioned, I first learn this function in step 7 tutorial. But the problem is that I want to integrate my vector-valued-solution to exact-vector-valued solution, while step 7 compares scalar-solution. Or is there any

[deal.II] Re: Can we say that the higher order method, the more accurate?

I also wonder whether their is a limitation in higher order method that deal offer. If no, then we can use Q10 method, just putting degree value as 10? -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see

[deal.II] Re: Can we say that the higher order method, the more accurate?

This is how my mesh looks like. I am using half domain of 2D sphere assuming that it would symmetrical along phi domain. Tho my domain is 2d, I am

Re: [deal.II] VectorTools::integrate_difference function

egree+2) or so, which increases > the order of quadrature as the element order gets higher. > > Best, > Martin > > On 09/20/2016 04:46 PM, JAEKWANG KIM wrote: > > > > I wanted to evaluate Vectortool::integrate_difference to estimate my > erro

Re: [deal.II] Re: Can we say that the higher order method, the more accurate?

Can I ask more about "the limits of our implementation of our quadrature formulas?". I wonder when it usually happens. Once I calculated drag coefficient, with Q1, my error is 10% compare to exact solution. However, I can 1.7% error when I use Q2 which is significant decrease! At Q3, I get

[deal.II] Gauss-Lobatto node distribution on FESYSTEM?

I want to use Gauss-Lobatto node on my fesystem for vector valued problems. While seeing step-48, I see that how Gauss-Lobatto node is used but I couldn't figure out how I should use it for vector value problems. To be specific, In step 48, fe is declared in main class as follow

[deal.II] Re: How to use sacado package

to > configure Trilinos without fortran using > -DTrilinos_ENABLE_Fortran:BOOL=OFF > > After that it's the usual make install [1]. > > Best, > Daniel > > [1] https://www.dealii.org/developer/external-libs/trilinos.html > > Am Samstag, 22. Oktober 2016 01:44:26 UT

Re: [deal.II] Re: How to use sacado package

Thanks for the comment. I'd like to follow the first option to install MPI on my system. It seems that there are different options that I can choose. Ready page in deal.ii linked me to the site "http://mpi-forum.org; and it seems that there are different packages of mpi. Does it matter to

Re: [deal.II] Re: How to use sacado package

Always thank you for the comments!! I'd like to follow the first choice, installing MPI on my machine. I installed openmpi-2.0.1. Coming to the original problem, I tried install trillions 12.4.2 again. but now I received this error message while i am installing trilinos. what might be the

[deal.II] Re: Easy way to calculate second-deviatoric tensor

Thank you all guys who replied here!! -- 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

[deal.II] How to use sacado package

I'd like to use sacado package for automatic differentiation. However, I have installed deal.ii on my computer with trilinos configuration set off. I looked back deal.ii installation manual and got to know that I have to install trilnos package first. However, when I try to install

[deal.II] Some questions on Nonlinear Problems

error However, if I use enough small tolerance between abs(u_{k}-u_{k+1}), will the error will also show ~O(h^p)? Always thank you all and I am learning a lot form this amazing community!! Jaekwang Kim -- The deal.II project is located at http://www.dealii.org/ For mailing list

[deal.II] Very basic question on setting boundary id

2016-11-13%2Bat%2B12.59.38%2BPM.png> Always thank you all you for fast response! Jaekwang Kim -- 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 subscr

[deal.II] Re: Very basic question on setting boundary id

Thank you for quick reply!!! I tried your suggestions but I couldn't resolve the problem. The more weird thing is if I commands as.. GridGenerator::hyper_cube (triangulation, 0, 1); triangulation.refine_global (); for (typename Triangulation::active_cell_iterator

[deal.II] Re: Very basic question on setting boundary id

Thank you for being kindness! I followed your suggestions but I couldn't resolve the problem yet. Or did I make other mistake in another places? >From a slight modification from STEP-4 tutorial , I have been trying solve Poisson Equation. Actually, I have manufactured solution to see my

[deal.II] Re: Very basic question on setting boundary id

Problem Solved! Thank you!! 2016년 11월 13일 일요일 오후 5시 25분 20초 UTC-6, Daniel Arndt 님의 말: > > Jaekwang, > > It seems that you are missing > constraints.distribute(solution); > after solving the linear system. The constrained dofs are set to zero > during assembly and > this call is necessary to

Re: [deal.II] Re: hp convergence for non-linear solver.

Do you have any idea about this? Always thank you for your time and advice! Best regards, Jaekwang Kim 2016년 11월 2일 수요일 오후 8시 30분 19초 UTC-5, Wolfgang Bangerth 님의 말: > > On 11/02/2016 07:17 PM, Jaekwang Kim wrote: > > > > I'd like check the visualized error per cell. &

[deal.II] hp convergence for non-linear solver.

couldn't get such convergence for the case p=2, and 3. I have tried to change tolerance number in iterative scheme, but it didn't work. Are there other things that I failed to consider? Or did I make wrong way to test the code? Always thank you for all smart guys here!! Jaekwang Kim

[deal.II] Re: hp convergence for non-linear solver.

t; > On Tuesday, November 1, 2016 at 5:43:36 PM UTC-4, Jaekwang Kim wrote: >> >> I think that I didn't make any mistake on my code lines because at least >> my numerical solution is converging for the case p=1. >> >> but I couldn't get such convergence for the case p=2, an

[deal.II] Re: hp convergence for non-linear solver.

(), difference_per_cell,QGauss(degree+2),VectorTools::L2_norm); How can I use non-default Mapping in here? (I wonder specific lines of code...) Thank for giving me a help! Regards, Jaekwang Kim 2016년 11월 1일 화요일 오후 5시 19분 35초 UTC-5, Daniel Arndt 님의 말

[deal.II] Easy way to calculate second-deviatoric tensor

SAa2IYsdK4NCMOuY3QXJydmwCLcB/s1600/Screen%2BShot%2B2016-10-16%2Bat%2B6.42.40%2BPM.png> Regards, Jaekwang Kim -- 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 mes

[deal.II] Re: Postprocessing on Vector Valued Problem

You were right Oh I just got the result! I didn't know that dof_handler can automatically matching the sequence of the cell. Thank you for your time Jaekwang Kim -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d

[deal.II] Vector Value Problem Solution initialization

Hi all. I have a question on the solution initialization for iterative methods, in vector valued problem. For example, I want to solve 2D flow, where my solution consist of two block, one is for velocities and the other is for pressure. First, I construct my initial solution as below

Re: [deal.II] Re: AMR , how to pass solution vector to refined mesh

Thank you!! Just got what've meant. and fixed my code properly Thanks ! Jaekwang Kim 2017년 3월 21일 화요일 오후 12시 55분 7초 UTC-4, Jean-Paul Pelteret 님의 말: > > Hi Jaekwang, > > No, you've misunderstood in that I was just showing that to you as an > example. You don't need to use th

[deal.II] is there a way to refine mesh only in one direction

nks. Bests, Jaekwang Kim <https://lh3.googleusercontent.com/-QyyWHbmWMC4/WNwDRpHKeEI/BAQ/3NlDrVNVqM0KsSdDDdu4XgePeI846Zj1gCLcB/s1600/Picture1.png> -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/de

[deal.II] Re: Timo Heister and I on a podcast, talking about deal.II

Nice interview. I really enjoyed the cast. -- 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.

Re: [deal.II] Re: is there a way to refine mesh only in one direction

Bruno Turcksin 님의 말: > > Jaekwang, > > 2017-03-30 17:39 GMT-04:00 Jaekwang Kim <jaekw...@gmail.com >: > > > //This one is problematic > > triangulation.prepare_coarsening_and_refinement (); > It could be a bug in the library, most

Re: [deal.II] Re: AMR , how to pass solution vector to refined mesh

problem solution (u,v,p) each by each to refined mesh. How can I do this...? (Is there a short way that I can assess each solution vector that live in BlockVector ? Always thank you. Jaekwang Kim -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options

[deal.II] Re: AMR , how to pass solution vector to refined mesh

Thanks for the response! I tried to put "Assert (system_matrix.n() == dof_handler.n_dofs(), ExcInternalError()); " when I assemble system , but i get the same error message, regardless of it *An error occurred in line <1668> of file in function* *void

Re: [deal.II] Re: AMR , how to pass solution vector to refined mesh

,* */*keep_constrained_dofs = */ false);* 2017년 3월 20일 월요일 오후 1시 3분 50초 UTC-4, Wolfgang Bangerth 님의 말: > > On 03/20/2017 10:58 AM, Jaekwang Kim wrote: > > > > I tried to assert assemble function as you suggested, but I didn't > > receive any error

Re: [deal.II] Re: AMR , how to pass solution vector to refined mesh

Thanks for your response. I tried to assert assemble function as you suggested, but I didn't receive any error message. What make me more confuse me is the code is working with global refinement. if I input 100% cells to be refined.. GridRefinement::refine_and_coarsen_fixed_number

[deal.II] Re: AMR , how to pass solution vector to refined mesh

"Set up system finished" << std::endl; } } I think this two function do everything I needed... but when I run my code, I run into error message when I assemble system in second time. what might be the problem.? Always thank you!! Jaekwang Kim *An error occurred in line <166

[deal.II] AMR , how to pass solution vector to refined mesh

is refined is there anyone who has idea on this ? Thanks, Jaekwang Kim template void nonlinear::run () { unsigned int cycle_control=4; for (unsigned int cycle=0; cycle<cycle_control; ++cycle) { std::cout << "Cycle " <<

[deal.II] any suggestion how to run same code repeatedly with different boundary condition ?

Hi all, I want to run same problem for different boundary condition values that will come globally... for example, if i want to run step-4 code with different boundary condition, #include #include #include #include #include #include #include #include #include #include

[deal.II] Re: any suggestion how to run same code repeatedly with different boundary condition ?

Thank you !! solve the problem with the second method! I really appreciate this! Jaekwang Kim 2017년 3월 4일 토요일 오후 1시 37분 40초 UTC-6, Jean-Paul Pelteret 님의 말: > > Hi Jaekwang, > > The error message is pointing you to the fact that when you call run() a > second time you end up tr

Re: [deal.II] Re: is there a way to refine mesh only in one direction

Bangerth 님의 말: > > On 03/31/2017 01:30 PM, Jaekwang Kim wrote: > > > > I found that what was the problem. > > > > I was using 'maximum' Meshsmoothing which does not allow > > > > So, It was not a bug of deal.ii but I was using wrong smoothing method

Re: [deal.II] Re: AMR , how to pass solution vector to refined mesh

efinement(x_temperature); > stokes_trans.prepare_for_coarsening_and_refinement(x_stokes); > > I hope this helps, > Jean-Paul > > On Tuesday, March 21, 2017 at 12:26:59 AM UTC+1, Jaekwang Kim wrote: >> >> Dear, Dr. Bangerth >> >> Adding up to this, may I ask yo

Re: [deal.II] Distribute constraints when using Workstream class

I just understand what you meant ! Thanks ! 2017년 4월 6일 목요일 오후 7시 24분 1초 UTC-5, Wolfgang Bangerth 님의 말: > > On 04/06/2017 02:55 PM, Jaekwang Kim wrote: > > template > > > > void > > > > Step4::copy_local_to_global (constAssemblyCopyData _data) > &g

[deal.II] Re: advection equation weak form in tutorial 9

Thanks, now I clearly understand what's going on! Jaekwang On Monday, August 14, 2017 at 3:08:18 PM UTC-5, Jaekwang Kim wrote: > > Hi, all > > I am little bit confused by weak form for advection problem written in > step-9 tutorial > > in the its tutorial, in the middle

[deal.II] advection equation weak form in tutorial 9

Hi, all I am little bit confused by weak form for advection problem written in step-9 tutorial in the its tutorial, in the middle of the line, it says that... I am little bit unsure the red-marked form, >From the whole boundary integral, we put inflow case, to the right hand side, and

Re: [deal.II] how to use post_processor for FE_System

value I am receiving error as... *grad_u: 0 , 0* *make[3]: *** [CMakeFiles/run] Segmentation fault: 11* *make[2]: *** [CMakeFiles/run.dir/all] Error 2* *make[1]: *** [CMakeFiles/run.dir/rule] Error 2* *make: *** [run] Error 2* *Thanks...* *Jaekwang * On Wednesday, July 19, 2017 at 6:02:14 PM UTC-5,

[deal.II] Datapostprocessor first derivative based on what?

meet. (My taylor element is degree 2 for velocity and degree 1 for pressure. ) If then, on those vertexes, how the Datapostprocessor decides which element to refer for the first, or second derivatives on that dof? Thank you Jaekwang Kim -- The deal.II project is located at http

[deal.II] how to use post_processor for FE_System

Hi, all I was practicing using post_processor, referring step-29 and step 31. If I want to use post_processor for vector-valued problem, where my solution is (u,v,p) and I want to compute magnitude of velocity, |u^2+v^2| ignoring pressure parts. i wrote as step 29; template class

Re: [deal.II] how to use post_processor for FE_System

MagnitudeILi2EE is still used by 1 other objects.* *(Additional information: )* *See the entry in the Frequently Asked Questions of deal.II (linked to from http://www.dealii.org/) for a lot more information on what this error means and how to fix programs in which it happens.* On Wednesd

Re: [deal.II] Re: Can we say that the higher order method, the more accurate?

e >> discretization order and refine the mesh. >> >> >> <https://lh3.googleusercontent.com/-Gz5932Zt6e0/WYqvbO3P-4I/AAAABSM/6hBHn1FO5W0P7X3SPfQ4iyRaldDYzOt3QCLcBGAs/s1600/Image.png> >> >> As is shown in the figure, the Cd increased as the cells increased for

[deal.II] Ideas on accuracy of Streamline diffusion method

order for integrating the weak form.. Any one has idea on possible point I will have to check to attain (p+1/2) error behavior ? Always thank you, Regards, Jaekwang Kim -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https

[deal.II] Weak form for advection - diffusion equation

Since since my question include some mathematical expression, I attached figure for this...

[deal.II] Re: weak form for advection-diffusion equation

Adding up to this, I was trying as... scratch_data.fe_values.shape_grad_grad (j,q_point) but keep receiving errors */Users/kimjaekwang/repos/advection_diffusion/step9/step-9.cc:487:53: **error: **no member named 'shape_grad_grad' in* * 'dealii::FEValues<2, 2>'* *

[deal.II] Re: weak form for advection-diffusion equation

tion term. > > Best, > > Bruno > > On Wednesday, June 14, 2017 at 3:33:18 PM UTC-4, Jaekwang Kim wrote: >> >> >> <https://lh3.googleusercontent.com/-qsDNe2Fge5w/WUGPdbVCfRI/BCI/egHPb83UKLIjnTGxmMfHiI4ZKGgXq4FhgCLcBGAs/s1600/Screen%2BShot%2B2017-06-14%2Bat%2B2

[deal.II] Distribute constraints on PETSc wrapper

Hi, all I was practicing using Petsc library through deal.ii wrappers, While, step 17 gives me clear explanation, I am having trouble in distributing hanging node constraints on Petsc system matrix, I was modifying step-9 hyperbolic equation solvers with PETSc matrix, but when I compile my

[deal.II] Step-9, iterator over the cell

Hi, all. I have a questions on an algorithm in step-9 that does some iteration over the cell with its neighbors. I tried to read explanation carefully, but

[deal.II] Re: Relation between Solution Error Behavior and Polynomial Approximation Degree

.ii to generate the figure. If not, what did you use? > > Thanks, > Jiaqi > > 在 2016年9月29日星期四 UTC-4上午11:41:48，Jaekwang Kim写道： >> >> >> <https://lh3.googleusercontent.com/-dkfwJvgPEmc/V-00Y-c4tzI/A5s/pb0Z-IBRe-An_S39dSDzxWg1k5GJlEbKQCLcB/s1600/step7plotting.jpg>

[deal.II] use of "meshworker" for FESystem

Hi, all. I am curious about whether i can use MeshWorker: class for FESystem. I am trying to solve system, that I need to solve 'stokes flow' and then use velocity solution to solve other scalar transport equation using DG method that looks like

[deal.II] Re: Use of 'VectorTools::interpolate_boundary_values' using numerical solution.

me this function needs to know the solution, > triangulation etc. Then you would have to create some references or dealii > SmartPointers in this class as well. > > I hope this helps. > > Jie > > On Monday, April 9, 2018 at 6:32:28 PM UTC-4, Jaekwang Kim wrote: >>

[deal.II] Use of 'VectorTools::interpolate_boundary_values' using numerical solution.

Hi, all I am always thank you for all guys here in this community. I have a question on the use of VectorTools::interpolate_boundary_values functions Typically we use this function in when we set-up system as in the form VectorTools::interpolate_boundary_values (dof_handler,

[deal.II] Adding rows on system matrix connecting to dof

Hi, all I have a question on using deal.ii for solving unknown boundary condition problem. I have thought the most minimized related to my question. where I consider a navies slip boundary condition, and I have tried to solve this with finite difference. For more demonstration please see

Re: [deal.II] interpolating boundary values from a previously found FE solution

Dr. Bangerth Is there any example tutorial related to the first approach? *"You can build a ConstraintMatrix object that sets boundary DoFs of V to the corresponding values of U. This should be easy to do." * On Friday, March 30, 2018 at 8:15:23 AM UTC-5, Wolfgang Bangerth wrote: > > On

Re: [deal.II] Re: use of "meshworker" for FESystem

Thanks for the attetion. I overcame this problem, though I don’t recognize what was solution now. It turns out to be possible… Thank you for your time Jaekwang On Sun, Mar 25, 2018 at 3:13 AM, Daniel Arndt < daniel.ar...@iwr.uni-heidelberg.de> wrote: > Jaekwang, > > do you have a minimal

Re: [deal.II] Discrete Galekin with periodic boundary condition

itbucket.org/cpraveen/deal_ii/src/master/dg/1d_scalar_legendre/ > > Best > praveen > > On 08-Nov-2018, at 8:31 AM, Jaekwang Kim > > wrote: > > > > I am trying to solve advection equation using DG method with Meshworker. > > Yet, I have a problem in implem

[deal.II] Re: Discrete Galekin with periodic boundary condition

Thanks, for the reply.. Hi am using following function to construct periodicity when I set up dof DoFTools::make_hanging_node_constraints (dof_handler, constraints); DoFTools::make_periodicity_constraints(dof_handler,6,5 //Boundary ID 6 is outlet and 5 is inlet

[deal.II] Re: Discrete Galekin with periodic boundary condition

It seems that I don't need to add_periodicity to the triangulation unless it is parallel distributed triangulation? On Friday, November 9, 2018 at 9:45:53 AM UTC-6, Jaekwang Kim wrote: > > Thanks! Now it becomes clear to me. > > By the way, when I 'add_periodicity', I have co

[deal.II] Re: Discrete Galekin with periodic boundary condition

You were right. I couldn't get the correct results without Triangulation::add_periodicity(). So, it seems that I should compile deal.ii with MPI,P4EST,Trilinos to use that function(Triangulation::add_periodicity()) since my triangulation should be parallel distributed... Right? Best,

[deal.II] Re: Discrete Galekin with periodic boundary condition

No, you were right, I couldn't get the correct results without 'Triangulation::add_periodicity()'. The current my local environment does not configure deal.ii with MPI and Trilinos as the example step-45, which shows how this function is used. But I guess I can still use the "

[deal.II] Re: Discrete Galekin with periodic boundary condition

Thanks! Now it becomes clear to me. By the way, when I 'add_periodicity', I have compiler error on the red line of the code. Can you give me a clue to resolve this error? //original argument std::vector::cell_iterator>> pairs;

[deal.II] Stokes flow with Periodic boundary condition

Hi, all. I'd like ask a question implementing periodic boundary condition. Modifying step-22 tutorial code, I want to solve a stokes flow at Couette configuration I plan to solve it at (x,y) in [0,5] x [0,1], rectangular domain with boundary condition \vec{u} (x,1)= (1,0) (moving

Re: [deal.II] Re: Install error when compile with Trilinos

Thanks I complied it successfully On Sun, Feb 24, 2019 at 10:55 PM Bruno Turcksin wrote: > Jaekwang, > > This is the problem: > > Le dim. 24 févr. 2019 à 22:10, Jaekwang Kim a > écrit : > > > -- Include > /Users/jaekwangkim/Program/dealii-9.0.1/cmake/config

[deal.II] Re: Install error when compile with Trilinos

I resummarized screen messages that seems to be a problem ...~~ *-- Performing Test DEAL_II_HAVE_AVX* *-- Performing Test DEAL_II_HAVE_AVX - Failed* *-- Performing Test DEAL_II_HAVE_AVX512* *-- Performing Test DEAL_II_HAVE_AVX512 - Failed*

[deal.II] Install error when compile with Trilinos

Hi, all. I have used deal.ii for year, and recently I need to compiled deal.ii with Trilinos I installed Trilinos on my local enabling only serial (non MPI) However, after the installation, when I tried to reinstall with deal.ii enabling Trilinos packages, I got error so that I cannot

[deal.II] ILU preconditioner and mesh-scale dependency?

solution.block(2), system_rhs.block(2), preconditioner); // constraints.distribute (solution); } I wonder now ... 1. why smaller scale mesh only results such error, 2. If possible, is there any choice of better preconditionner for my system other than SparseILU to go around th

Re: [deal.II] ILU preconditioner and mesh-scale dependency?

Thanks for sharing insight! On Thursday, March 14, 2019 at 12:18:30 PM UTC-5, Wolfgang Bangerth wrote: > > On 3/13/19 8:18 PM, Jaekwang Kim wrote: > > > > > > I wonder now ... > > > > 1. why smaller scale mesh only results such error, > > Y