Hi,
I'm studying Wolfgang's lecture 31.65 that deals with picard iteration.
In the lecture, previous_solution_gradients[] is used. But, I also want to
know How I can use just previous_solution[].
At first, I used previous_solution[q_index] in assemble_system as in
previous_solution_gradients[
I've just updated to 8.4.2 and repeated the change to
trilinos_sparse_matrix.cc, and the matrices appear the same as my previous
post. So I've at least ruled out that it's not an issue specific to my
version number.
--
The deal.II project is located at http://www.dealii.org/
For mailing list/f
Hey there, thanks for the response. I made the recommended change to
trilinos_sparse_matrix.cc, and the corresponding test code returns the
appropriate output. However, I'm still getting an issue in my particular
case. It's improved in that the resulting matrix is now square and the
dimensions
Juan,
[...]
> 1) My first question is about how the "direction" works, is it 1
> correspond to "x" and 2 corresponds to "y"?
>
No, direction 0/1/2 means that the support_points of the dofs to be
identified just differ in the x-/y-/z-component.
What you are doing should be perfectly fine.
> 2)
Here it is:
andrew@andrew-apple-laptop:~/dealii-8.3.0/build$ ls -l /usr/lib/liblapack*
lrwxrwxrwx 1 root root 29 Sep 21 11:54 /usr/lib/liblapack.a ->
/etc/alternatives/liblapack.a
lrwxrwxrwx 1 root root 28 Nov 23 2015 /usr/lib/liblapack_atlas.a ->
atlas-base/liblapack_atlas.a
lrwxrwxrwx 1 root
On 09/22/2016 11:13 AM, Andrew Lambe wrote:
I'm trying to include the bundled UMFPACK package in my deal.II
installation, but when I switch the options DEAL_II_WITH_LAPACK and
DEAL_II_WITH_UMFPACK on, I get a weird error from CMake that the LAPACK
library is not detected. I say weird because I ha
Chenchen,
If you just need to simulate a multlayer sphere (or a part of it) and the
bounding box should not have to be a box (e.g. a free space problem) you
can also take a look at my attempt for this
here https://groups.google.com/d/msg/dealii/pvqNQUM3eyM/ftqQK1ypBQAJ Just
set all raddii of t
Chenchen,
2016-09-22 10:54 GMT-04:00 Chenchen Liu :
> Thank you for your sharing. To my mind, after we generate Q2 mesh in Gmsh,
> we have all the information about vertices and cells. We can read these
> information to deal.ii as point list and cell list, and then generate mesh
> like step-14 and
Thank you for your sharing. To my mind, after we generate Q2 mesh in Gmsh,
we have all the information about vertices and cells. We can read these
information to deal.ii as point list and cell list, and then generate mesh
like step-14 and 49. But still, I am not sure whether deal.ii supports
hi
I see your point. Thank you! But it seems that hanging nodes is really
annoying :) I will have a try.
在 2016年9月22日星期四 UTC-4上午7:38:35，Konstantin Ladutenko写道：
>
> BTW, the element in bottom left corner is too triangle for good FEM. To
> solve this providing an additional element inside the curved
Thank you! I almost done with this.
在 2016年9月21日星期三 UTC-4下午12:04:19，Bruno Turcksin写道：
>
> Chenchen,
>
> On Wednesday, September 21, 2016 at 11:56:31 AM UTC-4, Daniel Arndt wrote:
>>
>> We only support to read in meshes that are not refined. Translated to
>> your situation this means that you are
Dear Luca
I wanted to use MappingManifold for sphere, but it looks like this is not
possible either with PolarManifold or SphericalManifold.
I am doing this for a shallow water model on sphere, where I want to do
long time integration, and look at conservation of mass, energy, etc.
Having an exac
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
BTW, the element in bottom left corner is too triangle for good FEM. To
solve this providing an additional element inside the curved boundary looks
to be a common approach, you can find more details
here http://dealii.org/developer/doxygen/deal.II/step_6.html#Abettermesh
среда, 21 сентября 201
Hello
This is a bit old question, and I would like to pose it again. The first
post is here
https://groups.google.com/d/msg/dealii/pB-s0JELJy0/2vBXef_IBQAJ
I am using PolarManifold and a 90 deg rotated PolarManifold to cover the
surface of sphere.
As I posted in this thread, some cells at the c
Thanks Luca. I have been looking at tests also.
For CompositionManifold, if
H = CompositionManifold(F, G)
then
H_pullback = F_pullback * G_pullback
H_pushforward = G_pushforward * F_pushforward
I think adding this to documentation will improve the description.
Also this line
dim1 The dimens
16 matches
Mail list logo
