what happened when you
first built the matrix and the original right hand side.
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http
to see how far it runs, or (better!) by running the program in a
debugger.
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu
a deal.II-based program or
later on, it's not manageable one way or the other.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.ed
d consider the resulting matrix as
block diagonal.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II
(history_dof_handler.n_dofs());
}
This does not look crazy, unless you have a large number of degrees of
freedom. How large is history_dof_handler.n_dofs() in your case?
Best
W.
--
Wolfgang Bangerth email
uto i : locally_owned_dofs)
outloop << testvec(i);
or similar.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
mentioned before: If you
have a large number of processes (say, 1000), then you have one process
doing a lot of work (1000x as much as necessary) and 999 doing nothing.
This is bound to take a long time.
Best
W.
--
Wolfgang
if you wrote a patch to add the copy constructor!
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II pro
chine with more memory, (ii) use distributed computing where the 6M
unknowns are distributed across more than just one machine.
Best
W.
--
----
Wolfgang Bangerth email: bange...@col
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,
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
s does not have a unique solution: If u(x,y)
is a solution, then u(x,y)+c is also a solution for any constant c. You
shouldn't expect the solver to provide you with a solution.
Best
W.
--
----
Wolfgang Bangerth
ve to
assign that property the value of your solution at that point (or
whatever other quantity you want visualized at this point). The final
step is to use the particles::DataOut class to output the information at
the particle locations.
Best
W.
--
-
data during
re-partitioning. Data transfer is separate from partitioning in deal.II.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colos
can I?
I am pretty sure that is no longer correct. Can you check in step-19 or
one of the other particle programs?
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
about the case here.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http
tions that can translate
in both directions.
I hope this helps!
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
-
;
spacedim > > > & ,
is it?
If you're working on a sequential triangulation, you just pass an (outer)
vector of length 1. (I believe at least.)
Best
W.
--
----
Wolfgang Bangerth email:
Finite elements and fast
iterative solvers". The canonical book on iterative solvers is the one by
Saad, but it only has a few pages about block methods.
Best
W.
--
----
Wolfgang Bangerth email: ba
r objects to use the same triangulation. But
knowing the implementation, it would make sense to me if they did. If that
were true, then the last question above would has "yes" as the answer.
Best
W.
--
----
Wo
GMSH writes. If you created the mesh in the x1-x2
plane, then those will be the coordinates in Triangulation<2>.
Or maybe I don't understand the question?
Best
W.
--
----
Wolfgang Bangerth email:
s I have in the
email you quote. Either do everything within Trilinos, or don't use
Trilinos at all if you want to hand your matrix to some other package.
Mixing and matching is always difficult... "
Best
W.
--
-
wn to a smaller part of the code, and the
questions above should help you with it!
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate
num is declared inside KINSOL::AdditionalData.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II proj
IDs.
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
using a
backtrace) and why that triggers the assertion. Knowing how to work with
a debugger is a skill you will find is useful for the rest of your life.
Best
W.
--
----
Wolfgang Bangerth email:
various
project() functions to see whether one of them works in parallel. If
none do, you can always use VectorTools::interpolate().
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
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
f. If you'd like to explore this
and write a patch, let me know and I'll walk you through how one would
approach that!
Best
Wolfgang
On 9/9/22 01:08, Abbas wrote:
*** Caution: EXTERNAL Sender ***
Yes
On Friday, September 9, 2022 at 12:06:23 AM UTC+3 Wolfgang Bangerth wrote:
ogram in a debugger and see where the segfault
happens by looking at the backtrace. Then look around what expression causes
the segfault, and 75% of the time you will already see what the problem is.
Best
W.
--
----
Wolfgan
to
each.
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:/
.
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
state.edu%7Cfbdc82ddb44644ad4f1c08da9ea1fb5d%7Cafb58802ff7a4bb1ab21367ff2ecfc8b%7C0%7C0%7C637996715835974004%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=5%2BXzdxaM3tXgOcs7LHyQXFK4fog6nXaSktC%2F%2Bzq%2FdlM%3D&reserved=0>
--
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
er it failed on GMRES solver. Does deal.II built-in GMRES solver
only supprts double type?
Yes. You probably want to take a look at step-58 as well, along with the
discussion in the results section there.
Best
W.
--
----
Wo
PETSc with the flag necessary to support complex arithmetic?
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is
kes if
you work in complex, rather than separate real/imaginary, arithmetic. But it
can be done correctly, as shown in step-58.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.
tructure you find convenient, but it can also be done
with standard MPI calls.
This isn't particularly convenient, but it is the best anyone can offer short
of writing the necessary functionality in PETSc itself.
Best
W.
--
----
example.
Best
Wolfgang
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.
tion.
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
, the latter is called GridGenerator::merge_triangulation().
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The
/dealii/dealii/wiki/Frequently-Asked-Questions#dealii-programs-behave-differently-in-1d-than-in-23d
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http
.
--
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
s well.
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/f
stablished researchers who seek geographic relocation or temporal
employment.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangert
f the text surrounding step-44, would you be interested in extending the
documentation of that program?
Best
Wolfgang
--
----
Wolfgang Bangerth email: bange...@colostate.edu
point number.
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
want to look at the step-79
that does something similar.
Best
W>
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
, to a vtk file.
Did you use .setprecision(...) also on the stream that you give to
data_out.write_vtk()? You need to set the precision separately for each
stream on which you want to output.
Best
W.
--
Wolfgang Bangerth
un] Aborted (core
dumped)
make[2]: *** [CMakeFiles/Makefile2:116: CMakeFiles/run.dir/all] Error 2
make[1]: *** [CMakeFiles/Makefile2:123: CMakeFiles/run.dir/rule] Error 2
make: *** [Makefile:137: run] Error 2
Is it because I am applying Dirichlet bc using two objects, constraints and
boundary
dure.
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
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/dealii?hl
of the current J^T J once.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is
on A^*, B^*,
you can back out to x^*, y^*.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project
to error that "point is not on
manifold" even tolerance is 10 times of original iges tolerance.
In 3d, we first refine the edges and then the faces. You need manifolds for
both 1d and 2d objects.
Best
W.
--
-----
e is that
the scaled version improves (amongs others) the condition number of J and may
lead to
a better convergence of the optimsation algorithm, right?
Yes.
Best
W.
--
----
Wolfgang Bangerth email:
s
Wasim Niyaz
Research scholar
CE Dept.
IITM
On Tue, 11 Oct, 2022, 7:03 pm Wolfgang Bangerth,
wrote:
Wasim,
I imagine you are calling your setup_boundary_values() function before
you
call dof_handler.distribute_dofs()?
Best
l) or cell->face(f)->line(l) gives you what you want.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
-
the other FE field, and then call
cell_2->set_dof_values()
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
-
be a different story.
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at
s are kept.
This exceeds my knowledge of meshing and geometry software. Maybe others here
can help with those questions?
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
,
and also
> fix the lower
> > > right node in y.
> > > The lower right node is free to move in x.
> > > The problem I am facing is that both these nodes are
on the
> same face,
>
.
--
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
should be using
FEPointEvaluation instead.)
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
.
--
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
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
r next. I mainly wanted to let you
know that the problem is not forgotten.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.e
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
terms what your goal is?
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at
question.
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
compile; it does not
imply that the library exists on the system on which you want to do
computations.
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http
Uxk%2F2nsPd2CuK%2FyGJ7LrZTB6CFLod9UdkcTIzCff4%2Bk%3D&reserved=0>.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is locate
--
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
t is
already fixed.
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/
r the computed solution is correct?
Best
W.
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is locat
should work.
Best
WB
--
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http
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
to read individual vector elements also for ghost entries, and
so then you need the locally relevant dofs. Vectors with ghost elements are
read-only.
Best
W.
--
----
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,
And at the other functions that are in that namespace, for example the
hyper_rectangle function mentioned in the documentation linked to above?
Best
W.
--
Wolfgang Bangerth email: bange
ny other format?
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 mailin
out ("my_output_file.vtk");
out.precision(16);
data_out.write_vtk(out);
I haven't tested this, but believe that it should work.
Best
W.
--
----
Wolfgang Bangerth email: bange..
a specific file format,
you have to state which one you mean and how it is defined.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.mat
--
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
ular.
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/dealii
mpile everything.
Alternatively, you can install deal.II through your package manager if you are
on a supported system, and then the deal.II installation should just work.
Best
W.
--
Wolfgang Bangerth
/science/article/pii/S0997753822002753
If you want a pdf version of it, you can find it here:
https://arxiv.org/abs/2209.04198
We hope that the paper is useful to many of you as well!
Best
Wolfgang
--
----
Wolfgang Bangerth email:
es or why, having an understanding
what part of the function it is is a good start.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.c
of it, you can find it here:
https://arxiv.org/abs/2209.04198
Best
WB
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~banger
e_vtk(...)). That function does not support tensor output,
either because we didn't implement it or because the VTK file format
doesn't support tensors -- I don't recall. In any case, switch to
write_vtu() instead and everything should work.
Best
W.
--
---
me polynomial degree on all cells, this should not
even be particularly difficult and I would be happy to walk you through
the necessary steps if you were interested in writing up a patch!
Best
Wolfgang
--
Wolfgang Banger
yet, so this will have to wait
for a few years :-)
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is lo
.
--
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
of us, so we don't really know what is specific about
your situation to help. Can you run the program in a debugger and obtain a
backtrace?
Best
Wolfgang
--
----
Wolfgang Bangerth email: bange..
ndary conditions. The form of this
integral does not depend on whether you are in 2d or 3d.
Best
W>
--
----
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colos
const doublephi_i
= fe_face_values[displacement].shape_value(i,q);
rhs_values[i] += (n*p*phi_i) * fe_face_values.JxW(q);
}
... copy local to global ...
Does this make sense?
Best
W.
--
-
way as the computation of the stress tensor is currently in the
documentation of DataPostprocessTensor.)
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
unction.
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
uation you care about.
The final step is to translate that in the same way as what the program
currently does.
Best
W.
--
----
Wolfgang Bangerth email: bange...@colostate.edu
ient.
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/
901 - 1000 of 3019 matches
Mail list logo