individual header files to the
compiler, just include paths. So if the compiler can find one Trilinos header
file, it should find them all.
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
th of them, and that
guarantees that the elements of vectors match.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangert
oking forward to the patch! :-)
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is locat
index() const;
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/foru
).
Sticking with the old release would require you to write work-arounds for
functionality that exists today, and to stick with these work-arounds for the
indefinite future.
Best
W.
--
Wolfgang Bangerth email
tutorial programs (e.g., step-9,
step-21, or step-15's and step-19's use of the corresponding
get_function_gradients() function).
Best
W.
--
----
Wolfgang Bangerth email: bange...@col
Here you multiply tensor times vector times tensor. It's not clear to me what
that is supposed to represent, but that's at least where the error is from.
Best
W.
--
----
Wolfgang Bangerth
rs!
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
On 5/31/21 9:36 PM, Praveen C wrote:
It this how it should be used ?
Something like this. You might want to look at step-47 which assembles face
terms for interior faces.
Best
W.
--
Wolfgang Bangerth email
lf in your home directory and take it from there?
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II projec
.
--
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
ering internally -- i.e., there is
no need to do that on the user side.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~ba
contiguous in memory? And why do you care
about the cell-local order of DoFs, and not the global DoF indices?
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
renumbering scheme is called anyway?
Correct.
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located
index );
++index;
}
There is no need for the whole *local_dof_indices*, if I know the starting
index on cell, it is enough.
But you can achieve this using global DoF index renumbering. That's the way
I'd go.
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
rdless, I think that the keywords you are looking for are
"checkpoint/restart". If you search the archives of this forum, you will find
numerous mentions of this topic.
Best
W.
--
----
Wolfgang Bangerth
s the correct solution to the wrong
linear system.
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II pr
r we can't fix it anyway!
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/
Fo
.
(You could also look at the new step-72 and step-77 programs to see how
certain aspects of this could be simplified.)
Best
Wolfgang
--
Wolfgang Bangerth email: bange...@colostate.edu
admit that I have no intuition.
Do both schemes converge? To the same solution?
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http
ould you use FE_FaceQ?
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
around node i.
Best
W>
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at h
rm -r the entire build directory before you call cmake and see whether
that helps.
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www:
, i.e.,
||sigma|| / ||eps||
does not go to zero or infinity, and is a smooth function of some sort.
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www:
.
--
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
ave zero boundary values, for example, then you'd just set
all DoFs to zero in that loop.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
;${_IMPORT_PREFIX}/lib/libsacado.so.12.8.1" )
Binary file libsacado.so matches
Binary file libsacado.so.12 matches
Binary file libsacado.so.12.8.1 matches
Best
W.
--
----
Wolfgang Bangerth email:
"before" and what is "after" this change?
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
DoFs
* Ones that have as many entries as there are active cells
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth
.
--
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
embles a rhs vector. As just one example, the new
step-77 program has a function compute_residual() that does essentially this.
Best
Wolfgang
--
----
Wolfgang Bangerth email: bange...@colos
a try -- it would solve your problem, and would be very useful
to others as well! We'd love walking you through the steps to make that happen!
Best
W.
--
----
Wolfgang Bangerth email:
On 6/10/21 7:10 AM, Wolfgang Bangerth wrote:
Indeed,
std::vector> velocity_qpt_vals(n_face_qpts,
Vector(dim));
fiv.get_fe_face_values(1).get_function_values(velocity_field,
velocity_qpt_vals);
accesses the values; the issue was I hadn't explicitly instanti
tions on faces
instead of edges/vertices.
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II proje
.org/participate.html
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
this. step-11, for example, constrains one node
as a linear combination of (many) others.
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http
things would further improve.
You might also want to try quadratic or cubic elements. For smooth solutions,
that often does wonders!
Best
W>
--
----
Wolfgang Bangerth email: bange...@colostate.
_2);
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 option
cell-based
vectors, we ignore the fact that the computation might be parallel.
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http
exists.
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
p-61, which does the kind of thing you're looking for, I think.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~b
atial
derivatives. How do you deal with this?
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II proje
close a list of pairs, which here has only one
element and that element is initialized by the pair u,v.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
ww
ata points what the rate is,
and it may well be possible that you can't ever observe the true rate of the
algorithm because you can't solve problems so large that the asymptotic rate
is attained.
Best
W.
--
-
omponents.
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
not.
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://gr
ther libraries deal.II depends on. All of
this is information that cmake exports and that your project would
automatically import if you used cmake.
This is wheel already invented. Don't waste your time trying to re-invent it :-)
Cheers
W.
--
--
d programs because all it requires is the solution of linear
systems -- something that you have likely already implemented.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
lator/nullspace.cc#L112-L209
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at
its degrees of freedom are, and then query the solution vector for these
components. Or just call
cell->get_dof_values(...)
right away for the values of degrees of freedom located on a cell.
Best
W.
--
----
Wolfgang Ba
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/forum/dealii?hl
be an easy fix, and a good first patch to
contribute to the library!
If you're curious about ho to write a patch, take a look at
https://github.com/dealii/dealii/wiki/Contributing
Best
W.
--
Wolfgang Ban
mp_points
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
oach using deal.II for your application very much
depends on the equations you want to solve (more than on the physics that
these equations describe).
Best
W.
--
----
Wolfgang Bangerth email:
%40colostate.edu%7Ce345d6aef515473cad2d08d93a3eb17f%7Cafb58802ff7a4bb1ab21367ff2ecfc8b%7C0%7C0%7C637604863211998906%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=XU5kLsJaryhDRV0CvR118q%2BrfgU0CjorvGPLFu3xhn4%3D&reserved=0>.
--
--
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 op
.
--
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
the
CellDataTransfer class, though I have no idea if anyone has ever tried that.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.ed
nt
to do.
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://group
mentations are
based on.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is l
t have written this
as
J(un) delta u_n
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.deal
h manifolds to the triangulation object that is being
reconstructed, before reconstruction.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colos
omehow to actually
connect a member function of Solid?
No, but you need to provide a 'this' pointer for the object you want to work
on. In the first case, this is pointer_M.get(). For a member function of the
current object, it would probably be 'this'.
Best
ne way or the other.
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.g
and
"captured variables". You might want to read up on them in a recent C++ book :-)
Cheers
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project
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 ma
e cells, including
ghost and artificial cells (for which vector entries are then just ignored). I
think step-47 shows this (whether in parallel or not doesn't matter in this
regard).
Best
W.
--
----
Wolfgang Bangerth
nt it out.
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/
Fo
rogram and add it to another one.
Best
Wolfgang
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project
!
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
in submitting a patch?
The right approach is to simply use the corresponding DoFHandler iterator type.
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www
dler());
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/
ectorTools::interpolate.
The problem is that FEFieldFunction is expensive. If you don't care about
speed, this is the way to go. If you do for whatever reason care about speed,
go with the approach in the previous paragraph :-)
Best
W.
--
-------
's behind
this is TriaActiveIterator, which indeed has the constructor you need here.
Best
W>
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math
point values that visualization programs then interpret
via bilinear/trilinear interpolation. There are numerous tutorial programs
that use DataPostprocessor. A particularly simple example is step-29.
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, see
https://groups.google.com/d/forum/de
.
--
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
oFHandler::active_cell_iterator
cell_2 (&cell_1->get_triangulation(), cell_1->level(), cell_1->index(),
&dof_handler_2);
Best
W.
--
----
Wolfgang Bangerth email:
h a way to
get what you are looking, assuming you really are only using Q1 elements.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.col
linear hex Q1 elements always?
Yes.
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is lo
or in which you just ignore whatever it gets as output,
and instead proceed to use as output whatever you query from your
std::vector> using the cell the evaluate_* function gets as
argument in the input object.)
Best
W.
--
----
tep will take many hours on a single-processor machine. If you have
enough cores, you can do
ctest -j80 -R debug
and it will finish in about 30 minutes, but you really do need 80 cores for
that.
Best
W.
--
----
Wolfgang
you have concrete evidence
that whatever performance you get, going with the easiest-to-implement
approach is always the right choice!
Best
Wolfgang
--
Wolfgang Bangerth email: bange
or of second derivatives
and sum over the appropriate set of elements.
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~ban
the second half of the pair.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://
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/f
nly one one side of the interface. The other
values (22116; 89130) shouldn't happen. If that's what you get, then there is
a bug somewhere.
Do you think you could construct a small testcase that illustrates what you
have?
Best
W.
--
---
e using DG elements. I tried to clarify this here:
https://github.com/dealii/dealii/pull/12595
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
p::d::T on
a single processor.
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is
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
ements].get_function_gradients(...);
i.e., you need to use an extractor in [square brackets].
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www
ck out the Utilities::MPI::RemotePointEvaluation class in
the latest (9.3) release for help with this.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
n memory?
Correct. They don't change.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II proje
to worry about "evaluating" the solution?
Yes, you can reload solutions onto the same Triangulation and DoFHandler. In
fact, the reload machinery even works if the number of processes is different.
Best
W.
--
---
in your .bashrc or, if you just
want to do it once, on the command line before running the program.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www:
fgang
--
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?
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/forum/dealii?
() function where that issue is even documented. You're
right that in principle, one would want to transform the manifold along with
the mesh, but nobody had an idea of how to do that in practice.
Best
W.
--
Wolfgang Ban
then output these values by hand.
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is
501 - 600 of 3019 matches
Mail list logo
