<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Mon, Aug 8, 2016 at 10:15 AM, Maximilian Hartig <span dir="ltr"><<a href="mailto:imilian.hartig@gmail.com" target="_blank">imilian.hartig@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">Thank you,<div><br></div><div>I used these routines to setup a CN type TS. The problem I face now is that I seem to be unable to access the temporal derivatives u_t[..] for the definition of the residuals. I get a segmentation violation inside the PetscFEIntegrateResidual routine whenever I try to. I have attached the error message below.</div><div>Also, I am wondering whether it is possible to update the neumann boundary condition of the FE object for each timestep. This would be useful for coupling purposes.</div></div></blockquote><div><br></div><div>Could you send you example so I can run it myself?</div><div><br></div><div> Thanks,</div><div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><div>Thank you,</div><div><br></div><div>Max</div><div><br></div><span><br></span><span>[0]PETSC ERROR: ------------------------------<wbr>------------------------------<wbr>------------<br>[0]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range<br>[0]PETSC ERROR: Try option -start_in_debugger or -on_error_attach_debugger<br>[0]PETSC ERROR: or see <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind" target="_blank">http://www.mcs.anl.gov/petsc/<wbr>documentation/faq.html#<wbr>valgrind</a><br>[0]PETSC ERROR: or try <a href="http://valgrind.org" target="_blank">http://valgrind.org</a> on GNU/linux and Apple Mac OS X to find memory corruption errors<br>[0]PETSC ERROR: likely location of problem given in stack below<br>[0]PETSC ERROR: --------------------- Stack Frames ------------------------------<wbr>------<br>[0]PETSC ERROR: Note: The EXACT line numbers in the stack are not available,<br>[0]PETSC ERROR: INSTEAD the line number of the start of the function<br>[0]PETSC ERROR: is given.<br>[0]PETSC ERROR: [0] PetscFEIntegrateResidual_Basic line 3503 /Users/maxhartig/PETSc/src/dm/<wbr>dt/interface/dtfe.c<br>[0]PETSC ERROR: [0] PetscFEIntegrateResidual line 5753 /Users/maxhartig/PETSc/src/dm/<wbr>dt/interface/dtfe.c<br>[0]PETSC ERROR: [0] DMPlexComputeResidual_Internal line 1706 /Users/maxhartig/PETSc/src/<wbr>snes/utils/dmplexsnes.c<br>[0]PETSC ERROR: [0] DMPlexSNESComputeResidualFEM line 2152 /Users/maxhartig/PETSc/src/<wbr>snes/utils/dmplexsnes.c<br>[0]PETSC ERROR: [0] SNESComputeFunction_DMLocal line 65 /Users/maxhartig/PETSc/src/<wbr>snes/utils/dmlocalsnes.c<br>[0]PETSC ERROR: [0] SNES user function line 2144 /Users/maxhartig/PETSc/src/<wbr>snes/interface/snes.c<br>[0]PETSC ERROR: [0] SNESComputeFunction line 2129 /Users/maxhartig/PETSc/src/<wbr>snes/interface/snes.c<br>[0]PETSC ERROR: [0] SNESSolve_NEWTONLS line 150 /Users/maxhartig/PETSc/src/<wbr>snes/impls/ls/ls.c<br>[0]PETSC ERROR: [0] SNESSolve line 3961 /Users/maxhartig/PETSc/src/<wbr>snes/interface/snes.c<br>[0]PETSC ERROR: [0] TS_SNESSolve line 188 /Users/maxhartig/PETSc/src/ts/<wbr>impls/implicit/theta/theta.c<br>[0]PETSC ERROR: [0] TSStep_Theta line 206 /Users/maxhartig/PETSc/src/ts/<wbr>impls/implicit/theta/theta.c<br>[0]PETSC ERROR: [0] TSStep line 3700 /Users/maxhartig/PETSc/src/ts/<wbr>interface/ts.c<br>[0]PETSC ERROR: [0] TSSolve line 3921 /Users/maxhartig/PETSc/src/ts/<wbr>interface/ts.c<br>[0]PETSC ERROR: --------------------- Error Message ------------------------------<wbr>------------------------------<wbr>--<br>[0]PETSC ERROR: Signal received<br>[0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html" target="_blank">http://www.mcs.anl.gov/petsc/<wbr>documentation/faq.html</a> for trouble shooting.<br>[0]PETSC ERROR: Petsc Release Version 3.7.2, unknown <br>[0]PETSC ERROR: Configure options --download-triangle<br>[0]PETSC ERROR: #1 User provided function() line 0 in unknown file<br>application called MPI_Abort(MPI_COMM_WORLD, 59) - process 0<br>[unset]: write_line error; fd=-1 buf=:cmd=abort exitcode=59<br></span><span><br></span><div><div><blockquote type="cite"><div>On 03 Aug 2016, at 16:44, Matthew Knepley <<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>> wrote:</div><br><div><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Tue, Aug 2, 2016 at 8:22 AM, Maximilian Hartig <span dir="ltr"><<a href="mailto:imilian.hartig@gmail.com" target="_blank">imilian.hartig@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex">Hello all,<br>
<br>
I would like to run a transient problem with PetscFE. Example ex11.c seems relevant since it uses the PestcFV context to create boundary conditions and RHS Functions for the TS.<br>
Is there an easy way to do transient analysis with TS and petscFE or do I have to code my own time-stepping routine?<br></blockquote><div><br></div><div>You can use</div><div><br></div><div><div> ierr = DMTSSetBoundaryLocal(<wbr>adaptedDM, DMPlexTSComputeBoundary, user);CHKERRQ(ierr);</div><div> ierr = DMTSSetIFunctionLocal(<wbr>adaptedDM, DMPlexTSComputeIFunctionFEM, user);CHKERRQ(ierr);</div><div> ierr = DMTSSetIJacobianLocal(<wbr>adaptedDM, DMPlexTSComputeIJacobianFEM, user);CHKERRQ(ierr);</div></div><div><br></div><div>I have been meaning to write a heat equation example, but I have not finished yet,</div><div><br></div><div> Thanks,</div><div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex">
Thanks,<br>
Max</blockquote></div><br><br clear="all"><span class="HOEnZb"><font color="#888888"><div><br></div>-- <br><div data-smartmail="gmail_signature">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</font></span></div></div>
</div></blockquote></div><br></div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</div></div>