Hi,
When you call reinit on hp::FEValues you can send in which of the elements
in the collection you want to reinitialize with.
That is, when you call reinit you call it as
fe_values.reinit(cell, q_index, mapping_index, fe_index);
where the additional integers specify which element in the colle
Well, I guess it's a question of how one should interpret this sentence.
Anyway, if you want to avoid calling reinit with extra parameters during
assembling (might be simpler anyway) you can probably achieve the same
thing by adding the same element twice to FECollection.
For example:
fe_colle
Hi.
Regarding this:
> On Thursday, April 6, 2017 at 1:51:41 PM UTC+2, Denis Davydov wrote:
>
> The issue with integrating and quadrature would still remain.
I've been implementing an algorithm to generate bulk and surface quadrature
given a level set function as boundary description. It's main
Hi,
In the stacktrace it looks like the error does not come from create_mass_matrix
but from VectorTools::project. It complains that the passed Function has 1
component when it should have dim components.
If you are projecting your own function, you need to make sure you set the
number of com
Hi,
From the error message, it looks like you don't have C++ compiler installed.
If you are using ubuntu, try
sudo apt-get install g++
and rerun cmake.
Best,
Simon
On 08/06/2022 16:32, Abhishek Nath Thakur wrote:
Respected deal.ii users,
I'm trying to install deal.ii 9.2.0 version but keep
Hi,
Just to elaborate slightly on what Chen said: If your installation stops at 52%, it is NOT because
of "dependencies not satisfied". This is not an error. Can you post the error you get
from calling make in the terminal? E.g. copy the output you get from cmake and "make" to
a txt-file and a
Hi,
This type of integral can be computed with the tools presented in the Step-85
tutorial:
https://www.dealii.org/developer/doxygen/deal.II/step_85.html
which you can use if you install the master branch of deal.II (or wait for the
soon-released 9.4 version).
In particular, take a look at t
Hi,
Since several XFEM papers formulate this as a Heaviside enrichment around the
crack, you might find this class useful:
https://www.dealii.org/current/doxygen/deal.II/classFE__Enriched.html
Relating to Wolfgangs answer, note in particular the section "Using enriched and
non-enriched FEs to
Hi,
Did you call
triangulation.execute_coarsening_and_refinement();
after setting the refinement flag?
Best,
Simon
On 23/06/2022 17:22, LY XXXiao wrote:
Hi there,
Has anyone experienced this issue in mesh refinement:
I have already a basic mesh with 120 cells, and I only want to refine one
Hi,
Unfortunately no-one has implemented FEInterfaceValues for hp yet (I think it
would be great if someone did). The only workaround I know is to NOT use
FEInterfaceValues and instead:
1. Use 2 hp::FEFaceValues objects, one initialized with cell 0 and one
initialized with cell 1.
2. Split t
Hi,
I have done the mex-coupling with dealii before. As Bruno said, it is easiest
to mex with cmake. Example in attached zip. Hope it helps.
Best,
Simon Sticko
On 15/08/2022 15:05, Bruno Turcksin wrote:
Simon
Le lun. 15 août 2022 à 08:52, Simon Wiesheier mailto:simon.wieshe...@gmail.com
Hi,
I think it's just the documentation that can be slightly improved here.
The function get_active_child_cells has a similar signature but it has
this note in its documentation:
* @note Since in C++ the MeshType template argument can not be deduced from
* a function call, you will have to speci
Hi,
Fixed in #9108, already merged in. I should probably have written that here for
completeness. Sorry.
Best,
Simon
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message bec
Hi,
Just to elaborate on what Daniel already said, Oliver, note that it is the
constructor of the Step3-class that calls the constructor of FE_Q,
that is
Step3::Step3()
: fe(1)
, dof_handler(triangulation)
{}
If you miss this call in the Step3-constructor, you will get the error you
stated
Hi,
a better option than computing the inverse is to factorize the matrix. This
can be done using the SparseDirectUMFPACK solver:
https://www.dealii.org/current/doxygen/deal.II/classSparseDirectUMFPACK.html
You might want to take a look at step 22, which uses this solver:
https://www.dealii.org
Hi,
I suspect that the condition you have would be easier to enforce weakly by
adding terms to the weak formulation, as is done in e.g. discontinuous
Galerkin.
Nevertheless, if you still want the normals at the nodes you can create a
"dummy" quadrature which has the nodes on the face as quadrat
Hi,
from the error message you can see that the element you are using only has
1 component. You get an error because you are trying to access component 1,
which doesn't exist. Since your element should have dim components, there
is likely something wrong with how you create your element. It shou
{
> fe_collection.push_back(FE_Q(deg));
> quadrature_collection.push_back(QGauss(deg+1));
> face_quadrature_collection.push_back(QGauss(deg+1));
> }
>
> I tried to look into the fe_collection.n_components(), it yields 1.
> is that what you meant?
>
> BR,
> Ihsan
>
&
Hi,
Lex, given your previous question:
https://groups.google.com/forum/#!topic/dealii/xghVE7Km78o
If you are already using an FEFaceValues object with a dummy-quadrature
(created from FiniteElement::get_unit_face_support_points())
then you can use one of the following functions on FEFaceValues
>>
>> Doug
>>
>> On Monday, June 22, 2020 at 5:05:57 PM UTC-4, Simon Sticko wrote:
>>>
>>> Hi,
>>>
>>> Lex, given your previous question:
>>>
>>> https://groups.google.com/forum/#!topic/dealii/xghVE7Km78o
>>>
>&g
Hi,
I'm not sure I understand your question correctly, but if you want to
impose an inhomogeneous Neumann boundary condition you would typically
write a small class representing the boundary condition you want:
class NeumannBoundaryCondition : public Function<3>
{
public:
double
value(const
t;fe_face_values.JxW(q_point));
>}
> }
>
> }
> }
>
>
> This seems correct to me, but the code is showing a cyldiner that looks
> like this:
>
>
> [image: Screen Shot 2020-06-29 at 5.50.19 PM.png]
> which doesn't appear to be th
nsigned int j = 0; j < dofs_per_cell; ++j)
> system_matrix.add(local_dof_indices[i],
> local_dof_indices[j],
> cell_matrix(i, j));
> for (unsigned int i = 0; i < dofs_per_cell; ++i)
> system_rh
Hi,
You'll need to call cmake with the flag
-DDEAL_II_WITH_UMFPACK = ON
> Also, do I need to reinstall dealii everytime I realise that I have not
configured dealii with a particular package?
Yes. However, if you need several of the optional dependencies, there is a
script that can make the ins
Hi.
There are several point_value functions. Right now, you are calling one that
only works for scalar elements:
https://www.dealii.org/current/doxygen/deal.II/namespaceVectorTools.html#a7be5c7eed52308898dfaad91c4cff204
You need to call the more general function:
https://www.dealii.org/curren
Hi.
From your attached files it looks like you have saved the matrix and rhs with
quite few digits. Do you get the same result if you increase the number of
digits you write to your std::ofstream object?
Compare:
https://en.cppreference.com/w/cpp/io/manip/setprecision
Best,
Simon
On 29/07/20
Hi.
I have done this in the past. I dug up the wrapper function I had for it, in
case we decide that we should add it to the library:
https://github.com/dealii/dealii/pull/15838
The problem with this approach is that evaluating the reference space level set
function gets very expensive. Quadr
ly would need to interpolate
(evaluate) my real-space function at this one element and let the
class' implementation do the magic for me. This is probably cheaper
and about equally accurate?
/Anton
On Fri, Aug 4, 2023 at 12:07 PM Simon Sticko wrote:
Hi.
I have done this in the past. I d
Hi.
The function
GridIn::read_msh
uses the Gmsh-api internally to read a mesh from file. However, Gmsh is
distributed under GPL and not LGPL.
Does this mean that one can not use this function if one wants to use
deal.II under LGPL?
Best,
Simon
--
The deal.II project is located at http://www
Hi.
1. Is there a way to deal make a correct sparsity pattern accounting for
hanging nodes, FENothing elements and coupling of faces that have ghost penalty
or is this really a feature that isn't implemented in deal.ii ?
If you do not care so much about performance at this stage, you can jus
Hi.
1. Is there a way to deal make a correct sparsity pattern accounting for
hanging nodes, FENothing elements and coupling of faces that have ghost penalty
or is this really a feature that isn't implemented in deal.ii ?
If you do not care so much about performance at this stage, you can jus
31 matches
Mail list logo
