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
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
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
o 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 dug up the wrapper fu
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.
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
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:
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,
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
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
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
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
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
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,
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
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_rhs(local_dof_indices[i]
int));
>}
> }
>
> }
> }
>
>
> 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 the right solution. Even when I change
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
gt;> 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
>>>
>>> If
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
{
> 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,
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
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:
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
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
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
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
27 matches
Mail list logo
