Re: [deal.II] Re: installing dealii on a cluster

2017-10-05 Thread Anna Avdeeva
Dear Daniel, Thank you for you replies. The same error happens with the automatically configured version, but on 64%: [ 64%] Building CXX object source/CMakeFiles/deal_II.dir/dummy.cc.o Linking CXX shared library ../lib/libdeal_II.so collect2: error: ld terminated with signal 7 [Bus error], core

[deal.II] Copy example without compiling them

2017-10-05 Thread Praveen C
Dear all Is it possible to copy the deal.II example files during installation without compiling them. Since compiling examples takes some time on my laptop, I would like to avoid doing this. But I still want the example files to be installed to my deal.II installation directory. Thanks praveen -

[deal.II] Re: Vector value problem with FE_Enriched

2017-10-05 Thread Denis Davydov
Hi Edith, On Thursday, October 5, 2017 at 5:26:04 PM UTC+2, Edith Sotelo wrote: > > Hi deal ii community , > > I am trying to solve a similar problem as tutorial 29: complex valued > -Helmoholtz equation but in this case each component of my 2D vector (real, > imaginary ) is a FE-Enriched Fini

Re: [deal.II] 2D cantilever beam: reference to step-7,8,11,22

2017-10-05 Thread Wolfgang Bangerth
On 10/05/2017 12:19 AM, priyanka.nem...@iitgn.ac.in wrote: I am solving for a 2D cantiliever beam with a rectangular grid and referring to step-8. The problem is because boundary id is set to 0 by default, though I am appliying tensile force, could not obsrve any elongation towards right end.

Re: [deal.II] L2-norm gets decreased to certain value, but stays above threshold

2017-10-05 Thread Andrew McBride
I’m a bit confused by what you are asking. The entries in the residual vector corresponding to constrained dofs will converge to a non-zero number. Those corresponding to unconstrained dofs should converge to zero. As is stated in step-15 … we have to set to zero the elements of the residual ve

Re: [deal.II] L2-norm gets decreased to certain value, but stays above threshold

2017-10-05 Thread 'Maxi Miller' via deal.II User Group
Compared to the target value? Depending on the values between ~1e-6 and 2e-3, but it can be larger, too. I will check the latter function and compare it to the current solution, but I thought by using ZeroFunction as boundary condition that then all boundary_dofs are set to zero. Am Donnerstag,

Re: [deal.II] L2-norm gets decreased to certain value, but stays above threshold

2017-10-05 Thread Timo Heister
> Isn't the 'distribute'-line doing exactly that? Or is that a mistake? No, set_zero() will set the entries to zero, distribute() will set them to their correct nonzero entry. Also, how small are the errors that you are getting when it levels off? -- Timo Heister http://www.math.clemson.edu/~he

Re: [deal.II] L2-norm gets decreased to certain value, but stays above threshold

2017-10-05 Thread 'Maxi Miller' via deal.II User Group
Isn't the 'distribute'-line doing exactly that? Or is that a mistake? Am Donnerstag, 5. Oktober 2017 18:24:58 UTC+2 schrieb mac: > > I’m referring specifically to a line similar to the following where all > entries in the residual associated with constrained dofs are zeroed: > > for (unsigned int

Re: [deal.II] L2-norm gets decreased to certain value, but stays above threshold

2017-10-05 Thread Andrew McBride
I’m referring specifically to a line similar to the following where all entries in the residual associated with constrained dofs are zeroed: for (unsigned int i=0; ihttp://www.dealii.org/developer/doxygen/deal.II/classDoFHandler.html#a54e7270f2ba6206604f794114b39a2aa>(); ++i) if (boundary_dofs

Re: [deal.II] Installation: Candi on cluster (Platform)

2017-10-05 Thread 'Uwe Köcher' via deal . II User Group
Thanks for reporting. We still had some individual packages in the generic linux_cluster.platform file. I've created a pull request to remove them https://github.com/dealii/candi/pull/54 Kind regards Uwe On Thursday, 5 October 2017 17:01:00 UTC+2, Kartik Jujare wrote: > > That is exactly the

Re: [deal.II] L2-norm gets decreased to certain value, but stays above threshold

2017-10-05 Thread 'Maxi Miller' via deal.II User Group
I assume so. My approach to ensure that is: I create a ConstraintMatrix called "Newton_constraints", initialize it as following: newton_constraints.clear(); newton_constraints.reinit(solution_relevant_partitioning); DoFTools::make_hanging_node_constraints(dof_handler, newton_constraints); VectorT

Re: [deal.II] set default parameters

2017-10-05 Thread Wolfgang Bangerth
On 10/05/2017 09:39 AM, topologyd...@gmail.com wrote: But, for parameter type as patterns::list(), how can I set the value of it? Because there is no ParameterHandler::set() to set parameter which has a type of patterns::list(). If you guys know how to solve this question, please tell me, th

[deal.II] set default parameters

2017-10-05 Thread topologyding
Hi everyone, I am working on a application of dealii, and here is the question: I want to use ParameterHandler::set() function to overwrite the parameters read from .prm files or just set these parameters. Some types of parameters can be done in this way like double&, bool& cause dealii has alr

Re: [deal.II] L2-norm gets decreased to certain value, but stays above threshold

2017-10-05 Thread Andrew McBride
Are you ensuring that the entries in F(U) associated with the constrained dofs are set to zero? As shown in http://www.dealii.org/developer/doxygen/deal.II/step_15.html#MinimalSurfaceProblemcompute_residual

[deal.II] Vector value problem with FE_Enriched

2017-10-05 Thread Edith Sotelo
Hi deal ii community , I am trying to solve a similar problem as tutorial 29: complex valued -Helmoholtz equation but in this case each component of my 2D vector (real, imaginary ) is a FE-Enriched Finite element. e.g. my constructor is of the form FEsystem (fe_enriched,1,fe_enriched,1) and

Re: [deal.II] L2-norm gets decreased to certain value, but stays above threshold

2017-10-05 Thread 'Maxi Miller' via deal.II User Group
Yes, I forgot to clarify. I am calculating the l2-norm of the right hand side, i.e. which contains F(U) (when staying in the terminology of step-15). According to my knowledge that should go to zero, after that is the whole goal of the Newton iteration (if my understanding of the method is corr

Re: [deal.II] Re: installing dealii on a cluster

2017-10-05 Thread Wolfgang Bangerth
On 10/04/2017 11:11 PM, Anna Avdeeva wrote: I have tried to submit the job, but it stopped again at the same place: [ 57%] Linking CXX shared library ../lib/libdeal_II.so Attached is the error file. This looks like the linker runs out of either memory or disk space. For both, your sysadmins

Re: [deal.II] Installation: Candi on cluster (Platform)

2017-10-05 Thread Kartik Jujare
That is exactly the problem. The script tries to install zlib, bzip2 and boost even though they are switched off in the configuration file. Regards, Kartik Jujare On Wednesday, October 4, 2017 at 11:30:25 PM UTC+2, Uwe Köcher wrote: > > Ah okay, now I understand. Since the package boost (from ca

Re: [deal.II] L2-norm gets decreased to certain value, but stays above threshold

2017-10-05 Thread Wolfgang Bangerth
On 10/05/2017 08:52 AM, 'Maxi Miller' via deal.II User Group wrote: Having a nonlinear, time-dependent heat equation, I tried to implement it as in step-15 and in step-33. Everytime my l2-value goes below a certain threshold, The l2 norm of what? The solution? The residual? I make one step i

Re: [deal.II] Parallel implementation

2017-10-05 Thread Wolfgang Bangerth
On 10/05/2017 07:45 AM, Anna Avdeeva wrote: I do compute the solution at many points along the profile, and while Ex and Ey look reasonable, the Ez does not. Is it possible that there is a problem in point_value function for z component for Nedelec elements? It's possible (any software has bu

[deal.II] L2-norm gets decreased to certain value, but stays above threshold

2017-10-05 Thread 'Maxi Miller' via deal.II User Group
Having a nonlinear, time-dependent heat equation, I tried to implement it as in step-15 and in step-33. Everytime my l2-value goes below a certain threshold, I make one step in time (or I refine the mesh). At the beginning the L2-value goes down rather fast, but later it slows down, and even sto

Re: [deal.II] Re: installing dealii on a cluster

2017-10-05 Thread Timo Heister
Maybe you are running out of disk space (or some other filesystem quota). Please ask your admins for help. On Thu, Oct 5, 2017 at 9:35 AM, Daniel Arndt wrote: > Anna, > > the only other thing I could think of might be a problem with your > dependencies. > Does the problem persist if you try to ju

[deal.II] Re: Initial setup of dealII fails for Intel MPI, but not for OpenMPI

2017-10-05 Thread 'Maxi Miller' via deal.II User Group
When removing everything, and starting from a clean build-directory, I get LANG=eng cmake .. -DCMAKE_CXX_COMPILER=/share/apps/software/Compiler/intel/ 2017.1.132-GCC-6.3.0-2.27/impi/2017.1.132/bin64/mpicxx -DCMAKE_C_COMPILER= /share/apps/software/Compiler/intel/2017.1.132-GCC-6.3.0-2.27/impi/2017.1

Re: [deal.II] Parallel implementation

2017-10-05 Thread Anna Avdeeva
Dear Wolfgang, I do compute the solution at many points along the profile, and while Ex and Ey look reasonable, the Ez does not. Is it possible that there is a problem in point_value function for z component for Nedelec elements? Anna -- The deal.II project is located at http://www.dealii.or

Re: [deal.II] Parallel implementation

2017-10-05 Thread Wolfgang Bangerth
On 10/05/2017 03:05 AM, Anna Avdeeva wrote: * * 2) For the file with the solution values at the receiver locations I think I followed the approach implemented in ASPECT. So I go through all the receivers look for the point value and then check how many processors found this point and if more th

[deal.II] Re: Unit test in user project with .prm and .data input files to be read

2017-10-05 Thread Vinayak Gholap
Hello Daniel, as I understand this function is to include the .prm file input data. What I have is inside my .prm file: subsection somedata set input file = *xyz.data* end I read this input string from .prm and via my own implemented function read this *xyz.data* file. The problem occurs

Re: [deal.II] Parallel implementation

2017-10-05 Thread Anna Avdeeva
Dear Wolfgang, Thank you for your replies. I am still struggling with outputs of the solution to the file. I have two type of outputs: 1) output vtu files for each processor on the whole mesh and 2) creating a simple txt file with values of the solution at receiver locations 1) For the first

Re: [deal.II] SystemMatrix gets large fast, error or to be expected?

2017-10-05 Thread 'Maxi Miller' via deal.II User Group
Of course, having an (former) undetected bug in the equations itself also can lead to unwanted memory consumption... Thanks for the help! Am Mittwoch, 4. Oktober 2017 17:23:05 UTC+2 schrieb Wolfgang Bangerth: > > On 10/04/2017 01:03 AM, 'Maxi Miller' via deal.II User Group wrote: > > I rechecked

[deal.II] Initial setup of dealII fails for Intel MPI, but not for OpenMPI

2017-10-05 Thread 'Maxi Miller' via deal.II User Group
I tried to set up dealII with the Intel MPI compilers (using gcc 6.3.0 in the background), but already in the initial setup before I can choose any values I get the error Building CXX object CMakeFiles/cmTC_8341d.dir/testCXXCompiler.cxx.o /share/apps/software/Compiler/intel/2017.1.132-GCC-6.3.