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

is better than Q2. But Q4 or Q4 is not better than Q2 or Q3? Is there any reason on this? (To be specific, if i say degree 2 , that mean I use (2+1) for velocity, (2) for pressure, and (2+2) for Gauss integral Thank you Jaekwang Kim -- The deal.II project is located at http://www.

[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 describing

[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 https://groups.google.com/d/forum/dealii?hl

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 0

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

thank for the reply!! my fe degree is declared as "fe (FE_Q(degree+1), dim,FE_Q(degree), 1)" and I used "QGauss quadrature_formula(degree+2);" to calculate integral over the cell. 2016년 9월 20일 화요일 오전 9시 26분 6초 UTC-5, Praveen C 님의 말: > > > On Tue, Sep 20,

[deal.II] VectorTools::integrate_difference function

I wanted to evaluate Vectortool::integrate_difference to estimate my error. but I couldn't figure out what I need to enter for the third component. (I marked it as red) VectorTools::integrate_difference (dof_handler, solution,

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

l and make sure to use QGauss(fe.degree+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

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 way

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

tc.) Are the solver tolerances tight enough? > What happens if the mesh is refined? If this still does not help, it might > be good to share a small example where the issue is observed. > > Best, > Martin > > > On 09/20/2016 04:35 PM, JAEKWANG KIM wrote: > > thank for the

[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

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 ComponentSel

[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 templat

[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/

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

Yes! your explanation clear explains the situation ! Thank you!!! Jaekwang Kim 2016년 9월 29일 목요일 오전 11시 25분 28초 UTC-5, Wolfgang Bangerth 님의 말: > > On 09/29/2016 09:41 AM, JAEKWANG KIM wrote: > > > > I thought that the order of error is O(h^p) where h is a mesh-size and p

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

JydmwCLcB/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 message because you are subscribe

[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 uns

[deal.II] Some questions on Nonlinear Problems

ght also affect 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/ F

[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 tr

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

d be able 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

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 p

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

. but I 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

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

ekwang, > > 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

[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] Re: hp convergence for non-linear solver.

ave any other idea? Thanks again for your time Regards, Jaekwang Kim 2016년 11월 1일 화요일 오후 5시 19분 35초 UTC-5, Daniel Arndt 님의 말: > > Jaekwang, > > [...] >> >> VectorTools::integrate_difference >> (dof_handler,solution,Solution(), >> >>

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

Thank you for replying me! I'd like check the visualized error per cell. Is there any module in deal.ii that enables this...? 2016년 11월 2일 수요일 오전 7시 27분 56초 UTC-5, Wolfgang Bangerth 님의 말: > > On 11/01/2016 11:06 PM, Jaekwang Kim wrote: > > > > do you have any other

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 pe

[deal.II] Question on Vectortools::integrate_difference.

input of which guadrature to use to integrate this ...I don't understand where do we need this quadrature... what is the role of 6th argument of this function ? QGauss(degree+2) Thank you always ! Jaekwang Kim -- The deal.II project is located at http://www.dealii.org/ For mailin

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

wrong? I thought that the faces of cell is each edge of a cell, and examined whether their centers are x=0,1 or y=0,1. and If then, I set id on 1... <https://lh3.googleusercontent.com/-toZV8W-J-lw/WCi4Kldx3KI/A88/OF-tsoWgJpILwNta9WdYtUEyxMx77EBsQCLcB/s1600/Screen%2BShot%2B2016

[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 get

[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 *-

[deal.II] Postprocessing on Vector Valued Problem

k+=1; } } I could able to - shear_rate at each cell appropriately, but How can I relate my 'cellwise_shear_rate' with dof handler so that I finally able to see the values in 'vtk' format? for now, it just seems that no more than a list of arrays

[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.

[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 poin

[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 voi

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

you are right. Sorry about that. I may test and learn more by myself. 2017년 1월 23일 월요일 오후 1시 21분 56초 UTC-6, Bruno Turcksin 님의 말: > > Hi, > > On Monday, January 23, 2017 at 2:06:17 PM UTC-5, Jaekwang Kim wrote: >> >> I just wanted to make sure that whether I impleme

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

(const SymmetricTensor &t, *^* */Users/jaekwangjk/Programs/dealii/include/deal.II/base/symmetric_tensor.h:2791:1: note: * candidate template ignored: could not match 'SymmetricTensor' against 'Vector' operator * (const double factor, *^

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

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 visc

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 is

[deal.II] can we bring variables from outside class...?

template double BoundaryValues::value (const Point &p, const unsigned int /*component*/) const { double u = u_main; *<-obviously, this dose not **work,showing error message that use of undeclared identifier 'u_main'* retur

[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 #inc

[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

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

sh 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::active_cell_iterator cell=triangulation.begin_active();

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

} } 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 <1668> of file in function* *void dea

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

that I am reinitialize matrixes... but I am not sure about whether the redline is doing enough for reinitialization of sparsity pattern Jaekwang Kim 2017년 3월 17일 금요일 오전 8시 30분 21초 UTC-5, Wolfgang Bangerth 님의 말: > > On 03/16/2017 08:38 PM, Jaekwang Kim wrote: > > > >

[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 dealii::internals::dealiiSparseM

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 (tria

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

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

to transfer my fluid 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 mailin

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

); 2. solution_transfer.prepare_for_coarsening_and_refinement(solution.block(0 )); will there be other easy way I can transfer vector-valued solution without having this difficulty? I always appreciate your help ! Thanks Jaekwang Kim *An error occurred in line <253> of file in function*

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

://www.dealii.org/8.4.1/doxygen/deal.II/classTriangulation.html#ab9fa3177e0e43ab0cf243215d284a35a> > (); > temperature_trans.prepare_for_coarsening_and_refinement(x_temperature); > stokes_trans.prepare_for_coarsening_and_refinement(x_stokes); > > I hope this helps, > Jean-Paul > > On Tuesday,

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

[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: is there a way to refine mesh only in one direction

loper/doxygen/deal.II/classCellAccessor.html#afb6cc537720a5b6381c237abe0887de2) > > but only in serial. So depending if your code is parallel or not this may > be what you are looking for. > > Best, > > Bruno > > On Wednesday, March 29, 2017 at 3:01:49 PM UTC-4, Jaekwang Kim wrote: >> &

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

l.II, not in your own program. Try to come up with the smallest possible program that still demonstrates the error and contact the deal.II mailing lists with it to obtain help.* 2017년 3월 29일 수요일 오후 10시 20분 24초 UTC-5, Wolfgang Bangerth 님의 말: > > On 03/29/2017 06:37 PM, Jaekwang Kim wrote:

[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

37초 UTC-5, Bruno Turcksin 님의 말: > > Jaekwang, > > 2017-03-30 17:39 GMT-04:00 Jaekwang Kim >: > > > //This one is problematic > > triangulation.prepare_coarsening_and_refinement (); > It could be a bug in the library, most people don't use

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

fgang 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 smo

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] dealii, installation error with mpi

Hi all, I was trying to install deal.ii with MPI on cluster, but I met error. *-- Include /home/jk12/Programs/dealii-8.5.0/cmake/configure/configure_hdf5.cmake* *-- Found HDF5_INCLUDE_DIR* *-- Found HDF5_LIBRARY* *-- Found HDF5_HL_LIBRARY* *-- Found HDF5_PUBCONF* *-- HDF5_LIBRARIES: /u

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

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

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

-- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/

[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

on 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%2Ba

[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 Mag

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

) for a lot more information on what this error means and how to fix programs in which it happens.* On Wednesday, July 19, 2017 at 5:23:36 PM UTC-5, Wolfgang Bangerth wrote: > > On 07/19/2017 04:09 PM, Jaekwang Kim wrote: > > > > because my dof_handler has doc for pressure, s

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

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, Wolfgang Bangerth wrote: > > On 07/19/2017 05:00 PM, Jaekwang Kim wrote: >

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

t elements 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 locat

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

rease the >> 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

[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 rem

[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 m

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

used enough Gauss Quadrature 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

[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 some

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

oh Thanks, Now I understand what it had meant ! Thank you On Thursday, August 24, 2017 at 4:26:59 PM UTC-5, Bruno Turcksin wrote: > > Jaekwong > > On Thursday, August 24, 2017 at 5:01:24 PM UTC-4, Jaekwang Kim wrote: >> >> >> <https://lh3.googleusercon

[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] Re: Relation between Solution Error Behavior and Polynomial Approximation Degree

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: Issue with convergence of iterative linear solver for system matrix in modified step-57 with no-normal flux constraints

Hi, Bruno I am thinking of similar problem. Could I ask you one thing? When you allowed slip on the top and bottom wall, did you implement any specific form of tangential velocity profile? or or just only u_y=0 (no flux) is implied? Thanks, Jaekwang On Saturday, February 17, 2018 at 4:46

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 comp

[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] Re: Use of 'VectorTools::interpolate_boundary_values' using numerical solution.

his 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] 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 fol

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 03/29

[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 dragging

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

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; GridTools::collect_periodic_faces(dof_h

[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'

[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, Jaek

[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 " 'Triangulation

[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 'ma

[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_COMPILER_DEFAULTS_TO_CXX11_OR_

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] 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, > > You h