Re: [petsc-users] Periodic BC in petsc DMPlex / firedrake

_type dirichlet > -petscspace_degree 1 -vec_view vtk:test.vtu:vtk_vtu -f > ${PETSC_DIR}/share/petsc/datafiles/meshes/square_periodic.msh > -dm_plex_gmsh_periodic > > Long time ago I added support for viewing meshes with periodic vertices in > the VTK_VTU viewer, but I did not fix the

Re: [petsc-users] Periodic BC in petsc DMPlex / firedrake

ou should definitely subscribe! > > > On 24 Oct 2018, at 09:27, Maximilian Hartig > <mailto:imilian.har...@gmail.com>> wrote: > > > > Hello Lawrence, > > > > sorry to message you out of the blue. My name is Max and I found your post > > on GitHub (h

Re: [petsc-users] Hints for using petscfe for plasticity -- how to update/access internal variables?

> On 21. Sep 2017, at 17:33, Matthew Knepley <knep...@gmail.com> wrote: > > On Thu, Sep 21, 2017 at 11:28 AM, Maximilian Hartig <imilian.har...@gmail.com > <mailto:imilian.har...@gmail.com>> wrote: >> On 20. Sep 2017, at 22:51, Matthew Knepley <knep...@

Re: [petsc-users] Hints for using petscfe for plasticity -- how to update/access internal variables?

> On 21. Sep 2017, at 16:02, Mark Adams wrote: > >>> > > That depends on the type of plasticity. For a linear hardening formulation it > is correct that I could just apply a local operator. I’d be happy with that > for now. > > I would just do this for now and get it

Re: [petsc-users] Hints for using petscfe for plasticity -- how to update/access internal variables?

> On 20. Sep 2017, at 22:51, Matthew Knepley <knep...@gmail.com> wrote: > > On Wed, Sep 20, 2017 at 3:46 PM, Maximilian Hartig <imilian.har...@gmail.com > <mailto:imilian.har...@gmail.com>> wrote: > >> On 20. Sep 2017, at 19:05, Matthew Knepley <knep.

Re: [petsc-users] Hints for using petscfe for plasticity -- how to update/access internal variables?

> On 20. Sep 2017, at 19:05, Matthew Knepley <knep...@gmail.com> wrote: > > On Wed, Sep 20, 2017 at 12:57 PM, Maximilian Hartig <imilian.har...@gmail.com > <mailto:imilian.har...@gmail.com>> wrote: >> On 20. Sep 2017, at 18:17, Matthew Knepley <knep...@

Re: [petsc-users] Hints for using petscfe for plasticity -- how to update/access internal variables?

> On 20. Sep 2017, at 18:17, Matthew Knepley <knep...@gmail.com> wrote: > > On Wed, Sep 20, 2017 at 11:46 AM, Maximilian Hartig <imilian.har...@gmail.com > <mailto:imilian.har...@gmail.com>> wrote: > Hello, > > I’m trying to implement plasticity u

[petsc-users] Hints for using petscfe for plasticity -- how to update/access internal variables?

Hello, I’m trying to implement plasticity using petscFE but I am quite stuck since a while. Here’s what I’m trying to do: I have a TS which solves the following equation: gradient(stress) +Forces = density*acceleration where at the moment stress is a linear function of the strain and hence the

Re: [petsc-users] Gradient in DMProjectField

On 12. Sep 2017, at 21:26, Matthew Knepley <knep...@gmail.com> wrote:On Tue, Sep 12, 2017 at 12:52 PM, Maximilian Hartig <imilian.har...@gmail.com> wrote:Hello, in my attempt to use an incremental formulation for elastoplasticity I wanted to create an auxiliary field in which the s

[petsc-users] Gradient in DMProjectField

Hello, in my attempt to use an incremental formulation for elastoplasticity I wanted to create an auxiliary field in which the stresses from the past tilmestep are stored. I managed to create an auxiliary DM on which I have a field for the displacement, one for the stress and one for the inner

Re: [petsc-users] accessing fields from previous, converged solution.

> On 28. Aug 2017, at 12:31, Matthew Knepley <knep...@gmail.com> wrote: > > On Mon, Aug 28, 2017 at 5:15 AM, Maximilian Hartig <imilian.har...@gmail.com > <mailto:imilian.har...@gmail.com>> wrote: > Sorry to bother you again, but I have checked and rechecked my

Re: [petsc-users] accessing fields from previous, converged solution.

it with auxiliary fields that get updated every timestep but I gather that is not a good idea. Could you please point me in some direction? I appreciate you patience with my questioning. Thanks, Max > On 22. Aug 2017, at 13:52, Maximilian Hartig <imilian.har...@gmail.com> wrote: > >>

Re: [petsc-users] accessing fields from previous, converged solution.

> On 17. Aug 2017, at 13:51, Matthew Knepley <knep...@gmail.com> wrote: > > On Thu, Aug 17, 2017 at 3:13 AM, Maximilian Hartig <imilian.har...@gmail.com > <mailto:imilian.har...@gmail.com>> wrote: > Thanks for your help Matt. >> On 16. Aug 2017, at 16

Re: [petsc-users] accessing fields from previous, converged solution.

Thanks for your help Matt. > On 16. Aug 2017, at 16:17, Matthew Knepley <knep...@gmail.com> wrote: > > On Wed, Aug 16, 2017 at 8:56 AM, Maximilian Hartig <imilian.har...@gmail.com > <mailto:imilian.har...@gmail.com>> wrote: > Hello, > > I have a p

[petsc-users] accessing fields from previous, converged solution.

Hello, I have a problem with several fields that I solve with PetscFE and TS. I now need to access the solution from the previous timestep to compute the residual for the current timestep. I tried a TSMonitor with the following code in it: TSGetDM(ts,); DMClone(dm,_aux); DMGetDS(dm,_aux);

Re: [petsc-users] -snes_fd for problems with residuals with non-continuous first derivative?

an Auxiliary field. Thanks, Max 2017-07-05 18:46 GMT+02:00 Matthew Knepley <knep...@gmail.com>: > On Wed, Jul 5, 2017 at 9:39 AM, Maximilian Hartig < > imilian.har...@gmail.com> wrote: > >> I do not clearly understand the discrimination between local and global >> pl

Re: [petsc-users] -snes_fd for problems with residuals with non-continuous first derivative?

> wrote: > > Let me suggest that you grab a hold of Simo and Hughes, Computational > Inelasticity, Springer-Verlag (1998). It explains a lot about how to set > up this problem -- in particular Chapter 1 gives a comprehensive > one-dimensional tutorial on everything you need to kn

Re: [petsc-users] -snes_fd for problems with residuals with non-continuous first derivative?

he solution process. > > > > > > On Jul 5, 2017, at 11:24 AM, Maximilian Hartig <imilian.har...@gmail.com> > wrote: > > > > The solver goes on for some timesteps, then it gets hung up in step 10. > The output is as follows: > > > > 28 SNES Function norm

Re: [petsc-users] -snes_fd for problems with residuals with non-continuous first derivative?

R H J Peerlings > > Best regards, > Luv > > On 30 Jun 2017, at 11:52, Maximilian Hartig <imilian.har...@gmail.com> > wrote: > > Hi Luv, > > I’m modelling linear hardening(sigma_y = sigma_y0 + > K_iso*epsilon_plast_eqiv) with isotropic plasticity only. So I s

Re: [petsc-users] -snes_fd for problems with residuals with non-continuous first derivative?

sidual -pc_type lu > -ksp_converged_reason -snes_converged_reason > > > > On Jun 30, 2017, at 3:45 AM, Maximilian Hartig <imilian.har...@gmail.com> > wrote: > > > > Hello, > > > > I’m trying to implement plasticity and have problems getting the Petsc > SNES to conve

Re: [petsc-users] -snes_fd for problems with residuals with non-continuous first derivative?

ons of plasticity but I cannot seem to find my mistake. Thanks, Max > On 30. Jun 2017, at 11:09, Luv Sharma <luvsharm...@gmail.com> wrote: > > Hi Max, > > Is your field solver not converging or the material point solver ;)? > > Best regards, > Luv >> On 3

Re: [petsc-users] Problems imposing boundary conditions

the Jacobian and the iFunction? Thanks, Max > On 7 Mar 2017, at 18:21, Matthew Knepley <knep...@gmail.com> wrote: > > On Tue, Mar 7, 2017 at 11:11 AM, Maximilian Hartig <imilian.har...@gmail.com > <mailto:imilian.har...@gmail.com>> wrote: > >> On 7 Ma

Re: [petsc-users] Problems imposing boundary conditions

> On 7 Mar 2017, at 16:29, Matthew Knepley <knep...@gmail.com> wrote: > > On Tue, Mar 7, 2017 at 3:28 AM, Maximilian Hartig <imilian.har...@gmail.com > <mailto:imilian.har...@gmail.com>> wrote: > It seems you are correct. In theory, the problem should not

Re: [petsc-users] Problems imposing boundary conditions

t; > On Mon, Mar 6, 2017 at 8:38 AM, Maximilian Hartig <imilian.har...@gmail.com > <mailto:imilian.har...@gmail.com>> wrote: > Of course, please find the source as well as the mesh attached below. I run > with: > > -def_petscspace_order 2 -vel_petscspace_order

Re: [petsc-users] Problems imposing boundary conditions

.@gmail.com > <mailto:lukas.drinkt.t...@gmail.com>> wrote: > You have apparently preallocated the non-zeroes of you matrix, and the room > was insufficient to accommodate all your equations. > > What happened after you tried: > > MatSetOption(A, MAT_NEW_NON

[petsc-users] Problems imposing boundary conditions

Hello, I am working on a transient structural FEM code with PETSc. I managed to create a slow but functioning program with the use of petscFE and a TS solver. The code runs fine until I try to restrict movement in all three spatial directions for one face. I then get the error which is

Re: [petsc-users] TS and petscFE -unable to access u_t

Gladly. I also send the simple gmsh I am using. cimply.c Description: Binary data testmesh_2D_box_quad.msh Description: Binary data Thank you,MaxOn 08 Aug 2016, at 17:23, Matthew Knepley <knep...@gmail.com> wrote:On Mon, Aug 8, 2016 at 10:15 AM, Maximilian Hartig <imilian.har...@

Re: [petsc-users] TS and petscFE -unable to access u_t

provided function() line 0 in unknown file application called MPI_Abort(MPI_COMM_WORLD, 59) - process 0 [unset]: write_line error; fd=-1 buf=:cmd=abort exitcode=59 > On 03 Aug 2016, at 16:44, Matthew Knepley <knep...@gmail.com> wrote: > > On Tue, Aug 2, 2016 at 8:22 AM, Ma

[petsc-users] TS and petscFE

Hello all, 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. Is there an easy way to do transient analysis with TS and petscFE or do I have to code my own