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
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?
&
> 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
>
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
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 p
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_lambda = [&
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 = [&polar](Vector a, Vectorb) {return
> polar.Callback(a,b);};
> opt.solve ( my_l
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,
>>
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 exam
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 ent
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
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
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
(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ó:
> >
> >
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 ‘MPI_Comm_create_group
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} \
; 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
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 is
> Daniel
>
> Am Samstag, 9. Dezember 2017 00:48:04 UTC+1 schrieb Juan Carlos Araujo
> Cabarcas:
>>
>> Thanks again,
>>
>> Now it passes the *interpolate_boundary_values *call.
Thanks again,
Now it passes the *interpolate_boundary_values *call. However, it is
conflicting now with:
MatrixTools::apply_boundary_values (bval,
system_matrix,
lect2: 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 Bangerth
escribió:
>
> On 12/05/2017 05:47 PM,
å Universitet
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_b
ollect2: 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, 15:49:48 (UTC-5), Wolfgang Bangerth
escribió:
>
> On 12/01/2017 01:29 PM, Juan Carlo
e 2017, 15:06:29 (UTC-5), Bruno Turcksin
escribió:
>
>
> 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 gav
ngerth
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
> =
> > complex) and using MappingCollection as discussed before. I now
> try:
> >
> > |
>
ngerth
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
> =
> > complex) and using MappingCollection as discussed before. I now
> try:
> >
> > |
>
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 Carlos Araujo Cabarcas wrote:
> > Yes, it explains the Function
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:
> >
>
: *** [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 Carlos Araujo Cabarcas wrote:
> >
> > I would like to reproduce step-27 but with curved boundaries with the
> > use o
Dear all,
As described in the forum:
https://groups.google.com/forum/?utm_source=digest&utm_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
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
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:
https://www.dealii.org/8.5.0/doxygen/deal.II/namespaceVectorTools.html#a0
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 co
Dear all,
I have been designing 2D grids according to the standards given in:
https://www.dealii.org/8.4.1/doxygen/deal.II/structGeometryInfo.html
The procedure goes as follows,
1) we introduce a set of vertices as points (x_i,y_i), numbered by the
index "i",
2) connect them forming quads/ce
n Debug mode but not in Release mode.
>
> praveen
>
> On Tue, Jun 6, 2017 at 8:07 AM, Bruno Turcksin > wrote:
>
>> Juan Carlos
>>
>> On Tuesday, June 6, 2017 at 5:34:09 AM UTC-4, Juan Carlos Araujo Cabarcas
>> wrote:
>>>
>>> The war
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 profi
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) Fi
Thanks a lot Martin,
I had been waiting for this type of solution in dealii for a long time!
I also read in another thread about using MappingQGeneric instead of
MappingQ, so I updated my test code as below, and everything worked like a
charm!
PolarManifold boundary;
functionality for
> PETScWrappers::MatrixBase.
>
> On Wednesday, May 24, 2017 at 5:18:46 AM UTC-4, Juan Carlos Araujo
> Cabarcas wrote:
>>
>> Dear all,
>>
>> I have a matrix written by using PETSc in a Binary file, and I am trying
>> work with it in my de
{
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å Universitet
El miércoles, 24
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,&fd);
MatCreate(PETSC_COMM_WORLD,
Thanks, Denis and Wolfgang for the quick reply.
I will give it a try and report back whatever issue I encounter.
El miércoles, 24 de mayo de 2017, 7:26:42 (UTC+2), Denis Davydov escribió:
>
>
> On Tuesday, May 23, 2017 at 5:22:44 PM UTC+2, Juan Carlos Araujo Cabarcas
> wrote:
&g
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 be
(UTC+2), Juan Carlos Araujo
Cabarcas escribió:
>
> 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 escrib
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 wrote:
>
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: Something w
arlos 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 optimizati
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 Blen
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) M
Thanks for the fast reply, this is exactly what I needed!
El miércoles, 22 de junio de 2016, 19:38:35 (UTC+2), bangerth escribió:
>
>
> > I am solving an eigenvalue problem similar than step-36. After solving
> for the
> > eigenpairs, I evaluate the eigenfunctions in the standard way:
> >
> >
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,
aster/cdr
>
> Let me know if this helps.
>
> Thanks,
> David Wells
>
> On Monday, May 30, 2016 at 10:35:14 AM UTC-4, Juan Carlos Araujo Cabarcas
> wrote:
>>
>> Dear all,
>>
>> I have several independent classes that use deal.II objects, and
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:
http://www.cprogramming.com/tutorial/shared-librar
53 matches
Mail list logo
