Krishanu
Le mar. 16 juil. 2019 à 16:00, Krishanu Sen a écrit :
>
> I am using deal.ii for a nonlinear elasticity problem. And the change in the
> code that is causing this error message is:
>
>
> template
> ElasticProblem::ElasticProblem ()
> :
> dof_handler (triangulation),
>
Hello Bruno,
Here is the error message that I am getting:
An error occurred in line <388> of file
in function
dealii::Tensor::value_type& dealii::Tensor::operator[](unsigned int) [with int rank_ = 2; int dim = 3;
Number = double;
Krishanu,
It's impossible to help you with so little information. Somewhere you set
something to the wrong size. Can you give us a backtrace? We need to see
the error message and the code that produces it.
Best,
Bruno
On Tuesday, July 16, 2019 at 3:08:42 PM UTC-4, Krishanu Sen wrote:
>
> I
Ramprasad,
On Tuesday, July 16, 2019 at 11:51:00 AM UTC-4, Ramprasad R wrote:
>
> The version of deal.ii is 9.0.1 and the version of PETSc is 3.10.5. I also
> changed the config file of deal.ii to have MPI PETSc on.
>
What do you mean by "I changed the config file of deal.II"? You should not
I am trying to run a code using higher order elements. The code runs fine
when I am using linear [FE_Q(1)] or quadratic [FE_Q(2)] elements.
But when I tried to use cubic elements [FE_Q(3)], the code stops with
an error "ExcIndexRange", which says that the code is trying to access an
index of a
A short comment:
On Tue, Jul 16, 2019, at 11:18 CDT, Matthias Maier wrote:
> struct LeftVector {
> };
> struct RightVector {
> };
These two classes are of course just decoration (showing the minimal
interface a vector has to possess). There is usually no need to define
custom Vector clases.
Hi,
the following example (which is a simplified version of test
lac/linear_operator_01) should get you started. In short: Simply create
an empty LinearOperator object (with the appropriate template
parameters) and populate the corresponding std::function objects.
Best,
Matthias
#include
I tried to follow the examples for setting up a LinearOperator (such as
shown in Step-22 and Step-20), but instead of a matrix I wanted to provide
my own function which should serve as vmult-function (i.e. a matrix-free
linear operator). Is that possible? My initial test was the following code:
Dear Daniel,
The version of deal.ii is 9.0.1 and the version of PETSc is 3.10.5. I also
changed the config file of deal.ii to have MPI PETSc on.
Thank you.
Best regards,
Ramprasad
On Tue, Jul 16, 2019 at 5:12 PM Daniel Arndt wrote:
> Ramprasad,
>
> Let's first rule out some simple problems.
If I understand it correctly, that is equivalent to solving Ax = b, with A
and b given, and x returned?
Am Montag, 15. Juli 2019 22:03:46 UTC+2 schrieb Daniel Arndt:
>
> Maxi,
>
> The interface for TimeStepping::evolve_one_time_step is:
>
> virtual double
> evolve_one_time_step(
>
Ramprasad,
Let's first rule out some simple problems. Which version of deal.II are you
using? Did you compile deal.II with MPI and PETSc support?
We only define the namecpacePETScWrappers when DEAL_II_WITH_PETSC=ON. Can
you provide us with the content of the detailed.log file in your build
Hello All,
I am getting an error when I try to compile the step-18. It states
"PETScWrappers does not name a type". I wanted to know if anyone else has
faced the same problem. And if there is a solution to this.
Thank you.
--
The deal.II project is located at http://www.dealii.org/
For
Franco,
On Tuesday, July 16, 2019 at 3:54:34 AM UTC-4, Franco Milicchio wrote:
>
>
> Can I ask why there are some discrepancies between those two wrappers?
> Apart from an obvious one, that they're different libraries, but uniformity
> may help a lot IMO.
>
Because the wrappers have been
>
>
> step-40 tries to allow switching between Trilinos and PETSc as easy as
> possible by using the deal.II/lac/generic_linear_algebra.h header and the
> defining the LA::MPI namespace.
>
Thanks for your answer, Daniel. Even with that example I find it a little
harder than I expected to use
14 matches
Mail list logo