Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing list/for
.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum
VERBOSE=ON
and it should tell you the command lines being used.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
T
at
https://github.com/dealii/dealii/wiki/Frequently-Asked-Questions#im-stuck
Extra question: Would my work be of interest for your code gallery maybe in
the near future?
If you have good questions, always!
Best
W.
--
----
Wolfgan
the Laplace equation, but if you see it
happen there, you will know how it works for the elasticity equation as well.)
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
ing the local
right hand side you may want to see steps 3 and 4, as well as 8. For how
exactly this interacts with hanging nodes, you may want to see step-6.
The point is that you need to expect to see similar mathematics, but
different terms than you are used to.
Best
W.
--
-
ere is truth to that. So,
*unless you have measured that it is a problem*, then the simplest approach is
the way to go.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
iner mesh -- maybe a
triangular mesh, whereas the other one is computed on a quadrilateral mesh.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://ww
tion.
What version of the TBB are you using? And can you look up the declaration of
the spawn() function in task.h in your version of the TBB?
Best
W.
--
----
Wolfgang Bangerth email: b
refine.
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http
uch appreciated!
Best
Wolfgang
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.
:-)
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing list/
Or wait for the patch I'm going to submit in a minute :-)
Specifically, I would appreciate if you could see whether this helps:
https://github.com/dealii/dealii/pull/3889
Best
W.
--
Wolfgang Bangerth
Is
this in debug mode?
Best
Wolfgang
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options
never impose any pressure
boundary conditions.
There is a description of boundary conditions for the Stokes equations
in the introduction of step-22. Take a look there to see what you can
and cannot impose for the Stokes equations.
Best
W.
--
---
ve to multiply it by invJ
again. This is the source of your factor of two.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.ed
columns of that operator.
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is locat
aside, both the x and y components of
the velocity are 0.0 ON the boundaries and 1.0 everywhere else.
Are you calling constraints.distribute on the solution vector after
solving your linear system?
Best
W.
--
Wolfgang
appen if you
compile with -j1? The kind of internal compiler error you have here
seems to happen most often if you run out of memory.
Best
W.
--
Wolfgang Bangerth email: bange...@colostat
er than its source
files) but can't read successfully.
The solution is to do 'make clean' or just to blow away the build
directory and start again.
Best
W.
--
----
Wolfgang Bangerth email
/dir/, which is a path you want to
replace by the one where you want to install deal.II in :-)
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http
mesh plus
one layer of ghosts, and the derived class is responsible for the
communication? And then derive another class from DoFHandlerPolicy to deal
with this triangulation?
Best
W.
--
----
Wolfgang Bangerth
0sec output, and 80 sec
create_tria, and 45 sec assembly. Each core requires 9GB memory. This is why I
want to reduce memory usage.
Yes, I can see why the reorder step is annoying then.
Best
Wolfgang
--
----
Wolfgang
want to multiply with. (Hint: The correct component is
fe.system_to_component_index(i).first).
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http
s why that is so -- is your rhs correct, for example?
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
T
.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
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
ed, but take a look at lectures 21.6 and 21.65
here:
http://www.math.colostate.edu/~bangerth/videos.html
Best
Wolfgang
--
----
Wolfgang Bangerth email: bange...@colostate.edu
e purposes,
but you need to know what you are doing.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.I
set of delta functions).
I also attach the plot for d^2(solution)/dx^2 that look weird to me...
I think it is plausible. It is at least not obviously wrong in my mind.
Best
W.
--
Wolfgang Bangerth email
read.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum opti
a default value), or indeed a positive integer.
I'm going to note that figuring out that a string is indeed a positive
integer (and nothing else) is not trivial.
Best
W.
--
Wolfgang Bangerth email:
that's where I
would look things up.)
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II projec
imately (and caught), and in those cases just
`continue` in gdb, until gdb stops in the place you care about in fe.cc,
line 859.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
By adding catches on assertion,
It took me a few years until I learned this trick (the `catch throw`
thing) but it turns out to be really useful.
Cheers
W.
--
Wolfgang Bangerth email: bange
ument you provide.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://ww
invertible. This, as well, can be tested with
a direct solver.
Go through this list of ideas and you'll be closer to finding the root
cause!
Best
W.
--
----
Wolfgang Bangerth email: bange..
et this to work.
Can you open an issue on the github forum at
https://github.com/dealii/dealii/issues
It would be great to have a minimal example that shows the problem.
Best
W.
--
----
Wolfgang Bangerth em
ate the 'build' directory?
Parallel to the directory into which you unpacked deal.II?
You will also want to replace /path/to/install/dir by a directory of your
choice.
Best
W.
--
----
Wolfgang Bangerth
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing
n be a second bug that it doesn't use the generalized support
points. We'd take a (separate) patch for that as well.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
ve elements. But the patch is definitely not going to make anything
worse. Give it a try (in a separate pull request).
Cheers
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
plot
(neither of which are equal to the P2 function in deal.II).
You can probably get a better picture by providing a positive argument to the
build_patches() function.
Best
W.
--
Wolfgang Bangerth
is used internally in deal.II,
and (ii) what Visit plots is not what you output. But with higher arguments to
build_patches, you get closer in both regards.
Best
W.
--
Wolfgang Bangerth email
d 1d things?
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
of
unknowns?
Best
W.
PS: I do like this approach to testing. It shows great maturity in designing
code and testing numerical methods to think this way. Well done!
--
----
Wolfgang Bangerth email:
or a penalty term) to force one
function to be constant in y-direction.
But why care? 2d computations are so cheap these days. Can't you just
assume that everything depends on x and y?
Best
W.
--
----
n your
equation.
And how do the two meshes differ?
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal
gy norm error of the respective
solutions should be better, but we don't know whether that's true for
the L2 norm.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
read user indices. We do read material_ids, however.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal
h should not be
such a big "deal" in deal.II.
Because it's *not* easy. It's only easy if you consider simple situation
like yours. But if you want to do complicated things, they generally
turn out to be complicated.
Best
W.
--
--
ximum among all processors, for example using Utilities::MPI::max().
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
for (l=0...) {
...
output << w << ' ' << val << std::endl;
}
and then visualize the resulting data using gnuplot's "plot" command.
Best
W.
--
Wolfgang
gested
here, and see if you can get a backtrace? Alternatively, can you create a
small, self-contained testcase that shows the problem?
Best
W.
--
----
Wolfgang Bangerth email: bange...@co
rce the constraint, due to time
constraints.
But there is an H-div conforming projection function in the same namespace
that should work for your case.
Best
W.
--
----
Wolfgang Bangerth email: bange...@
the wrong place.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http:/
id, do you end up with a Makefile file in your build directory? If so,
I suspect you can just ignore the error.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
are weak, you can just omit the term
in the assembly.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at
not there.
I would suggest to go the route that Uwe suggested: use candi.
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu
at you are trying to do here may or
may not make sense for the equations you are considering.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.m
thing else that may seem appropriate.
Thanks!
Wolfgang
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II proje
st
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/
, the resizing is unnecessary because the system has the same
size. But you could try to just set matrices and vectors to zero, and see
whether that makes a difference.
Best
W.
--
Wolfgang Bangerth
That doesn't trigger anything for me :-(
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II projec
polynomial arc.
A similar effect is mentioned in one of the tutorial programs that deal with
mappings (9 and 10? or 10 and 11?).
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
m is
going to consume at least as much memory as the 450k artificial cells.
So I think that 500k coarse cells is sort of the limit of what makes sense,
but it's not in itself a death sentence yet that derails the whole approach.
Best
W.
--
---
27;t understand what it
computes. Learning strategies to debug software is the only way you can learn
to write good software.
Best
WB
--
----
Wolfgang Bangerth email: bange...@colostate.edu
reduce.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing lis
e to set it up once, but that is again a skill that you will need
one way or the other at one point as your programs become more complicated.
Best
W.
--
----
Wolfgang Bangerth email:
this happens to work, please submit this as a patch in general -- we
should try to preserve error messages.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
know that
I'm going to spend an afternoon writing stuff on a piece of paper and
comparing with what I get, but it's also often helped me debug things
with which I've banged my head against the wall for too long already.
Best
W.
--
----
may go wrong in the process.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is loca
question and set
u_i = u_42
This is a straightforward constraint to add to the ConstraintMatrix class.
Take a look at step-11 to see how constraints are added.
Best
W.
--
Wolfgang Bangerth email
ution might
different when mesh is refined
is there anyone who has idea on this ?
Take a look at step-15 -- it does exactly this, using the SolutionTransfer
class.
Best
W.
--
Wolfgang Bangerth
does not exist. Did you re-initialize the sparsity pattern and matrix after
you refined the mesh?
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http
ure on as coarse a mesh as you can use to demonstrate
the error?
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
a 2x2 or 3x3 mesh, to make sure it really is because of the
constraints? For small problems, you can often figure out on a piece of paper
what the matrix should be, and compare with what you get from your program.
Best
W.
--
-
ry.
And I solved AU_0=F-AG_tilt by...
Vector tmp;
tmp.reinit(solution.size());
system_matrix.vmult(tmp,G_tilt);
system_rhs -= tmp;
solve ();
solution += G_tilt;
This looks correct.
Best
W.
--
---
.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options,
taOut output this.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d
node constraints? Take a look how that is done in step-6, for example.
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~ban
hen convert the result of that function into a vector.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The dea
.
It's a well done podcast I'm happy to recommend if any of you have long
commutes during which you like to listen to something.
Cheers
W.
--
----
Wolfgang Bangerth email: ba
.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
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
uot;. Can you elaborate?
Cheers
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum optio
his is issue 4092 now.
Great, thanks. I can reproduce this, but I don't know when I (or anyone else)
may be able to get around to fixing this :-(
Best
W.
--
----
Wolfgang Bangerth email:
h you only want to simulate in r-z space?
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II proje
edges via the
boundary_id tag. The function you call only works with the latter. But you use
the former. You will want to use the HyperBallBoundary class for the boundary.
Best
W.
--
----
Wolfgang Bangerth email:
o.h:293:7:
/**/note: /*
/ 'RefinementCase' declared here/
/class RefinementCase : public RefinementPossibilities/
You should read through step-30 to see how this is done.
Best
W.
--
----
Wolfgang Bangerth email:
ration.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.or
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing list/for
ingQEulerian
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing li
, and why that entry has not been allocated before.
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II
others that are also good.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is
LinearAlgebraTrilinos::MPI::BlockVector class.
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II
see how
this function is supposed to look.
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II
sh
step-47. Since Thomas Carraro essentially finished the program in the ANS
article, I may even propose removing step-47: it's never going to be finished
as far as I know, and a completed version exists.
Best
W.
--
----
W
where this would fit? For the ellipse,
maybe as an example in the documentation of the ChartManifold class?
We'd be happy to add it there or somewhere else!
Best
W.
--
Wolfgang Bangerth email:
mance.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mai
ks
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing
.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d
301 - 400 of 3019 matches
Mail list logo
