Felix,
the usual approach for verifying your code is the method of manufactured
solution. Choose functions that are elements of your ansatz space for all
solution variables
(possibly chaningyour grid to a cartesian one) and compute the
corresponding right-hand side. Check that the error
> Would that work together with the GMG-preconditioner?
You mean you want to use the linear operator as a smoother for the GMG? Or
where in the GMG do you want to use it?
> Will test it there as
> soon as my problem in my second question could be solved.
> There I have the following code:
> |
Would that work together with the GMG-preconditioner? Will test it there as
soon as my problem in my second question could be solved.
There I have the following code:
SolverControlcoarse_solver_control(1000, 1e-10, false, false);
SolverGMRES
On 7/19/19 9:59 AM, 'Maxi Miller' via deal.II User Group wrote:
> What do the parallel multigrid functions do? Could you just have your if()
> inside that function -- i.e., it is called on all cells, but does nothing
> if
> the cell doesn't allow you to do anything?
>
>
> That's
>
> What do the parallel multigrid functions do? Could you just have your if()
> inside that function -- i.e., it is called on all cells, but does nothing
> if
> the cell doesn't allow you to do anything?
>
That's exactly what my if()-conditions should do (but is not doing)...
Am Freitag,
On 7/19/19 9:49 AM, 'Maxi Miller' via deal.II User Group wrote:
> I can not use that function, else I will get
> An error occurred in line <3543> of file
> in function
> bool dealii::CellAccessor::is_locally_owned() const [with
> int dim = 2; int spacedim = 2]
> The violated condition was:
On 7/19/19 9:47 AM, 'Maxi Miller' via deal.II User Group wrote:
> What I did was the replacement of system_matrix with a function, providing a
> vmult()-function. Now I would like to add a preconditioner based on that
> function to reduce the necessary GMRES-iterations. But until now most of the
On 7/18/19 8:22 AM, Pawan Kumar wrote:
>
> I have recently updated my Ubuntu, deal.II and eclipse, I am getting some
> warnings (red colored underlined as in attached screenshot). The eclipse
> project was made using cmake and the paths are also seem correct(another
> screenshot attached).
I can not use that function, else I will get
An error occurred in line <3543> of file
in function
bool dealii::CellAccessor::is_locally_owned() const
[with int dim = 2; int spacedim = 2]
The violated condition was:
this->active()
Additional information:
is_locally_owned() can
What I did was the replacement of system_matrix with a function, providing
a vmult()-function. Now I would like to add a preconditioner based on that
function to reduce the necessary GMRES-iterations. But until now most of
the preconditioners require a matrix for initialization. Thus I wanted
On 7/19/19 9:34 AM, 'Maxi Miller' via deal.II User Group wrote:
> I.e. (based on step-37) I can write something like
> |
> LinearOperator jacobian_operator;
> //vmult-operation is added later
>
> using SystemMatrixType = jacobian_operator;
> SystemMatrixType system_matrix;
> |
>
> and follow the
On 7/19/19 9:31 AM, 'Maxi Miller' via deal.II User Group wrote:
> Based on the stacktrace the problematic access happens within
> mg_cell_worker():
> ...
> which is guarded by the if-condition.
Then we have a miracle :-) You'll have to check in a debugger how you got into
this place. If you
I.e. (based on step-37) I can write something like
LinearOperator jacobian_operator;
//vmult-operation is added later
using SystemMatrixType = jacobian_operator;
SystemMatrixType system_matrix;
and follow the example accordingly?
Am Freitag, 19. Juli 2019 17:26:31 UTC+2 schrieb Wolfgang
Muhammad,
> Hi! I am successful in running the code and it
> is without a doubt a nice addition to deal.ii code gallery. It reproduces the
> expected result for 3d Beam case now.
Great. Can you say what you had to change to make things run? We should apply
Based on the stacktrace the problematic access happens within
mg_cell_worker():
Stacktrace:
---
#0 /opt/dealii/lib/libdeal_II.g.so.9.2.0-pre: void
dealii::DoFCellAccessor,
true>::get_dof_values(dealii::TrilinosWrappers::MPI::Vector const&, double*, double*)
const
#1
On 7/19/19 9:18 AM, 'Maxi Miller' via deal.II User Group wrote:
> Is it then possible to use this LinearOperator in a preconditioner? The code
> works without one (i.e. PreconditionIdentity()), but the Iterations are
> increasing quite fast, thus negating all speedup gained from the
>
On 7/19/19 9:15 AM, 'Maxi Miller' via deal.II User Group wrote:
> I am trying to port example 15 to a MPI-supported version, using the
> geometric
> multigrid method as preconditioner (to get more familiar with the
> implementation). I followed example 50 while doing that, but currently I get
Is it then possible to use this LinearOperator in a preconditioner? The
code works without one (i.e. PreconditionIdentity()), but the Iterations
are increasing quite fast, thus negating all speedup gained from the
LinearOperator.
I tried to find something in the examples, but only found data
I am trying to port example 15 to a MPI-supported version, using the
geometric multigrid method as preconditioner (to get more familiar with the
implementation). I followed example 50 while doing that, but currently I
get the error
An error occurred in line <3866> of file
Franco
Le ven. 19 juil. 2019 à 10:04, Franco Milicchio
a écrit :
> Is there a small documentation that I can read? For instance, I see there are
> several namespaces like LinearAlgebra, PETSc/Epetra/Trilinos/CUDAWrappers,
> LinearAlgebraDealii and so on, that might be restructured without
>
> > Thanks for your answer. By the way, would a wrapper for the Eigen
> library be appealing to anyone?
> I am sure it would be useful!
>
>
Good to know, Bruno. I will start looking into those wrappers and I hope
some student of mine will be helping (fingers crossed).
Is there a small
> > Thanks for your answer. By the way, would a wrapper for the Eigen
> library be appealing to anyone?
> I am sure it would be useful!
>
>
Good to know, Bruno. I will start looking into those wrappers and I hope
some student of mine will be helping (fingers crossed).
Is there a small
Hi Daniel,
The problem is now solved. The issue was that, the bash rc did not have the
location of PETSc. This line in the bash file solved the problem "export
PETSC_DIR=/path/to/petsc-x-y-z""
Regards,
Ramprasad
On Wed, Jul 17, 2019 at 5:31 PM Bruno Turcksin
wrote:
> Ramprasad,
>
> Le mer. 17
23 matches
Mail list logo