Lose memory without lost contact

Hi everyone,

I have written a contact law using memory (history_offset = hsetup->add_history_value("memory", "1");), as in tangential_model_history.h. I observe memory was reset whereas contact was not lost.

I write a minimal normal law and example which lose memory. Law is like hooke/stiffness, with a dozen of line modified. I add a memory, which make nothing, not used, just printed, and initialized to 1 at first time. When memory is lost, this parameter is reinitialized to 0 (because it is no more the first time). The output file is what I see in my terminal (printing are from "memory is 1" for each contact to "memory is 0" when losing memory).

Note bigger neighbor is, more conserved is the memory. And less is timestep, more conserved is the memory. I have verify, contact is *not* lost in this example (but memory is), particles are indented during all simulation.

Does someone have a solution, or the same problem?

Thanks!
Étienne.