<html><head><meta http-equiv="Content-Type" content="text/html charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div class=""><br class=""><div class=""><br class=""></div><div class="">Thank you,</div><div class="">Max<br class=""><div><blockquote type="cite" class=""><div class="">On 08 Aug 2016, at 17:23, Matthew Knepley <<a href="mailto:knepley@gmail.com" class="">knepley@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><div class="gmail_extra"><div class="gmail_quote">On Mon, Aug 8, 2016 at 10:15 AM, Maximilian Hartig<span class="Apple-converted-space"> </span><span dir="ltr" class=""><<a href="mailto:imilian.hartig@gmail.com" target="_blank" class="">imilian.hartig@gmail.com</a>></span><span class="Apple-converted-space"> </span>wrote:<br class=""><blockquote class="gmail_quote" style="margin: 0px 0px 0px 0.8ex; border-left-width: 1px; border-left-color: rgb(204, 204, 204); border-left-style: solid; padding-left: 1ex;"><div style="word-wrap: break-word;" class="">Thank you,<div class=""><br class=""></div><div class="">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 class="">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 class=""><br class=""></div><div class="">Could you send you example so I can run it myself?</div><div class=""><br class=""></div><div class=""> Thanks,</div><div class=""><br class=""></div><div class=""> Matt</div><div class=""> </div><blockquote class="gmail_quote" style="margin: 0px 0px 0px 0.8ex; border-left-width: 1px; border-left-color: rgb(204, 204, 204); border-left-style: solid; padding-left: 1ex;"><div style="word-wrap: break-word;" class=""><div class="">Thank you,</div><div class=""><br class=""></div><div class="">Max</div><div class=""><br class=""></div><span class=""><br class=""></span><span class="">[0]PETSC ERROR: ------------------------------<wbr class="">------------------------------<wbr class="">------------<br class="">[0]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range<br class="">[0]PETSC ERROR: Try option -start_in_debugger or -on_error_attach_debugger<br class="">[0]PETSC ERROR: or see<span class="Apple-converted-space"> </span><a href="http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind" target="_blank" class="">http://www.mcs.anl.gov/petsc/<wbr class="">documentation/faq.html#<wbr class="">valgrind</a><br class="">[0]PETSC ERROR: or try<span class="Apple-converted-space"> </span><a href="http://valgrind.org/" target="_blank" class="">http://valgrind.org</a><span class="Apple-converted-space"> </span>on GNU/linux and Apple Mac OS X to find memory corruption errors<br class="">[0]PETSC ERROR: likely location of problem given in stack below<br class="">[0]PETSC ERROR: --------------------- Stack Frames ------------------------------<wbr class="">------<br class="">[0]PETSC ERROR: Note: The EXACT line numbers in the stack are not available,<br class="">[0]PETSC ERROR: INSTEAD the line number of the start of the function<br class="">[0]PETSC ERROR: is given.<br class="">[0]PETSC ERROR: [0] PetscFEIntegrateResidual_Basic line 3503 /Users/maxhartig/PETSc/src/dm/<wbr class="">dt/interface/dtfe.c<br class="">[0]PETSC ERROR: [0] PetscFEIntegrateResidual line 5753 /Users/maxhartig/PETSc/src/dm/<wbr class="">dt/interface/dtfe.c<br class="">[0]PETSC ERROR: [0] DMPlexComputeResidual_Internal line 1706 /Users/maxhartig/PETSc/src/<wbr class="">snes/utils/dmplexsnes.c<br class="">[0]PETSC ERROR: [0] DMPlexSNESComputeResidualFEM line 2152 /Users/maxhartig/PETSc/src/<wbr class="">snes/utils/dmplexsnes.c<br class="">[0]PETSC ERROR: [0] SNESComputeFunction_DMLocal line 65 /Users/maxhartig/PETSc/src/<wbr class="">snes/utils/dmlocalsnes.c<br class="">[0]PETSC ERROR: [0] SNES user function line 2144 /Users/maxhartig/PETSc/src/<wbr class="">snes/interface/snes.c<br class="">[0]PETSC ERROR: [0] SNESComputeFunction line 2129 /Users/maxhartig/PETSc/src/<wbr class="">snes/interface/snes.c<br class="">[0]PETSC ERROR: [0] SNESSolve_NEWTONLS line 150 /Users/maxhartig/PETSc/src/<wbr class="">snes/impls/ls/ls.c<br class="">[0]PETSC ERROR: [0] SNESSolve line 3961 /Users/maxhartig/PETSc/src/<wbr class="">snes/interface/snes.c<br class="">[0]PETSC ERROR: [0] TS_SNESSolve line 188 /Users/maxhartig/PETSc/src/ts/<wbr class="">impls/implicit/theta/theta.c<br class="">[0]PETSC ERROR: [0] TSStep_Theta line 206 /Users/maxhartig/PETSc/src/ts/<wbr class="">impls/implicit/theta/theta.c<br class="">[0]PETSC ERROR: [0] TSStep line 3700 /Users/maxhartig/PETSc/src/ts/<wbr class="">interface/ts.c<br class="">[0]PETSC ERROR: [0] TSSolve line 3921 /Users/maxhartig/PETSc/src/ts/<wbr class="">interface/ts.c<br class="">[0]PETSC ERROR: --------------------- Error Message ------------------------------<wbr class="">------------------------------<wbr class="">--<br class="">[0]PETSC ERROR: Signal received<br class="">[0]PETSC ERROR: See<span class="Apple-converted-space"> </span><a href="http://www.mcs.anl.gov/petsc/documentation/faq.html" target="_blank" class="">http://www.mcs.anl.gov/petsc/<wbr class="">documentation/faq.html</a><span class="Apple-converted-space"> </span>for trouble shooting.<br class="">[0]PETSC ERROR: Petsc Release Version 3.7.2, unknown <br class="">[0]PETSC ERROR: Configure options --download-triangle<br class="">[0]PETSC ERROR: #1 User provided function() line 0 in unknown file<br class="">application called MPI_Abort(MPI_COMM_WORLD, 59) - process 0<br class="">[unset]: write_line error; fd=-1 buf=:cmd=abort exitcode=59<br class=""></span><span class=""><br class=""></span><div class=""><div class=""><blockquote type="cite" class=""><div class="">On 03 Aug 2016, at 16:44, Matthew Knepley <<a href="mailto:knepley@gmail.com" target="_blank" class="">knepley@gmail.com</a>> wrote:</div><br class=""><div class=""><div dir="ltr" class=""><div class="gmail_extra"><div class="gmail_quote">On Tue, Aug 2, 2016 at 8:22 AM, Maximilian Hartig<span class="Apple-converted-space"> </span><span dir="ltr" class=""><<a href="mailto:imilian.hartig@gmail.com" target="_blank" class="">imilian.hartig@gmail.com</a>></span><span class="Apple-converted-space"> </span>wrote:<br class=""><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 class=""><br class="">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 class="">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 class=""></blockquote><div class=""><br class=""></div><div class="">You can use</div><div class=""><br class=""></div><div class=""><div class=""> <span class="Apple-converted-space"> </span>ierr = DMTSSetBoundaryLocal(<wbr class="">adaptedDM, DMPlexTSComputeBoundary, user);CHKERRQ(ierr);</div><div class=""> <span class="Apple-converted-space"> </span>ierr = DMTSSetIFunctionLocal(<wbr class="">adaptedDM, DMPlexTSComputeIFunctionFEM, user);CHKERRQ(ierr);</div><div class=""> <span class="Apple-converted-space"> </span>ierr = DMTSSetIJacobianLocal(<wbr class="">adaptedDM, DMPlexTSComputeIJacobianFEM, user);CHKERRQ(ierr);</div></div><div class=""><br class=""></div><div class="">I have been meaning to write a heat equation example, but I have not finished yet,</div><div class=""><br class=""></div><div class=""> <span class="Apple-converted-space"> </span>Thanks,</div><div class=""><br class=""></div><div class=""> Matt</div><div class=""> </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 class="">Max</blockquote></div><br class=""><br clear="all" class=""><span class="HOEnZb"><font color="#888888" class=""><div class=""><br class=""></div>--<span class="Apple-converted-space"> </span><br class=""><div data-smartmail="gmail_signature" class="">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br class="">-- Norbert Wiener</div></font></span></div></div></div></blockquote></div><br class=""></div></div></blockquote></div><br class=""><br clear="all" class=""><div class=""><br class=""></div>--<span class="Apple-converted-space"> </span><br class=""><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 class="">-- Norbert Wiener</div></div></div></div></blockquote></div><br class=""></div></div></body></html>