Praveen,

>From what I can discern, it seems that you've posted this twice:
https://groups.google.com/forum/#!topic/dealii/4TRWnmMywRI

Do you have any objections to me deleting this post?
J-P

On Friday, June 24, 2016 at 8:37:45 AM UTC+2, Praveen C wrote:
>
> Dear all
>
> I am using the new TrilinosWrappers::SolverDirect which allows two step 
> solution: initialize and solve. Since I want to calculate LU decomposition 
> only once, I do something like this
>
>       static int first_time = 1;
>
>       
>       static TrilinosWrappers::SolverDirect::AdditionalData data (false, 
> "Amesos_Mumps");
>       static SolverControl solver_control (1, 0);
>       static TrilinosWrappers::SolverDirect direct (solver_control, data);
>
>       
>       // If it is first time, compute LU decomposition
>       if(first_time)
>       {
>          pcout << "Performing LU decomposition\n";
>          direct.initialize (mass_matrix);
>          first_time = 0;
>       }
>
>       
>       // solve for ax
>       {
>          TrilinosWrappers::MPI::Vector tmp (locally_owned_dofs, 
> mpi_communicator);
>          direct.solve (tmp, rhs_ax);
>          constraints.distribute (tmp);
>          ax = tmp;
>       }
>
> This works but then program finishes, I get this message.
>
>
> *** The MPI_Comm_f2c() function was called after MPI_FINALIZE was invoked.
> *** This is disallowed by the MPI standard.
> *** Your MPI job will now abort.
> [(null):4853] Local abort after MPI_FINALIZE completed successfully; not 
> able to aggregate error messages, and not able to guarantee that all other 
> processes were killed!
> *** The MPI_Comm_f2c() function was called after MPI_FINALIZE was invoked.
> *** This is disallowed by the MPI standard.
> *** Your MPI job will now abort.
> [(null):4854] Local abort after MPI_FINALIZE completed successfully; not 
> able to aggregate error messages, and not able to guarantee that all other 
> processes were killed!
> -------------------------------------------------------
> Primary job  terminated normally, but 1 process returned
> a non-zero exit code.. Per user-direction, the job has been aborted.
> -------------------------------------------------------
> --------------------------------------------------------------------------
> mpirun detected that one or more processes exited with non-zero status, 
> thus causing
> the job to be terminated. The first process to do so was:
>
>   Process name: [[14617,1],0]
>   Exit code:    1
> —————————————————————————————————————
>
> Removing “static” for the solver and recomputing it everytime gets rid of 
> this error. So I am guessing it is the “static” qualifier that is causing 
> problem. How can I fix this issue ?
>
> Thanks
> praveen
>

-- 
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 received this message because you are subscribed to the Google Groups 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

Reply via email to