Re: [deal.II] Different shape representations with manifolds on the same triangulation

2020-01-20 Thread Juan Carlos Araujo Cabarcas
equired during the assembly process. > > In summary, I can’t think of anything useful to help reduce the amount of > computation effort during the assembly process for the example you > previously described. I guess the logical question for me to ask is, why do > you want t

Re: [deal.II] Different shape representations with manifolds on the same triangulation

2020-01-20 Thread Juan Carlos Araujo Cabarcas
seful to help reduce the amount of > computation effort during the assembly process for the example you > previously described. I guess the logical question for me to ask is, why do > you want to do this in the first place? Is the assembly definitely a > bottle-neck for you? &

Re: [deal.II] Different shape representations with manifolds on the same triangulation

2020-01-20 Thread Juan Carlos Araujo Cabarcas
> would be required during the assembly process. > > In summary, I can’t think of anything useful to help reduce the amount of > computation effort during the assembly process for the example you > previously described. I guess the logical question for me to ask is, why do >

Re: [deal.II] Different shape representations with manifolds on the same triangulation

2020-01-16 Thread Juan Carlos Araujo Cabarcas
elf as part of a > pre-processing step? The GridTools::laplace_transform() > <https://www.dealii.org/current/doxygen/deal.II/namespaceGridTools.html#a7ed2aaa1aea3ac22b1e1807ce6d0b5f3> > function > could be helpful for this purpose. This, though, seems less elegant than > the f

[deal.II] Different shape representations with manifolds on the same triangulation

2020-01-15 Thread Juan Carlos Araujo Cabarcas
Dear all, I would like your guidance on how to perform the assembly of different shape representations on the same triangulation and on the same loop through cells. Let me try to explain a bit more. I have designed a grid containing two concentric squares (or circles). Additionally, I have a

Re: [deal.II] Re: optimization, SolverBFGS

2019-10-28 Thread Juan Carlos Araujo Cabarcas
Thanks for the help, now it is running nicely! On Monday, 28 October 2019 15:05:19 UTC+1, Bruno Turcksin wrote: > > There is a small bug in the lambda: > > Le lun. 28 oct. 2019 à 09:36, Juan Carlos Araujo Cabarcas > > a écrit : > > auto my_lambd

Re: [deal.II] Re: optimization, SolverBFGS

2019-10-28 Thread Juan Carlos Araujo Cabarcas
28, 2019 at 1:24:22 PM UTC+1, Bruno Turcksin wrote: > > Here is a simple example: > https://wandbox.org/permlink/fbcldo9PIiHl6tfI In your case, try > something like > > auto my_lambda = [](Vector a, Vectorb) {return > polar.Callback(a,b);}; > opt.solve ( my_lambda, X ); &

[deal.II] Re: optimization, SolverBFGS

2019-10-28 Thread Juan Carlos Araujo Cabarcas
y-bind-a-member-function-to-an-stdfunction-in-visual-studio-11 > > I advise using a lambda like the second reply suggests > > Best, > > Bruno > > On Friday, October 25, 2019 at 10:53:08 AM UTC-4, Juan Carlos Araujo > Cabarcas wrote: >> >> Dear all, >>

[deal.II] optimization, SolverBFGS

2019-10-25 Thread Juan Carlos Araujo Cabarcas
Dear all, I would like to use SolverBFGS from deal.II/optimization/ in my project, and I try to follow the documentation in https://www.dealii.org/current/doxygen/deal.II/classSolverBFGS.html First, I would like to mention that the documentation would really improve by adding a minimal

Re: [deal.II] How to manually create sparsity pattern for PETSc sparsity matrix in parallel

2019-04-03 Thread Juan Carlos Araujo Cabarcas
Hey Pai Liu, At some point in my research I had to do something similar to what you want to do. In the application of the so-called Dirichlet-to-Neumann maps for Helmholtz problems, one has to assemble a dense block corresponding to the boundary degrees of freedom, so I had to include extra

Re: [deal.II] Re: shape_value, shape_grad in arbitrary points within quad that are not quadrature points

2018-08-17 Thread Juan Carlos Araujo Cabarcas
Thanks a lot for the pointers! El viernes, 17 de agosto de 2018, 22:57:24 (UTC+2), Wolfgang Bangerth escribió: > > On 08/17/2018 02:56 PM, Juan Carlos Araujo Cabarcas wrote: > > Hej, thanks for your prompt answer! > > Yes, your interpretation of my poor description is

[deal.II] Re: shape_value, shape_grad in arbitrary points within quad that are not quadrature points

2018-08-17 Thread Juan Carlos Araujo Cabarcas
Hej, thanks for your prompt answer! Yes, your interpretation of my poor description is correct. So then I just need to create a new FEValues within the cell loop containing the points where I want to evaluate the shape functions, right? How do I add this points? Also if the points are in the

Re: [deal.II] Re: installation error

2018-01-26 Thread Juan Carlos Araujo Cabarcas
es, 25 de enero de 2018, 12:29:53 (UTC-5), Wolfgang Bangerth escribió: > > On 01/24/2018 07:34 AM, Juan Carlos Araujo Cabarcas wrote: > > Please find the file: detailed.log attached. > > > > El martes, 23 de enero de 2018, 17:02:14 (UTC-5), Wolfgang Bangerth > > e

Re: [deal.II] Re: installation error

2018-01-26 Thread Juan Carlos Araujo Cabarcas
(UTC-5), Wolfgang Bangerth escribió: > > On 01/24/2018 07:34 AM, Juan Carlos Araujo Cabarcas wrote: > > Please find the file: detailed.log attached. > > > > El martes, 23 de enero de 2018, 17:02:14 (UTC-5), Wolfgang Bangerth > > escribió: > > > >

Re: [deal.II] Re: installation error

2018-01-24 Thread Juan Carlos Araujo Cabarcas
Please find the file: detailed.log attached. El martes, 23 de enero de 2018, 17:02:14 (UTC-5), Wolfgang Bangerth escribió: > > On 01/23/2018 02:13 PM, Bruno Turcksin wrote: > > > > mypath/dealii/source/lac/scalapack.cc:243:91: error: there are no > > arguments to

[deal.II] installation error

2018-01-23 Thread Juan Carlos Araujo Cabarcas
Dear all, I am trying to install deal.II from the GIT repository with the following features: petsc_ver='3.6.0'; trilinos_ver='12.4.2'; git clone https://github.com/dealii/dealii.git dealii cmake \ -DTRILINOS_DIR=${install_dir}/trilinos-${trilinos_ver}

[deal.II] Re: Error with Transfinite Interpolation

2018-01-18 Thread Juan Carlos Araujo Cabarcas
; converge. > > I got rid of the cube assumption. see > https://github.com/dealii/dealii/pull/5761 for more information. > > Do you mind if I use a modified version of your code for a new test in > deal.II? > > Thanks, > David Wells > > On Wednesday, January 17, 20

[deal.II] Error with Transfinite Interpolation

2018-01-17 Thread Juan Carlos Araujo Cabarcas
Dear all, Since the inclusion of Transfinite interpolation, I have been successful on working with this powerful technique in my research. I had coded a mesh implementing concentric circles, where the inner most is shifted a small distance s. All concentric circles are labeled 100+i, where i

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

2017-12-08 Thread Juan Carlos Araujo Cabarcas
Thanks again, Now it passes the *interpolate_boundary_values *call. However, it is conflicting now with: MatrixTools::apply_boundary_values (bval, system_matrix,

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

2017-12-06 Thread Juan Carlos Araujo Cabarcas
;, dealii::ComponentMask const&)' collect2: error: ld returned 1 exit status make[2]: *** [pFEM] Error 1 make[1]: *** [CMakeFiles/pFEM.dir/all] Error 2 make: *** [all] Error 2 Thanks in advance, Juan Carlos Araújo El martes, 5 de diciembre de 2017, 23:47:22 (UTC-5), Wolfgang Bangert

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

2017-12-05 Thread Juan Carlos Araujo Cabarcas
et El lunes, 4 de diciembre de 2017, 13:09:28 (UTC-5), Wolfgang Bangerth escribió: > > On 12/01/2017 08:36 PM, Juan Carlos Araujo Cabarcas wrote: > > > > path/petsc_eigs/pFEM.cc:370: error: undefined reference to 'void > > dealii::VectorTools::interpolate_boundary_values&l

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

2017-12-01 Thread Juan Carlos Araujo Cabarcas
locator<std::pair > > >&, dealii::ComponentMask const&)' collect2: error: ld returned 1 exit status make[2]: *** [pFEM] Error 1 make[1]: *** [CMakeFiles/pFEM.dir/all] Error 2 make: *** [all] Error 2 Any comments? thanks in advance, Juan Araújo El viernes, 1 de diciembre de 2017,

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

2017-12-01 Thread Juan Carlos Araujo Cabarcas
Any ideas? Thanks in advance. El viernes, 1 de diciembre de 2017, 8:09:26 (UTC-5), Wolfgang Bangerth escribió: > > On 11/30/2017 03:05 PM, Juan Carlos Araujo Cabarcas wrote: > > > > When using the PETSc wrappers as in step-36 (installed with PetscScalar > = > > comp

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

2017-12-01 Thread Juan Carlos Araujo Cabarcas
Any ideas? Thanks in advance. El viernes, 1 de diciembre de 2017, 8:09:26 (UTC-5), Wolfgang Bangerth escribió: > > On 11/30/2017 03:05 PM, Juan Carlos Araujo Cabarcas wrote: > > > > When using the PETSc wrappers as in step-36 (installed with PetscScalar > = > > comp

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

2017-11-30 Thread Juan Carlos Araujo Cabarcas
std::complex’) bval); ^ Errors: in attached file. Any ideas? Thanks in advance. El domingo, 12 de noviembre de 2017, 18:32:00 (UTC-5), Wolfgang Bangerth escribió: > > On 11/10/2017 08:15 AM, Juan Carl

Re: [deal.II] Using DataOut with MappingCollection

2017-11-28 Thread Juan Carlos Araujo Cabarcas
Thanks for your quick reply and all the support on this thread. Programming with templates it's not really my cup of tea! El martes, 28 de noviembre de 2017, 0:10:15 (UTC-5), Wolfgang Bangerth escribió: > > On 11/27/2017 08:13 PM, Juan Carlos Araujo Cabarcas wrote: > > > >

Re: [deal.II] Using DataOut with MappingCollection

2017-11-27 Thread Juan Carlos Araujo Cabarcas
t;2, 2> >::CurvedCellRegion’ make[2]: *** [CMakeFiles/pFEM.dir/pFEM.cc.o] Error 1 make[1]: *** [CMakeFiles/pFEM.dir/all] Error 2 make: *** [all] Error 2 Any ideas? El lunes, 27 de noviembre de 2017, 12:45:04 (UTC-5), Wolfgang Bangerth escribió: > > On 11/17/2017 12:37 PM, Juan

[deal.II] Using DataOut with MappingCollection

2017-11-17 Thread Juan Carlos Araujo Cabarcas
Dear all, As described in the forum: https://groups.google.com/forum/?utm_source=digest_medium=email#!topic/dealii/ijJFbOTZxkk I would like to reproduce step-27 but with curved boundaries with the use of MappingCollection. Everything seems to work fine, but I noticed that data_out does not

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

2017-10-23 Thread Juan Carlos Araujo Cabarcas
Thanks for the quick answer! My problem is the scalar wave equation and the boundary conditions are simple homogeneous Dirichlet. Before starting with hp I was successfully using ZeroFunction(). How would this "*map taking a types::boundary_id and a Function object*" look like? Any step in the

[deal.II] mapping_collection and step-27

2017-10-23 Thread Juan Carlos Araujo Cabarcas
Dear all, I would like to reproduce step-27 for a curved domain instead of the given domain. For this, I would like to make use of a MappingCollection. I am having problems using interpolate_boundary_values, suggested in:

Re: [deal.II] lexicographic ordering in 3D

2017-07-11 Thread Juan Carlos Araujo Cabarcas
Thanks for your quick and useful answer, and explanation ... GridReordering was what I needed! Now I see that I did not mean lexicographic order, that I had assumed from the begining in 2), but it didn't read like that in the email. I apologize for the wrong formulation of the question. *Some

Re: [deal.II] Re: How to initialize PETScWrappers::SparseMatrix from a PETSc Mat

2017-06-06 Thread Juan Carlos Araujo Cabarcas
Thanks for answering! The warnings come from my code, which is in developing stage! In the installation I had before, [-Wunused-variable],[-Wunused-parameter],[-Wunused-but-set-variable] did not show up. It's a bummer that changing deal.II version will make such a difference! As I am not

[deal.II] Re: How to initialize PETScWrappers::SparseMatrix from a PETSc Mat

2017-06-05 Thread Juan Carlos Araujo Cabarcas
Well, now I tried to update my code ... I was using: git clone https://github.com/davydden/dealii cd dealii git checkout branch_petscscalar_complex Now, because I would like to use a few recent additions from the development version, I do instead: git clone https://github.com/dealii/dealii 1)

[deal.II] Re: How to initialize PETScWrappers::SparseMatrix from a PETSc Mat

2017-05-24 Thread Juan Carlos Araujo Cabarcas
** \n\n" ); } else { printf("\n\n Mass matrix successfully read! \n\n" ); } Any simpler alternative? By the way, I noticed also the definition: MatrixBase::operator Mat () const { return matrix; } By I was unable to use it! Thanks, Juan Carlos Araújo, Umeå Uni

[deal.II] How to initialize PETScWrappers::SparseMatrix from a PETSc Mat

2017-05-24 Thread Juan Carlos Araujo Cabarcas
Dear all, I have a matrix written by using PETSc in a Binary file, and I am trying work with it in my dealii environment. For this, I read a system matrix from PETSc: Mat M_read; PetscViewerBinaryOpen(PETSC_COMM_WORLD,file.c_str (),FILE_MODE_READ,); MatCreate(PETSC_COMM_WORLD,

[deal.II] state of dealii in complex arithmetics

2017-05-23 Thread Juan Carlos Araujo Cabarcas
Dear all, There are features in SLEPc 3.7 that I am interested in, and will try to make a fresh re-installation soon. I have been using deal.II 8.3 with complex arithmetic’s from the branch: git checkout branch_petscscalar_complex and I wonder about the state of that branch, or if it has

Re: [deal.II] Re: Arbitrary mapping or Manifold?

2017-04-07 Thread Juan Carlos Araujo Cabarcas
I will see places where it may fit, I will come back to you. Hej, Toby! nice to hear that it may help, looking forward to see the result =) El viernes, 7 de abril de 2017, 15:42:13 (UTC+2), Wolfgang Bangerth escribió: > > On 04/07/2017 07:24 AM, Juan Carlos Araujo Cabarcas

[deal.II] Re: Arbitrary mapping or Manifold?

2017-04-07 Thread Juan Carlos Araujo Cabarcas
In the same spirit I would like to share a code that uses ChartManifold, for generating a mapping of an ellipse applied to HyperBall. El miércoles, 5 de abril de 2017, 14:47:17 (UTC+2), Juan Carlos Araujo Cabarcas escribió: > > Dear all, > I recently found the thread: Someth

[deal.II] Re: Arbitrary mapping or Manifold?

2017-04-05 Thread Juan Carlos Araujo Cabarcas
Araújo, Umeå Universitet El martes, 9 de febrero de 2016, 17:10:51 (UTC+1), Juan Carlos Araujo Cabarcas escribió: > > Dear all, > > I am working with the wave equation with variable refractive indexes > within the domain. > I receive meshes that come from a shape optimization ro

[deal.II] MappingQ, and minimal test for reproducing Benchmark

2017-04-04 Thread Juan Carlos Araujo Cabarcas
Dear all, I'm am trying to reproduce with my implementation, the results in the Photonic Crystal computations performed in [1]. Here, the author uses a grid with an inner disk with radius R=0.475, and for FEM it is used the software Concepts [2] that implements curvilinear elements denoted

[deal.II] Re: periodic boundary conditions

2016-09-23 Thread Juan Carlos Araujo Cabarcas
Hej, I found the issue ... I forgot the line constraints.distribute ( sol ); after the solution routine ... so it was just mainly a visualization issue. Thanks for the quick response! El jueves, 22 de septiembre de 2016, 19:26:28 (UTC+2), Daniel Arndt escribió: > > Juan, > > [...] >> 1)

[deal.II] Fast FE evaluation over several points and several functions

2016-06-22 Thread Juan Carlos Araujo Cabarcas
Dear all, I am solving an eigenvalue problem similar than step-36. After solving for the eigenpairs, I evaluate the eigenfunctions in the standard way: VectorTools::point_value ( mapping, dof_handler, efun[m], q_points[j], Uq ); where: efun[m] is the m-th eigenfunction from step-36,

[deal.II] precompiling independent classes as a library and link it to the main project

2016-05-30 Thread Juan Carlos Araujo Cabarcas
Dear all, I have several independent classes that use deal.II objects, and I would like to pre-compile them as shared libraries to be linked to the main executable. The idea is to be able to do something similar than what is described here: