Re: [deal.II] Re: Nedelec Elements for unstructured adaptively refined meshes

2018-05-10 16:52 GMT-04:00 Matthias Maier : > with the addition "in 2D" and for "lowest and second order" is true :-) > Yeah you are right, it's always in 3D that people have problems. Bruno -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https:/

Re: [deal.II] Re: Nedelec Elements for unstructured adaptively refined meshes

2018-05-10 16:35 GMT-04:00 Matthias Maier : > > We keep repeating that the current element only works on Cartesian > meshes, but I am not sure whether this is actually true. > Many people have reported problems see this thread for example https://groups.google.com/forum/#!searchin/dealii/nedelec%7

[deal.II] Re: Nedelec Elements for unstructured adaptively refined meshes

Artur, On Wednesday, May 9, 2018 at 8:51:37 PM UTC-4, Artur Safin wrote: > 1) Can FE_Nedelec > work > for unstructured adaptively refined meshes (for linear polynomials)? > 2) There is a new Nedelec FE class (pull 2240 >

Re: [deal.II] Re: Getting the maximum values for a vector containing multiple solutions

2018-05-04 10:42 GMT-04:00 'Maxi Miller' via deal.II User Group < dealii@googlegroups.com>: > Hmm, when doing that I only get the maximum value of the elements located > on one processor, but not for all nodes. > Once you get the maximum on each processor, use http://dealii.org/8.5.0/doxygen/deal.

Re: [deal.II] Re: Undefined Referece to `dealii::DataOut_DoFData, 2, 2>::get_dataset_names[abi:cxx11]() const'

2018-05-04 9:38 GMT-04:00 phillip mobley : > Doing it the normal way with cmake, step 6 is able to run properly without > any errors > You can do make VERBOSE=1 when you compile step-6 and you will see exactly which command is used to compile and link the tutorial. Then, you should be able to use

Re: [deal.II] Re: Undefined Referece to `dealii::DataOut_DoFData, 2, 2>::get_dataset_names[abi:cxx11]() const'

Can you run step-6 if you use cmake, i.e., if you use the recommended way? 2018-05-04 9:06 GMT-04:00 phillip mobley : > Yes, I do get the same error for step-6. > > Also, I tried the above fixes and I am still receiving the error > > On Friday, May 4, 2018 at 9:02:05 AM UTC-4,

Re: [deal.II] Re: Undefined Referece to `dealii::DataOut_DoFData, 2, 2>::get_dataset_names[abi:cxx11]() const'

2018-05-04 8:53 GMT-04:00 phillip mobley : > I as able to compile step-1 and run it within the IDE without any issues. > Maybe I am referencing something in my own code that is not being reference > in the step-1 > step-1 does not use DataOut or FE_Q. What happens with step-6? Do you have the same

[deal.II] Re: Undefined Referece to `dealii::DataOut_DoFData, 2, 2>::get_dataset_names[abi:cxx11]() const'

Philip, You can do what you are doing but this is really error prone. Two things that can help you: 1) don't move libdeal_II.so but add the path to your LD_LIBRARY_PATH: export LD_LIBRARY_PATH=/path/to/lib:$LD_LIBRARY_PATH 2) you want to use the flags in lib/cmake/deal.II/deal.IIConfig.cmake In

Re: [deal.II] confusing result when re-read a mesh in gmsh format

On Monday, April 30, 2018 at 10:45:24 AM UTC-4, zeng wrote: > > Thank you so much for your reply. > So in dealii we can't read a mesh with hanging nodes? Is there any other > way to read an adaptively refined mesh? Since the refinement happend inside > dealii, I think it should have the abilit

[deal.II] Re: Calculation fails with "AztecOO::Iterate error code -4: GMRES Hessenberg ill-conditioned"

Hi, On Friday, April 27, 2018 at 1:36:26 AM UTC-4, Maxi Miller wrote: > What can I do to find and fix the problem for that error? Or should I try > to switch to a direct solver (like UMFPACK)? > You need to find a preconditioner that works better for your problem. You can use a direct solver bu

[deal.II] Re: Candi -ubuntu 17.10.1

Yaakov, Yes this should work fine. Best, Bruno On Tuesday, April 24, 2018 at 8:09:19 PM UTC-4, feap...@gmail.com wrote: > > Dear Prof. Bangerth, dear Dr. Köcher, dear Dealii Team, > > I would like to know, In Ubuntu *17.10.1*, can I use Candi. sh during > installation without any bugs? > > >

[deal.II] Re: Interior node locations for higher order elements

Shawn, On Friday, April 20, 2018 at 2:31:28 PM UTC-4, Yuxiang Wang wrote: > > Could you please help point to me where this is documented, about the rule > of where higher order interior nodes are located? > See here http://dealii.org/8.5.0/doxygen/deal.II/classFE__Q.html and here is the relevant

[deal.II] Re: Complex PetscScalar and distribute_local_to_global

Hi, On Thursday, April 19, 2018 at 9:12:06 AM UTC-4, Juan Carlos Araujo Cabarcas wrote: > > > I have reached a few missing definitions with *PetscScalar* in the past, > and I wonder if this could be a similar error: a missing > definition/declaration? > > This should be easy to check. What happ

Re: [deal.II] Re: Getting the maximum values for a vector containing multiple solutions

On Thursday, April 19, 2018 at 5:56:54 AM UTC-4, Maxi Miller wrote: > > Assumed I call that on a LinearAlgebraTrilinos::MPI::Vector, I get the > error message > The violated condition was: > local_index >= 0 > Additional information: > You tried to access element 0 of a distributed vect

[deal.II] Re: algorithm of finding cell in a mesh by a point

Ilya On Monday, April 16, 2018 at 11:53:20 AM UTC-4, Bryukhanov Ilya wrote: > > I heard about algorithms that produce mean asymptotic O(log(N)) or at > least O(square root (N)) > To do that you need to first build a tree and then search in it. This is not what is done the current function but if

[deal.II] Re: deal.II cannot found HDF5

Peng, On Saturday, April 7, 2018 at 7:29:31 AM UTC-4, Peng Ding wrote: > > Dear group: > I have compiled the HDF5 with "--enable-parallel". but deal.II still > cannot recognize it. > How to resolve this problem. > Regards > Like Daniel said you need to give us more information. One pos

Re: [deal.II] Re: Very minor bug in doxygen documentation

2018-04-06 11:26 GMT-04:00 Pratik : > I actually am still able to see it in line 296. > > ;) > Yes, you are right. I am not sure why this happens. The similar loop in 2170 does not have this problem. I don't kno

[deal.II] Re: Very minor bug in doxygen documentation

Pratik, You are right but it looks that it has been fixed in the developer version see around line 295 here Anyway thanks for reporting the problem!

Re: [deal.II] Debug error in Kelly Error Estimate function using parallel distributed vectors

On Wednesday, April 4, 2018 at 11:45:03 PM UTC-4, Wolfgang Bangerth wrote: > I don't know whether this function is supposed to initialize ghost entries > in > the vector -- maybe one of the MatrixFree folks could help out with the > answer? > initialize_dof_vector will create a vector that has

Re: [deal.II] Sparse ILU

Phani, On Monday, April 2, 2018 at 4:59:02 PM UTC-4, Phani Motamarri wrote: > > I am trying to modify the function "SparseILU::initialize" inside > sparse_ilu.templates.h to call the mkl function "dcsrilu0" which requires > the use of some of the data like ia, ja accessed inside "initialize" >

[deal.II] Re: dealii 8.5.1 installation issue (spack, docker, debian)

Marquito, On Sunday, March 11, 2018 at 9:27:40 AM UTC-4, Marquito Forrest wrote: > > > In summary the situation is: > - Debian testing image (to get recent gcc) > - Miniconda Python 3.6 (my own Dockerfile, cf. file 1) > - Latest spack from github (another Dockerfile derived from file 1, cf. > fil

[deal.II] Re: Getting the maximum values for a vector containing multiple solutions

Hi, On Friday, March 9, 2018 at 3:58:22 AM UTC-5, Maxi Miller wrote: > > I have a solution vector of the type LinearAlgebraTrilinos::MPI::Vector, > containing solutions for multiple equations. Now I would like to know the > maximum value for each of those solutions. When simply calling > vector

Re: [deal.II] Re: FE w/ Gauss Point DoFs

Jonathan, 2018-03-07 10:43 GMT-05:00 Jonathan Russ : > > Thank you for your fast reply. I'm sorry I'm not doing a great job of > explaining this. Here is an example from elasticity that illustrates what I > would like: > In this case there are displacement degrees of freedom at the nodes but > th

[deal.II] Re: FE w/ Gauss Point DoFs

Jonathan, On Wednesday, March 7, 2018 at 10:20:44 AM UTC-5, Jonathan Russ wrote: > > Is there a FiniteElement in deal.II that allows the addition of degrees of > freedom at gauss points so that the DoFHandler object also manages them? > For example, I would like to have just a standard Lagrangia

[deal.II] Re: Error in installing deal.ii with spack on linux openSuse

Alberto, On Monday, February 19, 2018 at 9:56:25 AM UTC-5, Alberto Salvadori wrote: > > dealii requires cmake version :3.9.99, but spec asked for 3.10.1 > deal.II 8.5 does not support cmake 3.10 so you need to force spack to use a lower version. Something like this should work spack dealii ^cmake

Re: [deal.II] Re: Problem compiling deal_ii with PETSc

Daipayan, 2018-02-15 11:17 GMT-05:00 Daipayan Sarkar : > Thanks for your help. I followed just as you mentioned in your reply and > now it works. For future, if I want to build any other external library > (e.g. Trillions) do I need to use cmake in the build directly to make a new > executable li

[deal.II] Re: Problem compiling deal_ii with PETSc

Hi, On Tuesday, February 13, 2018 at 8:55:34 PM UTC-5, Daipayan Sarkar wrote: > > At ~ 48% I get the following error message: > > *error: **"PETSc was configured with MPICH but now appears to be > compiling using a non-MPICH mpi.h"* > The problem is that you are mixing two different versions of

[deal.II] Re: Use CUDA or C++17 when compiling deal.II?

Hi, On Tuesday, February 13, 2018 at 7:43:36 AM UTC-5, Maxi Miller wrote: > > When configuring deal.II, I can either choose to use C++17 or CUDA. Which > advantages would I get if I use one over the other? > It depends if you plan on using CUDA or not. Right now, the support for CUDA in deal is

Re: [deal.II] Re: Heat equation calculation with narrow, high-intensity heat source

2018-01-25 9:05 GMT-05:00 'Maxi Miller' via deal.II User Group < dealii@googlegroups.com>: > > Which preconditioner could work better here? > > I don't know. You should look in the literature. Best, Bruno -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum option

[deal.II] Re: Heat equation calculation with narrow, high-intensity heat source

Hi, On Thursday, January 25, 2018 at 4:00:45 AM UTC-5, Maxi Miller wrote: > > I am trying to calculate a heat equation using a similar approach as in > step-33 (auto-diff). Now my heat source is very narrow and very intense, > resulting in a high gradient (1e6 within one or two cells), thus cras

Re: [deal.II] Re: Problem with the use of FE_Nedelec with a mesh file generated from Gmsh

Jianan, 2018-01-24 21:23 GMT-05:00 Jianan Zhang : > > Just want to inform you and the dealii users who may be interested that my > problem is solved perfectly by using Ross' new implementation FE_NedelecSZ. > The results given by brick elements and hexes are quite similar. Thanks a > lot for the

[deal.II] Re: installation error

Juan Carlos On Tuesday, January 23, 2018 at 3:12:20 PM UTC-5, Juan Carlos Araujo Cabarcas wrote: > > [ 50%] Building CXX object > source/fe/CMakeFiles/obj_fe_debug.dir/fe_poly.cc.o > mypath/dealii/source/lac/scalapack.cc: In member function ‘void > dealii::ScaLAPACKMatrix::copy_to(dealii::ScaLA

[deal.II] Re: Run function on one node, make result accessible for all nodes

Hi, On Tuesday, January 23, 2018 at 7:53:16 AM UTC-5, Maxi Miller wrote: > But now it looks like as if only the first node gets the result of the > calculations, but the others do not, instead defaulting to the default > values of the calculation function when not initialized. Is there a way I

Re: [deal.II] Re: Parallel DoFRenumbering does not work on mac

Jie, 2018-01-18 14:31 GMT-05:00 Jie Cheng : > Fatal error in MPI_Init_thread: Other MPI error, error stack: > MPIR_Init_thread(474): > MPID_Init(152)...: channel initialization failed > MPID_Init(426)...: PMI_Get_appnum returned -1 > [cli_0]: write_line error; fd=6 buf=:cmd=abort exitcode

[deal.II] Re: Parallel DoFRenumbering does not work on mac

Jie, On Thursday, January 18, 2018 at 11:17:55 AM UTC-5, Jie Cheng wrote: > > I found that DoFRenumbering with distributed dof_handler does not work on > my mac (Mac OS X 10.13.2). My dealii is built with PETSc and Step-40 runs > fine. But once I added DoFRenumbering::Cuthill_McKee(dof_handler)

[deal.II] Re: Writing parallel distributed vectors

Krishnendu, On Saturday, January 13, 2018 at 7:29:26 PM UTC-5, Krishnendu Ghosh wrote: > I would like to write dealii parallel distributed vectors on files in > order to implement restart capability. Please let me know what will be the > best way to do that. > Take a look here https://groups.

[deal.II] Re: Undefined reference to p::d::SolutionTransfer with hp::DoFHandler

Claire, On Friday, January 12, 2018 at 1:55:31 PM UTC-5, Claire wrote: > > I tried to implement the new feature made available recently in the > development version of dealii that allows to use a p::d::Triangulation > together with a hp::DoFHandler. > This is not working yet. To quote Wolfgang

Re: [deal.II] Re: Set inner boundary conditions

Qing, 2018-01-12 13:00 GMT-05:00 Qing Yin : > OK, let's see the simple example in step-3. From a physical viewpoint, > this is a steady-state heat equation without a heat source. Then we hope > there is a particular line area inside, where the temperature is, for > example, 1. In other place, the

[deal.II] Re: Set inner boundary conditions

Qing, What do you want to do? It's hard to help you if you don't tell us why you want to set an inner boundary conditions. Best, Bruno On Thursday, January 11, 2018 at 6:30:11 PM UTC-5, Qing Yin wrote: > > Dear deal.II community, > > I am wondering if it is possible to set inner boundary condi

[deal.II] Re: dof renumbering on distributed meshes

Marek, Take a look here https://github.com/dealii/dealii/issues/5407 and here https://github.com/dealii/dealii/pull/5475#commitcomment-25928779 I think this is similar to what you are trying to do. Note that you need the development version of deal.II Best, Bruno On Wednesday, January 10, 20

Re: [deal.II] Re: Deal.II with cmake 3.10 and boost-1.66

2018-01-10 3:20 GMT-05:00 'Maxi Miller' via deal.II User Group < dealii@googlegroups.com>: > I definitely implemented that pull-request, but I still get the same > error. Furthermore I do not see those switches coming up in cmake-gui, even > in advanced mode. > When looking into the pull request,

Re: [deal.II] Re: Deal.II with cmake 3.10 and boost-1.66

2018-01-09 9:44 GMT-05:00 'Maxi Miller' via deal.II User Group < dealii@googlegroups.com>: > How can I implement that when using cmake-gui? Then the command is not > working. > What do you mean? It is a variable created by cmake when you use FindMPI. I can see it using the gui when I toggle the ad

Re: [deal.II] Re: Deal.II with cmake 3.10 and boost-1.66

I meant path to mpicxx not mpicc 2018-01-09 8:49 GMT-05:00 Bruno Turcksin : > > 2018-01-09 8:31 GMT-05:00 'Maxi Miller' via deal.II User Group < > dealii@googlegroups.com>: > >> Thanks, I already found that. The bigger problem for me atm. is the >> coup

Re: [deal.II] Re: Deal.II with cmake 3.10 and boost-1.66

2018-01-09 8:31 GMT-05:00 'Maxi Miller' via deal.II User Group < dealii@googlegroups.com>: > Thanks, I already found that. The bigger problem for me atm. is the > coupling between MPI and deal.II. Is there a possible fix for me, too? > Thanks! > Hmmm are you sure you implemented the PR correctly?

[deal.II] Re: Deal.II with cmake 3.10 and boost-1.66

Hi, This is is not a real bug. The problem is that cmake check that it knows the version of boost that you are using and so every time that boost creates a new release, the cmake guys need to update the list of supported version. Since boost 1.66 was release after cmake 3.10, it is not in the

Re: [deal.II] Re: Tutorial fail to run in my enviroment and i need help ....

2018-01-05 10:11 GMT-05:00 Hungi : > /usr/local/opt/suite-sparse/lib/ , I do not have > libsuitesparseconfig.4.5.5.dylib but libsuitesparseconfig.5.1.0.dylib and > libsuitesparseconfig.5.dylib , libsuitesparseconfig.dylib instead ? > I don't know how you installed deal.II but you did something wr

[deal.II] Re: Tutorial fail to run in my enviroment and i need help ....

Hi, On Friday, January 5, 2018 at 9:29:45 AM UTC-5, Hungi wrote: > > after I updateing sundials through Homebrew : failure were changed : > > bash-3.2$ ./step-1 > > dyld: Library not loaded: > /usr/local/opt/suite-sparse/lib/libsuitesparseconfig.4.5.5.dylib > > Referenced from: > /Users/hicho

Re: [deal.II] Re: Problem with the use of FE_Nedelec with a mesh file generated from Gmsh

Jianan 2017-12-28 10:40 GMT-05:00 Jianan Zhang : > > So, can you help me contact the people working on the PR to try this out > to see if my code works with the new implementation? > If you just want to try the new implementation you can take Ross' branch

[deal.II] Re: Ask for deal.II method to write a template class in separate files

Qing On Friday, December 22, 2017 at 1:44:29 PM UTC-5, Qing Yin wrote: > > I am asking for the method to declare and define a dimension-independent > template class in separate files . Because from all the book/tutorial on > templates I have read, it recommends that template declaration and >

Re: [deal.II] Re: Problem with the use of FE_Nedelec with a mesh file generated from Gmsh

Jianan, 2017-12-20 17:21 GMT-05:00 Jianan Zhang : > Thanks for the prompt reply. I will try the implementation as you > suggested, hope I can make it work with it. Otherwise, I will just go with > brick elements, although the EM structure I can model may be limited. > According to the comments in

Re: [deal.II] Re: Problem with the use of FE_Nedelec with a mesh file generated from Gmsh

Jianan, 2017-12-20 14:54 GMT-05:00 Jianan Zhang : > > What I am really expecting is that someone who has some experience of > using FE_Nedelec to solve EM problems can give me some general advice. For > example, can we use FE_Nedelec on a mesh imported from other mesh > generators instead of the

[deal.II] Re: Problem with the use of FE_Nedelec with a mesh file generated from Gmsh

Jianan, On Wednesday, December 20, 2017 at 12:56:11 AM UTC-5, Jianan Zhang wrote: > > I've spent like ten days finding the reason, and I am almost 100% sure > that my exact solution is given correctly and the assembly is correct. The > only thing I changed is the exact solution, and the differe

Re: [deal.II] Re: Moving Mesh when using Fe_Q(p) p>1

Hi, 2017-12-19 15:59 GMT-05:00 RAJAT ARORA : > > So, going by what you described, if x2 is unknown, how can the mapping > class calculate the information? > > Is it because it uses the manifold information (which is a straight > boundary by default) attached to it to get the nodal coordinates of

[deal.II] Re: Why does the L2-norm increase when the number of nodes becomes too large?

Hi, On Tuesday, December 19, 2017 at 2:58:20 PM UTC-5, jie liu wrote: > > > This phenomenon is not desired. Different solvers are adopted, but there > is no improvement on it. Can anybody give some clues? Thank you. > You have a bug in your code. What you should see is that the error decreases a

Re: [deal.II] Re: Cmake-MPI

g trilinos 12-10-1.* > Many thanks for your help again! > > Warm regards, > Yaakov > > On Thursday, December 14, 2017 at 5:29:31 PM UTC+1, Bruno Turcksin wrote: >> >> Yaakov, >> >> 2017-12-14 11:23 GMT-05:00 : >> >>> *Error: Error: p4e

Re: [deal.II] Re: Cmake-MPI

Yaakov, 2017-12-14 11:23 GMT-05:00 : > *Error: Error: p4est couldn't find zlib* > > how can I solve this problem in candi? > sudo apt install zlib1g-dev Best, Bruno -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/for

Re: [deal.II] Re: Cmake-MPI

; $ P4EST_DIR="..." cmake <...> > $ cmake -DP4EST_DIR="..." <...> > > or set the relevant variables by hand in ccmake. > > Call Stack (most recent call first): > cmake/macros/macro_configure_feature.cmake:269 (FEATURE_ERROR_MESSAGE) >

Re: [deal.II] Re: Cmake-MPI

Yaakov, 2017-12-11 11:00 GMT-05:00 : > > for mpi installation. I only use command : *sudo apt-get install mpich* > Yet you have openmpi installed: > > -- Found MPI_MPI_H > -- MPI_VERSION: 3.1 > -- MPI_LIBRARIES: /usr/lib/x86_64-linux-gnu/openmpi/lib/libmpi_cxx.so;/ > usr/lib/libmpi.so;/usr/l

[deal.II] Re: independent triangulations on different mpi processes

Jose, On Monday, December 11, 2017 at 5:36:10 AM UTC-5, Jose Javier Muñoz Criollo wrote: > > My current problem involves the use of a relatively big main domain that > is distributed accross a number of mpi process. The solution on this domain > is coupled with the solutions on smaller domains.

[deal.II] Re: Moving Mesh when using Fe_Q(p) p>1

Hi, On Sunday, December 10, 2017 at 9:55:49 PM UTC-5, RAJAT ARORA wrote: > My question is: > Is the mesh movement independent of the order of the polynomial used to > interpolate the solution? > > What is confusing me is that I learned that there are nodes (other than > vertices) in the element

[deal.II] Re: Parallelization of two step method

Mustafa, On Sunday, December 10, 2017 at 4:06:21 PM UTC-5, Mustafa Aggul wrote: > > Here, run_FS(n_time_steps) is independent from CS, so I can run it a > little bit earlier but parallel to run_CS(n_time_steps). I think I can take > care of running it earlier part with mutex, but I couldn't crea

[deal.II] Re: Heat equation (step-26): Negative values with small time step

Marc, On Thursday, December 7, 2017 at 11:52:39 AM UTC-5, Marc Fehling wrote: > I stumbled over some interesting behavior of the heat equation from > step-26. If I reduce the time step to a smaller value, let's say to 1e-6, I > observe negative values for the solution near the sources (where gr

[deal.II] Re: Some cells have negative volume error

Narendra, On Tuesday, December 5, 2017 at 11:27:15 AM UTC-5, Narendra Nanal wrote: > > > When the input file is imported in abaqus it is working fine. Please see > attached pictures.I am not sure if something is wrong with the mesh. Kindly > guide me on how to resolve this error. Thanks in advan

[deal.II] Re: Cmake-MPI

Yaakov, On Monday, December 4, 2017 at 2:47:00 AM UTC-5, feap...@gmail.com wrote: > [ 77%] Linking CXX executable ../bin/step-8.release > ../lib/libdeal_II.so.9.0.0-pre: error: undefined reference to > 'Epetra_MpiComm::Epetra_MpiComm(int)' > There is a problem with Trilinos. deal.II can't link a

Re: [deal.II] Re: mapping_collection and step-27

On Friday, December 1, 2017 at 1:47:44 PM UTC-5, Juan Carlos Araujo Cabarcas wrote: > I must admit that I am a bit lost here. However, I tried your suggestion > and used:FunctionMap::type, what gave the errors that > I attach in errors.txt. > > You have attached a file but it only contain

Re: [deal.II] Re: Error during p4est installation

Yaakov, 2017-12-01 8:26 GMT-05:00 : > one is defined order p4est > new one is p4est-build > p4est: directory where you wanted to install p4est p4est-build: directory where p4est is compiled. This can be removed once p4est is installed. It is common when you install a library to have three directo

Re: [deal.II] Re: Spurious eigenvalues

2017-12-01 8:50 GMT-05:00 Timo Heister : > 3. Did we rip out the support for removing the constrained entries > from the matrix completely? That was my plan at first using ConstraintMatix::condense but according to the documentation, the constrained entries are not removed. Best, Bruno -- The

Re: [deal.II] Re: Spurious eigenvalues

DoFs. > Then being Krylov-based method it should form subspaces {x, Ax, A^2x,...} > orthogonal to those constrained DoFs, so > you should not get any issues. > > Cheers, > Denis. > > > On Thursday, November 30, 2017 at 4:24:05 PM UTC+1, Bruno Turcksin wrote: >> >>

[deal.II] Spurious eigenvalues

Hi all, In step-36, there is an explanation on how Dirichlet boundary conditions introduce spurious eigenvalues because some dofs are constrained. However, there is no mention of hanging nodes. So I am wondering if I can treat them as shown for the Dirichlet boundary, i.e, the only difference

Re: [deal.II] Re: Conjugate Gradient for Schur complement, serial vs parallel discrepancy in solution and effect of tolerance.

Dimitris, 2017-11-27 17:10 GMT-05:00 Dimitris Ntogkas : > Thanks for your quick response! You are right about the l2 vs max norm. > However, the error is 1e-4 in the l2 norm too. Just a clarification to make > sure I understand your response. I was indeed thinking of the condition > number, that'

[deal.II] Re: Conjugate Gradient for Schur complement, serial vs parallel discrepancy in solution and effect of tolerance.

Dimitris, The same implementation of CG will give you (slightly) different results in serial an parallel because the round-off errors will be different. This round-off errors will be amplified if you have a large condition number (see https://en.wikipedia.org/wiki/Condition_number). So if you

[deal.II] Re: Isogeometric analysis

Yaakov, This is not my area of expertise but I guess you will want to take a look at http://dealii.org/8.5.0/doxygen/deal.II/step_54.html Best, Bruno On Tuesday, November 21, 2017 at 12:17:25 PM UTC-5, feap...@gmail.com wrote: > > Dear Deal.II Team & User, > > I would like to ask a question:

[deal.II] Re: Questions related to FE_Nothing

Hi, On Sunday, November 19, 2017 at 9:58:25 PM UTC-5, RAJAT ARORA wrote: > > > But I am still not sure about the answers to the other questions > (mentioned below). I will appreciate if anyone can help me with these. > > 1. Can Parallel::shared::triangulation be used with hp::dofhandler? I want

[deal.II] Re: what does it mean “the neighbor is coarser, but still has children which are finer than our current cell” in step 30？

Hi, Basically, what happen is that if you ask for the right neighbor of the cell # you will get the coarse right cell. The reason is that when you ask for the neighbor, you get a cell that has the same level of refinement or that is coarser. Now, when you look at cell # you could think that the

Re: [deal.II] Re: Read in tethex grid

Felix, I am glad that you found the problem. It's good to know that Tethex works correctly if you use double. Best, Bruno On Nov 6, 2017 6:20 AM, "Felix" wrote: The above code is correct. Tethex works fine if the values in your "iso.msh" are of type double. My used data is of type and so I c

Re: [deal.II] Re: Read in tethex grid

Felix, 2017-11-03 6:48 GMT-04:00 Felix : > For a refinement ref = 1 of the dealii grid (see 1.) it worked. But for > ref > 1 the above mentioned error happens. I do not really understand how > this not correctly oriented cells are created since I thought that tethex > takes care of that. Do you h

[deal.II] Re: VectorTools::point_value() throw exception point not locally owned

Yiyang, On Tuesday, October 31, 2017 at 10:56:00 AM UTC-4, Yiyang Zhang wrote: > The above code seems to output the "ERROR" statement in this code, which > means the GridTools::find_active_cell_around_point() will return a > different cell. > It is understandable that to return the cell that a

Re: [deal.II] Re: Read in tethex grid

Felix, 2017-10-31 7:07 GMT-04:00 Felix : > > However, what I would try is to use gmsh to directly produce a mesh of >> quadrangles. >> > > Do you mean that I should avoid this "iso.msh" (consist of triangles) and > directly construct the isosurface consisting of quadrangle and save it as > .msh?

[deal.II] Re: Read in tethex grid

Felix, On Monday, October 30, 2017 at 7:32:46 AM UTC-4, Felix wrote: > > > The name and call sequence of the exception was: > > ExcNonOrientableTriangulation() > > Additional Information: > > (none) > > > Has anyone an idea how to fix this? > The problem is that there is one or more cells tha

Re: [deal.II] Non-homogeneus boundary conditions with-matrix-free not working

Martin, Could you add this explanation to the documentation and/or to a tutorial. I think mentioning it step-37 would be great. Best, Bruno On Friday, October 27, 2017 at 3:57:09 AM UTC-4, Martin Kronbichler wrote: > > Dear Michal, > > This is expected: the matrix-free operator evaluation cann

[deal.II] Re: Trilinos::MPI::Vector and inhomogeneous constraints

Yiyang, On Wednesday, October 25, 2017 at 8:07:46 PM UTC-4, Yiyang Zhang wrote: > locally_relevant_sln_with_ghost_new = locally_relevant_sln_with_ghost_old; > //This and the following one are rarely seen in the tutorials. Does this > mean the vector is exactly copied? i.e. both locally_owned_do

[deal.II] Re: Installation error

Lev, How many processors did you use when you compiled deal.II? The compiler may run out of memory and hit an internal error if you use too many processor. Can you try again using only one processor. Best, Bruno On Monday, October 23, 2017 at 9:07:40 PM UTC-4, Lev Karatun wrote: > > Hi everyo

Re: [deal.II] Re: mesh generation in step-6 in parallel::distributed::Triangulation and Trilinos

Yiang, 2017-10-19 13:53 GMT-04:00 Yiyang Zhang : > > i) I checked this page. I am a little confused. As soon as > signal.post_refinement.connect() is called, I assume set_boundary_ids() > is always attached to the Triangulation object, and it will always refresh > the boundary IDs when a refineme

[deal.II] Re: mesh generation in step-6 in parallel::distributed::Triangulation and Trilinos

Hi, On Wednesday, October 18, 2017 at 4:42:03 PM UTC-4, Yiyang Zhang wrote: > > i) the signal.post_refinement_connect() need to be set only when the grid > is generated initially, then throughout the problem (even if we have AMR in > the code), we don't need to set it again. Right? > No you stil

Re: [deal.II] Re: Using values from the time step before while solving a heat equation fails

2017-10-18 8:54 GMT-04:00 'Maxi Miller' via deal.II User Group < dealii@googlegroups.com>: > I assume that my old values are correct, after I get a reasonable result > with them. > It doesn't mean it's right though. It maybe be slightly wrong and then, the error could be amplified when you are usi

Re: [deal.II] Re: SolutionTranfer with PETScWrappers::MPI::Vector

> #4 ./linear: Step27::LaplaceProblem<2>::run() > #5 ./linear: main > > > > > Il giorno mercoledì 18 ottobre 2017 14:20:52 UTC+2, Bruno Turcksin ha > scritto: >> >> Carlo, >> >> what i

[deal.II] Re: SolutionTranfer with PETScWrappers::MPI::Vector

Carlo, what is the error that you get? The program crashes is too vague. Best, Bruno On Wednesday, October 18, 2017 at 6:46:11 AM UTC-4, Carlo Marcati wrote: > > > Hi, > I am trying to use a SolutionTrasfer object with > PETScWrappers::MPI::Vector. When trying to do > solution_trasfer.interp

[deal.II] Re: Using values from the time step before while solving a heat equation fails

Hi, On Wednesday, October 18, 2017 at 3:55:05 AM UTC-4, Maxi Miller wrote: > Apparently my old values are correct, but my old gradients are not. Thus I > am wondering if i getting the gradient values from the old solution in an > incorrect way, or is it correct? > Well did you check the value

[deal.II] Re: Error when writing files using write_vtu_header

Hi, It could be a problem with the file system, it is hard to tell. The first thing you need to do is to be able to reproduce the error consistently. Try to change the size of the problem and the number of processor that you are using and see if the error happens more often (ideally every time)

Re: [deal.II] Reflect filter in ParaView

Chenchen, here is the function you want to use to do what Jean-Paul advise you to do http://dealii.org/8.5.0/doxygen/deal.II/classDataOut.html#a5eb51872b8736849bb7e8d2007fae086 You need to choose n_subdivisions > 1. Best, Bruno On Thursday, September 28, 2017 at 4:07:59 PM UTC-4, Jean-Paul P

[deal.II] Re: Parallel distributed triangulation: informations on part of the mesh which are not owned

Giovanni, It looks like you are doing a lot of things similar to what we do in this project https://github.com/ORNL-CEES/DataTransferKit In particular, one of our subroutine allows the user to ask for the value of the field at arbitrary points. We need our code to use GPUs, so we based our algo

Re: [deal.II] Re: Refining a cylinder only in x

2017-09-18 10:32 GMT-04:00 Lucas Campos : > Also, here is the new program, with the added if and > prepare_coarsening_and_refinement(). You will find the interesting part > between lines 85 and 104. > > You want to use prepare_coarsening_and_refinement after you have set the flags not before. On

[deal.II] Re: Refining a cylinder only in x

Lucas, On Monday, September 18, 2017 at 3:50:20 AM UTC-4, Lucas Campos wrote: > > I am trying to refine a cylinder in the x direction only, but it seems > that when I use > > > cell->set_refine_flag(RefinementCase<3>::cut_x); > > the whole cylinder is refined. Indeed, if I run the above line for

Re: [deal.II] Re: about deal.II with CUDA C programming acceleration

2017-09-17 1:22 GMT-04:00 Denis Davydov : > > I though you could also get binaries (nvcc) and the rest by compiling > CUDA, i.e. https://developer.nvidia.com/compute/cuda/8.0/Prod2/ > local_installers/cuda_8.0.61_375.26_linux-run > but of course one can also set it as an externally provided packag

Re: [deal.II] Re: about deal.II with CUDA C programming acceleration

Denis, 2017-09-16 18:30 GMT-04:00 Denis Davydov : > @Bruno I don't any experience with cuda, but do we need anything special > but > $spack install cuda > to use it? > > I wouldn't use spack to install cuda because there is nothing to compile. You need binaries and you need drivers that depends

Re: [deal.II] Re: about deal.II with CUDA C programming acceleration

Chih-Che, 2017-09-15 19:42 GMT-04:00 Chih-Che Chueh : > > Thanks for the reminder and helps. Indeed, our system administrator has > trouble with installing the newer version of GNU C compiler. Would you > provide detailed instructions on how to installing the GNU C compiler? > Thanks! > > The sys

Re: [deal.II] Re: about deal.II with CUDA C programming acceleration

Chih-Che, Did you get everything to work or you gave up? You are the first non-developer trying to install deal.II with CUDA, so I am really interested to know which problems you are encountering. Best, Bruno 2017-08-23 20:49 GMT-04:00 Bruno Turcksin : > Chih-Che, > > 2017-08-23 20:

[deal.II] Re: No namespace named 'distributed' in namespace 'dealii::parallel';

Chuks, it's hard to say what's going on with so little information but you should probably check your includes. Do the tutorials that are using trilinos and MPI work for you? What is the full error message? Best, Bruno On Thursday, September 14, 2017 at 1:50:22 PM UTC-4, Chuks wrote: > > Hi,

[deal.II] Re: Linking error when compiling with Trilinos

Hi, You have several problems. On Wednesday, September 13, 2017 at 4:58:10 AM UTC-4, Maxi Miller wrote: > > ~/local_opt/trilinos/lib/libmuelu-adapters.so: error: undefined reference > to 'dmumps_c' > You have compile deal or trilinos with MUMPS support but the linker does not find it /petsc/l

Re: [deal.II] Trouble getting UMFPACK to work

2017-09-12 9:49 GMT-04:00 曾元圆 <2012zg...@gmail.com>: > If I want to run all the 57 steps, should I install all the mentioned > libraries? How can I know in advance which library I may need? Perhaps now > I don't need PETSc, but maybe someday in the future I will need. > No you don't. You can use c

<    1   2   3   4   5   6   7   8   >