ght solution.
Want to try writing a patch for this? We'd be happy to merge it!
Best
Wolfgang
--
----
Wolfgang Bangerth email: bange...@colostate.edu
w
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 list/forum options, see
https://groups.google.com/d/forum/deal
->is_locally_owned)
data_array(index)=value;
|
or
|
for(;cell!=endc,++cell)
if(cell->is_locally_owned())
data_array(cell->active_cell_index())=value;
|
Both will work. They are identical.
Best
W.
--
just say so, and that's ok with me as well.
Best
Wolfgang
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
T
oftware are fully coupled and the boundary conditions can be set using the
equations from the other software etc.
That would make for a fantastic code gallery project!
Best
Wolfgang
--
----
Wolfgang Bangerth email:
gang
--
----
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/f
ated?
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/
eved using the Vector::scale() function.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is
become very dense and difficult to invert.
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is
re.
You can't expect a high accuracy for a boundary integral if the boundary is
not accurately represented.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
.
--
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
st
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, see
https://groups.google.com/d/forum/
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
.
--
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
ystem?
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
ed:
https://github.com/geodynamics/aspect/blob/master/source/postprocess/heat_flux_statistics.cc#L122
It loops over all faces and then over the quadrature points of each face.
Best
W.
--
----
Wolfgang Bangerth
time?
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
, send an
email to dealii+unsubscr...@googlegroups.com
<mailto:dealii+unsubscr...@googlegroups.com>.
For more options, visit https://groups.google.com/d/optout.
--
----
Wolfgang Bangerth email:
r more options, visit https://groups.google.com/d/optout.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is
all, self-contained testcase that shows the problem?
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at htt
ke it smaller and smaller until the issue goes away.
The part that you removed last was then the offending statement. Put it back
in, and that is your minimal testcase.
Best
W.
--
----
Wolfgang Bangerth ema
there can be no neighbor) and so no constraints
will be computed.
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth
On 09/27/2016 10:12 AM, Deepak Gupta wrote:
Thanks, after this explanation, it sounds extremely silly on my side :)
You're not the first one to make this mistake.
W.
--
Wolfgang Bangerth
would be a
worthwhile 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 mailing list/forum options, see
http
ors
(Vector::symmetric_gradient), or point out any material
(manual, video etc.)?
Can you be more specific which functions and which parts within them you
have trouble with?
Best
Wolfgang
--
----
Wolfgang Bangerth email:
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/fo
ause the Raviart-Thomas
element mixes both the x- and y-components of the velocity vector.
I think you should be able to use
VectorTools::project_boundary_values_div_conforming
instead to do what you want to do.
Best
W.
--
------
ss.
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, see
https:/
case?
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.o
go through to make
this work in all generality :-)
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The
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, see
https://groups.google.com/d/forum/deali
ction. Several example programs use it, so you will find use cases there.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~ban
-l' instead?
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
specific reason for that recommendation, Wolfgang?)
No, it will work with any reasonably new version of Trilinos. Go with
the latest one.
Best
W.
--
Wolfgang Bangerth email: bange
f you had a small program that contained nothing but
int main() {}
and which you compile via
c++ test.cc -l/usr/lib64/libblas.so.3.6.1
? What if this program contained a call to a BLAS function?
Best
Wolfgang
--
Wol
ge manager
until you figure this out. At least you now have a small testcase.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu
a field
(with sufficiently high n_subdivisions passed to DataOut::build_patches()).
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http
omatic partitioning in step-18. But in any case,
even the automatic partitioning relies on METIS if you run it in parallel.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colosta
obal, not
cell->distribute_local_to_gobal, right?
Meanwhile I am going to extend the simple example for a basic version of
the elastic problem I am trying to solve. Hope then I can figure out the
error.
Good plan.
Best
W.
--
---
rix may be indefinite, and CG can not solve it.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II
.
--
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
CG may not converge.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://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 options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You
things are as you describe.
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at
re complicated recently.
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
t the
problem is, and then fix it the right way. Don't paper over it.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangert
s not* s.p.d. Using a different solver might
allow you to solve the linear system, but since the linear system is
apparently not the correct one, the solution will also not be correct.
Best
W.
--
Wolfgang Bangerth
uted, then
Pi_gamma = std::sqrt (2*eps*eps);
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is
https://github.com/dealii/dealii/issues/3091
Best
Wolfgang
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II proje
l you need to do at
the symmetry axis.
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II proje
hieve O(h^p), you
will of course have to ensure that the iteration error also becomes
smaller and smaller as you refine the mesh, i.e., you have to reduce the
Newton tolerance with every mesh refinement.
Best
W.
--
----
today than in
the early days!
Enjoy!
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located
a previous email, you ought to
make sure that the matrix you build is the matrix you *want* to build.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
on of type ‘dealii::Triangulation<3, 3>’
You pass a 3d triangulation to a function that expects a
2d-triangulation-embedded-into-3d-space as argument.
Which function are you trying to call?
Best
W.
--
----
Wolfgang Bangert
On Thursday, October 20, 2016 at 9:48:55 PM UTC-5, Wolfgang Bangerth wrote:
What is the norm of the first right hand side? You now set the tolerance in
the first iteration to a fixed value of 1e-16, but how does this compare to
the previous value of 1e-12*system_rhs.l2_norm
you wanted to try this, we'd be very happy to take your
program and incorporate it into the tutorial or the code gallery!
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
rs
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/
viously have to have MPI installed on your system.
You have two options: You can either install the MPI package on your machine,
or you can tell Trilinos not to use MPI. That choice is yours.
Best
W.
--
----
Wolfgang Bangerth
symmetric.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.deali
that returns an IndexSet that contains all of the
elements locally stored.
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu
o keep two meshes around.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is locat
only determined up to a constant). What you are
asking is whether the prescribed boundary condition makes sense *physically*,
but that is a separate question.
Best
W.
--
Wolfgang Bangerth email
is functionality, and we'd be happy
to walk you through it.)
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is
entation of programs such as step-40 to understand which elements of a
parallel vector are stored on each processor.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
alii/issues/391
See the second comment there. It *should* be a vector of Tensor<1,dim>,
but it's not, unfortunately :-(
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
you traverse the cells, one element
per cell. You can use that for further computations. You can even pass this
vector to DataOut.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
st work because it can use cmake's Makefile generator in the Eclipse
project.
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The de
is
not periodic.
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/
only one rank. What happens if you output from rank 1? Do you
get any output at all?
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http
he implementation of the quarter circle to see how this is done!
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.o
le?
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
will also need to be
recompiled then ?
I don't know. You say that it suddenly stopped working. Did it work
correctly before? What has changed in the meantime?
Best
W.
--
Wolfgang Bangerth
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 op
p the projection and use the DataPostprocessor class
instead. This is certainly more efficient than what you do right now.
Best
W.
--
Wolfgang Bangerth email: bange...@colostat
.
--
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
nd sample the solution
at these points to get upper bounds for the minimum, and lower bounds for the
maximum of the solution function u_h(x).
Best
W.
--
----
Wolfgang Bangerth email: bange...@
.
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.c
On 11/01/2016 11:06 PM, Jaekwang Kim wrote:
do you have any other idea?
Output the errors per cell and visualize. Is there a pattern? Is the error
localized on the boundary, for example?
Best
W.
--
Wolfgang
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
On 11/02/2016 07:06 AM, hanks0...@gmail.com wrote:
Anyway, what I'm trying to do is to get error = |exact_solution -
numerical_solution| (absolute value of two Deal.II Vector) and then plot this
one using Visit.
But I don't know how to change all elements of error(Deal.II Vector) to
absolute va
herewith.
My question would be: how do you determine whether you have the
correct/wrong solution? What do you compare with?
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
that to the
DataOut object you have. step-27 shows how to do that, but it in essence works
just like for vectors that have as many entries as there are DoFs.
Best
W.
--
Wolfgang Bangerth email:
n of step-6.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.d
to the "Possibilities for extensions" of one of the tutorials.
Want to give this a try and see whether you can come up with some code that
does that?
Cheers
W.
--
----
Wolfgang Bangerth email:
ather the information
on cells that are ghosts on some other processor, and then send this
information to all relevant processors.
Best
W.
--
Wolfgang Bangerth email: bange...@colostat
On 11/04/2016 07:50 AM, Victor Eijkhout wrote:
You guys don’t run into problems with trilinos?
No, I install it all the time. It may be that we don't because we don't
usually use STK...
Best
W.
--
----
Wolfgan
already.
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
dynamics/aspect/blob/master/source/simulator/checkpoint_restart.cc
This is of course a much larger project, but it may still help you understand
how things are used in practice.
Best
W.
--
----
Wolfgang Bangerth email:
t is
true for everyone, not just you! These accounts could use more activity :-)
Cheers
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colosta
:-)
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
ar to others.
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 options, see
https://groups.google.com/d/fo
riangulation and all nodal and element information.
How to go about it?
That is where you should use the boost serialization functionality others have
already pointed you to.
Best
W.
--
----
Wolfgang Bangerth
nts in setup_system?
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 l
.
That said, you *do* need to apply boundary values in subsequent Newton
iterations. They just have to be zero displacement. Take a look at step-15.
Best
W.
--
Wolfgang Bangerth email: bange
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/
F
in between?
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.deali
?
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
them?
I believe so -- there should be iterators for that.
In parallel, you can of course only iterate over elements that are
stored locally on a processor.
Best
W.
--
Wolfgang Bangerth email
this is a symptom of a lose cable or a wrong
configuration of the MPI system -- both things we have seen in cluster
in the past.
Best
W.
--
Wolfgang Bangerth email: bange...@colostat
II. This is not going to work.
Best
W>
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located
101 - 200 of 3019 matches
Mail list logo
