Dear Praveen,
there is a bug in SphericalManifold for dimension 3, which makes it unusable so
far with MappingManifold.
The issue is with the computation of the tangent vector across the north and
south poles.
At the moment MappingManifold is the only class using this computation, so if
you
he firstChartManifold to the embedding
space of the second ChartManifold"
In your case, F is your spherical manifold, while your G is the rotation
manifold (to preserve periodicity of the spherical manifold).
Best,
Luca.
>
> Thanks
> praveen
>
> On Mon, May
he other faces.
>
> When I run my code with this, I find large errors arising in the region where
> north/south faces meet the other faces.
>
> Is there some incompatibility in my mapping ? Will the face normal vectors
> (which lie in tangent plane) be correct at
e not normal to sphere surface, but normal to edges on the sphere,
> which should be tangent to sphere surface.
>
> I get these normals from
>
> FEFaceValues<2,3> and using update_normal_vectors
>
> Best
> praveen
>
> On Tue, May 17, 2016 at 12:19 PM, luca.hel
rrors at the junction
> region of the two charts.
>
> Best
> praveen
>
> On Tue, May 17, 2016 at 12:30 PM, luca.heltai wrote:
> Then can you verify which ones are correct, by taking the scalar product with
> the quadrature points (assuming a sphere centered at the origin), and c
; Best
> praveen
>
> On Tue, May 17, 2016 at 1:55 PM, luca.heltai wrote:
> Dear Praveen,
>
> yes, I was expecting those numbers. Mapping Manifold should be machine
> precision. I’d expect Q9 to be a little better though… I’ll try to look into
> this (when I have som
Sam,
I’d suggest you to use a shared pointer as member variable:
std::shared_ptr > f_ptr;
Then the first time you know how to construct the function, you do this as
f_ptr = std::shared_ptr > (new FEFieldFunction<…>(here your
construction objects);
you access this as a standard pointer:
f
As long as the solution vector is the same (say current_solution), this will
work. Should you need to get another solution at next time step, just repeat
the
> f_ptr = std::shared_ptr > (new FEFieldFunction<…>(here
> your construction objects);
phase, and the previous one will be replaced wi
Sorry for replying late…
You could try this package:
https://github.com/luca-heltai/dealii/releases/tag/v8.5.0pre.v2
It is a deal.II library on steroids (650MB), with a full brew installation
inside. I installed the brew tree under
/Applications/deal.II.brew.app/Contents/Resources/brew/
which
David,
if you call
triangulation.set_manifold(id, new_manifold)
the old manifold is automatically discarded, unless there is more than one
manifold id using it.
Maybe you are setting the manifold inside a function, and you are exiting the
function (and destroying the manifold) before the tri
What J-P is suggesting is the following:
Have a
std_cxx11::shared_pointer solver
member in your class, then
if (!solver) {
solver = std_cxx11::shared_pointer(new TrilinosWrappers::SolverDirect direct(…) );
solver->initialize (mass_matrix);
}
this will be done only once, and
Dear Tom,
you’re not doing anything wrong. The Manifold/Boundary interfaces are
undergoing a lot of restructuring, and you hit a problem of backward
compatibility implementation…
Your issue is in the following lines of code, that were inserted to maintain
backward compatibility between Manifol
There is nothing wrong with what you are doing.
The problem is in the nature of your Manifold, since it contains 3 singular
points. The first is the center of the ellipsoid, while the second and third
are the north and south poles.
Try attaching a SphericalManifold to your deformed grid, and s
Dear Praven,
the issue seems to be the same as before, but we have not been able to pin
point what the problem could be.
In particular, the problem is not in Manifold *per se*. It may be related to
the way MappingQ and MappingQGeneric compute support points. I have verified
that all quadratur
n. Give me some places
> where I should start looking. I suppose I should start with MappingQ and
> MappingQGeneric.
>
> My shallow water project is stuck due to this, so I am really keen to resolve
> this issue.
>
> Best
> praveen
>
> On Wed, Aug 10, 2016 at 1:32 PM
For future reference, if you are in doubt on how to use a single
class/function, you can skim through the directory tests, i.e.,
tests/manifold/composition_manifold_*
would have shown to you an example which is identical to what you are trying to
achieve, and in general it contains very simple
Yes. This is a problem inherent with the fact that it is not possible to
approximate a sphere with a single chart without singularities.
On the other hand, the manifold you use to describe the top and bottom part of
the sphere are not continuous, so the triangulation gets confused when you add
Not necessarily. In fact the ordering is reshuffled so that vertex support
points comes first, then edge support points (3D), face support points, and
lastly cell support points.
To see the actual ordering, take a look at
fe.get_unit_support_points()
https://www.dealii.org/8.4.1/doxygen/deal.
> Now I have to figure out some other way to do this efficiently. I want to
> avoid use of get_function_values.
>
> Best
> praveen
>
> On Wed, Nov 23, 2016 at 3:27 PM, luca.heltai wrote:
> Not necessarily. In fact the ordering is reshuffled so that vertex support
>
Ciao Dirk,
unfortunately not. I’m forwarding your mail to the mailing list. Maybe someone
there has a better idea…
L.
> On 29 Nov 2016, at 19:19, Dirk Peschka wrote:
>
>
> Hey Luca,
>
> duing the petsc compilation, when making mumps, I get some strange
> errors and the installation fails.
On 2 Dec 2016, at 15:01, Bruno Turcksin wrote:
>
>
> 2 - Are there facilities of some kind that can help in generating the
> assembly code? In Fenics I just specified the weak formulations.
> You need to write the code in C++ and loop over cells, quadrature points, and
> basis functions yourse
Dear Anup,
one option you have is to use a Triangulation.
I would still use two scalar extractors:
const FEValuesExtractors::Scalar u_x
const FEValuesExtractors::Scalar u_y
Tensor<2,dim+1> gradient;
gradient[0] = fe_values_ref[u_x].gradient(k, q_point)
gradient[1] = fe_values_ref[u_y].gradien
Dear Joaquin,
as I wrote to you on one of my previous emails, the repository
https://github.com/luca-heltai/ans-ifem
contains a branch that compiles correctly, without those warnings, on the
development version of deal.II.
There are some more warnings due to the deprecation of the ‘sadd’ func
Actually, you have some alternatives…
Take a look at this class:
https://github.com/mathLab/deal2lkit/blob/master/include/deal2lkit/parsed_mapped_functions.h
It does something similar to what you are asking for.
(documentation here:
http://mathlab.github.io/deal2lkit/class_parsed_mapped_func
Dear Franco,
I think there is a problem in your formulation…
You are integrating
u_i n_i f_i
and while you can certainly come up with good reasons for this to work out, I’m
unsure about what you want to achieve. If your boundary conditions are
((grad(u)).n,v) = u_i,j n_j v_i
and you impos
> On 17 Feb 2017, at 14:46, Franco Milicchio wrote:
>
> I agree I have a translation problem here. In practice, I am integrating a
> weak problem as this (in Fenics' UFL terms, but easy to recognize):
>
>
> inner( sigma(u), eps(v) ) * dx = inner(f, v) * dx + inner(f1, v) * ds(1)
I’m gue
Hi Alex,
I’d use a vector of
std::vector > > v;
which are light objects, and can be resized and reshaped. Whenever you do a
push_back, you’d have to use
v.push_back(std_cxx11::shared_pointer(new
ParsedFunction(…) );
then
v[i]->value(…)
would work.
> On 22 Mar 2017, at 16:31, Alex Zimm
Resizing is fine. It is the equality that won’t work:
function_pointers[f] =
std::shared_ptr>(new
dealii::Functions::ParsedFunction()));
is the correct way to do it.
L.
> On 23 Mar 2017, at 15:09, Alex Zimmerman wrote:
>
> Also for the record, here's a copy of issue_10_push_back.cc from
Or, you could use as a stable pair of finite elements the FE_Q(k)-FE_DGP(k-1)
pair:
D. Boffi and L. Gastaldi. On the quadrilateral Q2-P1 element for the Stokes
problem. Int. J. Numer. Meth. Fluids, 39 (2002), 1001-1011.
This pair of finite elements is known to produce consistently better result
Same here. I can download the package both from home and from my office.
L.
> On 8 Apr 2017, at 18:35, Daniel Arndt
> wrote:
>
> Ana,
>
> I can't confirm your problem. Downloading works without issues. What browser
> are you using? Do you get any warnings?
>
> Best,
> Daniel
>
> Am Samsta
This is actually not a bug. A spherical manifold is undefined on the center of
the sphere/ball. You cannot attach set_all_manifold_ids of the spherical
manifold and attach a spherical manifold to it, since the central cell (the one
where the cell->center() coincides with the ball center) would n
What image are you using?
Try with
dealii/dealii:v8.5.0-gcc-mpi-fulldepscandi-debugrelease
my guess is that you are using
dealii/dealii:v8.5.0-gcc-mpi-fulldepsmanual-debugrelease
a bug with that image was reported by Timo, and I thought I had removed it from
the image list…
L.
> On 11 May
Did you install xcode and xcode-tools?
It looks like deal.II cannot find a lot of things…
The application opens a terminal and tries to figure out if you did. Does the
terminal open up with no issues?
L.
> On 5 Jun 2017, at 8:00, Sudarshan Kumar wrote:
>
> I have installed deal.II-8.5-br
No need to restore to the old version. I have 10.12.5 myself.
Can you try running xcode, and see if it asks you wether you want to upgrade
the command line tools?
L.
> On 5 Jun 2017, at 8:22, Sudarshan Kumar wrote:
>
>
>
> On Monday, June 5, 2017 at 11:39:28 AM UTC+5:30, Luca Heltai wrote:
If you don’t have XCode, then you don’t have XCode tools. You can download
XCode from the App store.
This error:
> The CMAKE_C_COMPILER:
>
>
>
>
> /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/cc
>
>
>
> is not a full path to an existing comp
Felix,
also, if the point you want is not a quadrature point, you can exploit the
following trick:
Point my_point_in_reference_cell;
std::vector > quad(1, my_point_in_reference_cell);
and then you can initialize an
FEValues fev(fe, quad, update_inverse_jacobians);
This will give you access
There may be a problem with the way boundary ids are set.
Can you try the following?
after creating the grid with removed cells, loop over all cells and all faces,
and if at boundary with boundary id == 1, then call
cell->face(f)->set_all_manifold_ids(1);
notice the “_all_”, i.e., it i
Dear Simon,
one solution you have is to pass a *Triangulation*, and then generate the dof
cell accessors using the constructor
typename DoFHandler::cell_iterator cell1(*tria_iterator, dh1);
typename DoFHandler::cell_iterator cell2(*tria_iterator, dh2);
Best,
Luca.
> On 30 Apr 2021, at 18:37
In addition to that, I’d suggest you use the
GridTools::Cache::get_vertex_to_cell_map class, which builds the map once, and
stores it for future reference. This one uses GridTools::vertex_to_cell_map,
which builds the map with one pass, instead of looping over all cells.
Indeed, if you want to
Alternatively, if you want to use deal.II in your daily terminal, you could
follow the instructions that come in the deal.II terminal, i.e., add this to
your .zshrc:
export DEAL_II_CONF_SILENT=ON
. /Applications/deal.II.app/Contents/MacOS/dealii.conf
Best,
Luca.
> On 28 Jun 2021, at 21:52, Je
If you want to use symbolic calculations, you could also leverage SymEngine,
and use Functions::SymbolicFunction
https://www.dealii.org/current/doxygen/deal.II/classFunctions_1_1SymbolicFunction.html
In your code, you could create such object by making a
std::unique_ptr> fun;
parse its expres
Dear Lucas,
you could use GridTools::remove_anisotropy
Best,
Luca.
> On 18 Sep 2017, at 16:29, Lucas Campos wrote:
>
> Dear Bruno,
>
> You will find attached the resulting grids. While I originally expected to
> find new divisions only along the x-direction, your previous answer tells me
Yes. It is trying to use the accelerate framework, but since it was installed
using 10.12, it links to that one.
If you install the 10.12 SDK you should have what you need.
L.
> On 11 Oct 2017, at 8:58, Alberto Salvadori wrote:
>
> Sure, I will. Am I understanding properly that deal.ii is tr
>> — Are there any applications with immersed FEM (cut-cell, etc.) with deal.II?
Yes. There are at least two open source codes, published on “Archive of
Numerical Software”:
http://journals.ub.uni-heidelberg.de/index.php/ans/issue/view/2930
"Carraro, Wetterauer: On the implementation of the eX
Oh, and naturally ’s/brew/spack/g’ in my previous email.
clang --version
Apple LLVM version 9.0.0 (clang-900.0.38)
Target: x86_64-apple-darwin17.0.0
Thread model: posix
InstalledDir:
/Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin
and
gfortran --version
It does not compile, anyway…
:D
Maybe Denis can help: I get the following when trying to compile pkg-config
using gcc installed with spack:
In file included from
/System/Library/Frameworks/Security.framework/Headers/AuthSession.h:32:0,
from
/System/Library/Frameworks/Securit
Dear Alberto,
> GridGenerator::hyper_ball ( triangulation, center, cell_radius );
>
> static SphericalManifold volume_boundary;
>
>
> triangulation.refine_global ( 3 );
>
> triangulation.set_all_manifold_ids_on_boundary(0);
>
> triangulation.set_manifold (0, volume_boundary);
Y
Dear Yaakov,
did you even try google?
:)
google search: dealii isogeometric analysis
https://www.google.it/search?q=dealii+isogeometric+analyis&gws_rd=cr&dcr=0&ei=vTkVWvurIcOUaPbVr6AO
gives you as a first result this repository:
https://github.com/mathLab/IGA-dealii
So far it supports singl
Dear all,
registration is now open for the Sixth deal.II Users and Developers Workshop,
to be held in Trieste@SISSA, from July 23 to July 28.
http://indico.sissa.it/e/sixth-dealii-workshop
The program includes invited talks and discussions by users and developers of
the deal.II library in are
M…
I’m actually not sure it’s a bug… The point where the error starts growing is
at 10e-12. This is basically machine precision, considering the fact that you
are computing the L2 error.
After you have reached machine precision, any addition you make to that, is
just roundoff error addin
If all you are solving is a two dimensional problem, you could encode your
“get_function_values” into a matrix vector multiplication to drastically
improve the situation.
I’m thinking of a matrix of size (n_quadrature_points x n_active_cells) x
n_dofs, and then you slice the results cellwise
multiplication (maybe with one identical Vandermonde matrix).
L.
> On 31 Dec 2017, at 6:09, Praveen C wrote:
>
>
>
>> On 30-Dec-2017, at 11:40 PM, luca.heltai wrote:
>>
>> I’m thinking of a matrix of size (n_quadrature_points x n_active_cells) x
>> n_d
Take a look at this:
https://github.com/luca-heltai/bare-dealii-app
it is a bare structure, for medium to large size problems, with separate
source/include/tests directories, that support travis CI integration.
L.
> On 11 Jan 2018, at 15:52, 'Maxi Miller' via deal.II User Group
> wrote:
>
You could use KDTree, based on nanoflann.
If you installed the development version of dealii with NANOFLANN support, then
you’d have access to algorithms that work in order log(N) complexity for each
search, and N log(N) to build the tree.
You’d have to massage it a little, as at the moment i
After you have created your coarse mesh, you could
GridTools::flatten_triangulation it. This will create a brand new coarse
triangulation, containing only all your active cells.
L.
> On 19 Feb 2018, at 13:35, Lucas Campos wrote:
>
> Dear all,
>
> I am creating a mesh using GridGenerator::hy
Kurt,
in principle deal.II does not require OCE. It can work also with OCCT (it will
use whatever it finds available, unless someone changed this behaviour).
When OCCT will be packaged into debian, I’m sure we can make dealii use that
instead of OCE.
Best,
Luca.
> On 11 Mar 2018, at 4:29, Ku
Dear all,
the following announcement may be of interest to you or to some of your
students. We are actively seeking for a deal.II experienced programmer ASAP!
Best,
Luca.
--
A one year postdoctoral position is available at SISSA @ mathLab, starting from
May 1, 2018
(see https://goo.gl/247WN
The problem is that you cannot read a mesh with hanging nodes. If you have a
mesh with hanging nodes, all faces with a hanging node are actually boundary
faces, so in your picture, all faces that are between a refined and an
unrefined cells, are actually boundary faces, i.e., they are topologica
Dear all,
I have just uploaded a new package for deal.II-9.0.0-rc2 here:
https://github.com/dealii/dealii/releases/edit/v9.0.0-rc2
This was compiled with Apple clang 9.1.0, and gfortran from gcc 7.3.
Please test, and let me know if everything is ok.
Best,
Luca.
--
The deal.II project is loc
ay, May 6, 2018 at 11:51:19 AM UTC+2, Luca Heltai wrote:
> Sorry. The address is the following:
>
> https://github.com/dealii/dealii/releases/v9.0.0-rc2
>
> Luca
>
> Il giorno 06 mag 2018, alle ore 11:30, luca.heltai ha
> scritto:
>
>> Dear all,
>>
>
The deal.II source directory can be found under
/Applications/deal.II-9.0.0.app/Contents/Resources/spack/src/deal.II-9.0.0
There you will find the examples directory.
In the new package (that I’m building right now), you’ll also have in the usual
location.
Best,
Luca.
> On 12 May 2018, at 1
Dear All,
the deadline for the submission of talks proposal for the “Sixth deal.II Users
and Developers Workshop” is on the 23 of June (one month from now).
You are invited to register and submit an abstract through the online procedure:
http://indico.sissa.it/e/sixth-dealii-workshop
If you n
Dear all,
thanks for the feedback. Here you’ll find a new version of the deal.II package,
made with spack, and compiled with clang@6.0.0, and gfortran@8.1:
https://github.com/luca-heltai/dealii/releases/tag/v9.1.0-pre-002bd21ce6
The package is independent on apple compilers, should work on mos
Dear Riku,
this is not (yet) supported by deal.II. You will have to attach manually a
CylindricalManifold after extrusion. The library does not know yet how to
compute the extrusion of a general manifold object, therefore it removes all
manifold objects from the extruded triangulation (as they
Dear Marek,
> On 29 Sep 2018, at 22:58, Marek Čapek wrote:
>
> Hi,
> I have a geometry in step file format. I would like to import it and generate
> from it a mesh digestible by dealii GridIn functions.
> I followed step-54, modifying thefunction read_domain
>
> https://www.dealii.org/8.4.1/do
Dear Rochan,
you can find prebuilt static binaries of deal.II 8.5.0 for windows, with no
external libraries, here:
https://github.com/luca-heltai/dealii/releases/tag/9.5.pre-windows
I also know of people that successfully compiled and used deal.II 9.0.0 on
windows.
L.
> On 12 Oct 2018, at 4
Dear Earl,
setting the Manifold should be done *before* you do refinment.
The line you added at the end, tria.execute_coarsening_and_refinement(), is not
doing any refinement, so you don’t see the effect of your set_manifold_id()
call.
What happens if you refine globally once more after the
Dear Praveen,
this is actually a feature, and not a bug. The reason is that, by default,
deal.II generates triangulations with boundary id = 0. If you don’t specify a
boundary id in the gmsh file, then all boundary faces get boundary id = 0.
This allows us to save space in the output grid, sin
Did you use the option “colorize=on”?
L.
> On 18 Nov 2018, at 10:45, Praveen C wrote:
>
>
>
>> On 18-Nov-2018, at 2:59 PM, luca.heltai wrote:
>>
>> Dear Praveen,
>>
>> this is actually a feature, and not a bug. The reason is that, by default,
I believe the first point is a bug of the read_ucd function.
I opened an issue for this (https://github.com/dealii/dealii/issues/7802)
The second point, however, is a feature. What should you do when two material
ids meet? I think it is the responsability of the user to make sure that the
corr
You may also find useful this link:
https://journals.ub.uni-heidelberg.de/index.php/ans/issue/view/2930
where the authors implement XFEM in deal.II to model interface problems.
L.
> On 18 Mar 2019, at 9:39, Daniel Arndt
> wrote:
>
> James,
>
> I am a graduate student who will be defending
Take a look at this PR for a few examples of usage of mesh_loop:
https://github.com/dealii/dealii/pull/7806
The function WorkStream::run
takes (a minimum of) 5 arguments:
WorkStream::run(cell, endc, cell_worker, copier, scratch, copy);
initial and final iterator, a worker function, a copier f
> Is there a reason then that there are several examples using
> integration_loop(), but (afaik) only one using mesh_loop?
Yes. mesh_loop is more recent w.r.t. integration_loop.
mesh_loop was introduced to address some of the oddities that are in
integration_loop, that make its use somewhat les
Nicola,
what would you think of a function in GridTools like
GridTools::propagate_internal_manifold_ids(tria, disambiguation_function)
that would loop over all cells, loop over all faces, check if the neighbor
manifold id is the same of this cell, and
i) the two ids are the same: assign the
You are integrating using two quadrature points per direction. Can you raise
that to (2*fe.degree+1)?
L.
> On 11 Apr 2019, at 10:58, bobspar...@gmail.com wrote:
>
> Dear all,
>
> First of all, thanks for the awesome library. I am starting to learn deal.II,
> and it is great! The documentation
If you plan to use any domain that is not a square (or an affine
transformation), you have to make sure you integrate exactly the product of two
polynomials of order degree and of the determinant of the Jacobian. This last
term is constant only for simple meshes, but it is the square root of a
Wolfgang, is that true also for mass matrices? I’d agree with you for stiffness
matrices, but I’d surprised this worked ok for mass matrices as well.
If so, I’ve always been over integrating in my life…
:)
L.
> On 12 Apr 2019, at 21:15, Wolfgang Bangerth wrote:
>
> On 4/12/19 8:41 AM, Robe
> On 16 May 2019, at 9:16, 'Maxi Miller' via deal.II User Group
> wrote:
>
> I would like to use the interpolate()-function to interpolate the correct
> function onto my mesh for comparison during some simple tests. This works
> fine when using FE_Q-elements, but fails when using FE_Bernstei
Did you look at this class?
https://www.dealii.org/developer/doxygen/deal.II/classEllipticalManifold.html
It is in deal.II 9.1.0…
L.
> On 3 Jun 2019, at 12:54, SebG wrote:
>
> Hi Thomas, hi Luca,
>
> I would like to re-generate the ellipsoidal with deal.ii 9.0.1 but the
> interfaces of the
Dear Sebastian,
have you tried using a simple SphericalManifold for your problem?
The way spherical manifold works is by interpolating *both* the radius and the
angles, so that points intermediate between two points with different radii get
an average radius between the two.
This may not be
Jhon,
the path is there, and it is correct. The message is telling you that you are
in a directory containing developement version of the examples.
If you inspect the CMakeLists.txt of the example you are trying to compile, it
tells you:
FIND_PACKAGE(deal.II 9.2.0 QUIET
HINTS ${deal.II_DIR
Dear all,
I’d like to advertise an opening for a Postdoctoral position (18 months) in my
group, to work with and on deal.II. The deadline is pretty tight… please send
me a private message if you are interested.
Luca.
---
Title of the position: Adaptive methods for solving multi-physics
prob
Dear All,
I would like to bring at your attention that the deadline to apply for a Phd
position in Mathematical Analysis, Modelling and Applications at SISSA, Scuola
Internazionale Superiore di Studi Avanzati Trieste, Italy is August 22nd, 2019
at noon (Rome time).
Applied mathematics activiti
Alberto,
be aware of the fact that, when you extract the boundary mesh, the orientation
of the cells may be different w.r.t. to the corresponding face on the
volumetric mesh. In other words, using
FEValues codim_one(…);
FEFaceValuescodim_zero(…);
codim_one.reinit(codim_one_cell);
co
Dear Felix,
by any chance, did you take a look at step 10?
https://www.dealii.org/current/doxygen/deal.II/step_10.html
This step explains a little bit what to do when you want to solve on curved
domains with high order finite elements. In particular, you need to ensure that
the mapping you ar
Dear Ester,
please take a look at the release page:
https://github.com/dealii/dealii/releases/tag/v9.1.1
There is a package for Mac OS Catalina. In Catalina, permissions are much more
restrictive, and you’ll have to explicitly enable the application to run (the
first time), by right clicking o
You can do
export PATH=$(spack location -i mpich)/bin:$PATH
to prepend the path of mpich to your PATH. This is not done automatically.
Alternatively, you could also run, as root:
cd /usr/local
spack view add . mpich
which aliases all mpi commands to /usr/local/bin
I usually do this in a cus
You don’t need to do anything special. Just create a MappingManifold, and pass
it where required. The manifold that is used is the one associated to the
triangulation. You don’t need to associate it to the Mapping, since the mapping
is used on *all* objects of the triangulation, and it will know
> CylindricalManifold<2> cylman(0);
> triangulation.set_manifold(numbers::flat_manifold_id, cylman);
The flat_manifold_id is usually reserved for the flat manifold, and yours is
not a flat manifold. Setting the flat manifold to the cylindrical manifold
implies
You should set the manifold id
ow recompile with: $ make
> ***
> [100%] Built target release
> bash-3.2$ make
> Scanning dependencies of target step-20
> [ 50%] Building CXX object CMakeFiles/step-20.dir/step-20.cc.o
> [100%] Linking CXX executable step-20
> [100%] Built target step-20
> bash-3.2$ ./step-20
tails, but I remember having some issues with XCode
> during the installation. I thought it had automatically updated correctly,
> but on closer inspection there was an error message (nothing an online search
> couldn't solve).
>
> Ester
>
>
> El dimarts, 3 desembre de 2019 11:
ectly,
> but on closer inspection there was an error message (nothing an online search
> couldn't solve).
>
> Ester
>
>
> El dimarts, 3 desembre de 2019 11:08:57 UTC-5, luca.heltai va escriure:
> I have not updated to the latest version yet. They changed clang, xcod
tel 030 3711239
> fax 030 3711312
>
> e-mail:
> alberto.salvad...@unibs.it
> web-page:
> http://m4lab.unibs.it/faculty.html
>
>
>
> On Fri, Dec 6, 2019 at 11:18 AM luca.heltai wrote:
> Alberto,
>
> this may be related to different hardware and diffe
Dear Juan,
The logic you should use is: when your neighbour is finer, you should not
assemble the full face term from the coarse neighbour. Instead, you have to
assemble your DG terms using a FESubfaceValues initialized with the correct
arguments on the coarse face, to match *each* of the fine
Dear Ellen,
you may want to compare with this:
https://github.com/luca-heltai/dealii/pull/91/files#diff-acfb3c7b43e4935be016fda6ebdc5881
It is a parallel version of step-38, that never got into the library, written
by one of our students, during a deal.II workshop in Trieste (held by Timo and
Nicola,
I think you are hitting a periodicity issue here.
The points with angle 7/4 pi and 5/4 pi are at a distance (from pi/4) of 6/4 pi
= -pi/2 (this is the shorter distance, and this is what is used in the code),
and pi = -pi (this could be either -pi or pi, and it is in fact not well
def
Nicola,
I would simply do this by
std::mapline(0)), unsigned int> line_count;
for(const auto &cell: tria.active_cell_iterators()) {
for(unsigned int i=0; i< GeometryInfo::lines_per_cell; ++i) {
const it = line_count.find(cell->line(i));
if(it != line_cou
Did you try the pre-compiled deal.II package?
Both
dealii-9.1.1-clang-9.0.0.dmg
and
dealii-9.1.1-catalina-haswell-ro.dmg.zip
from https://github.com/dealii/dealii/releases/tag/v9.1.1
should work.
Luca.
> On 25 Jan 2020, at 21:29, Ihar Suvorau wrote:
>
>
> Software:
> • macOS Ca
By the way, the examples in the docker image are located under
$(spack location -i dealii)/share/deal.II/examples
i.e., you can go to
cd $(spack location -i dealii)/share/deal.II/examples
and from there, compile and run. Notice that they are not pre-compiled (they
would just take up download
Bruno, it seems like you are attaching your faces to the *boundary* manifold.
Let me try to explain what is happening.
In the code for step-54, there is a wire that is used to describe the manifold
of the *boundary* of the surface (a curve of dimension one embedded in
dimension three). This is
1 - 100 of 138 matches
Mail list logo
