I wrote a “fix_axial” in the spirit of “fix_uniaxial” but to independently control any of the 3 axes which is periodic (not necessarily conserve volume). The idea is that I want to control sigma_xx and epsilon_yy, so I would apply a fix_nph on x-axis to keep sigma_xx around zero and use my fix to expand/contract y-axis. (I live in 2D)
My “fix_axial” works perfectly fine by itself, but when I try to use it along with fix_nph to keep sigma_xx=0, the x-axis doesn’t seem to do any rescaling even when sigma_xx starts to increase.
Any gotchas anyone can think of? Looks like fix_nph should leave domain->boxylo, etc untouched and my fix_axial leaves domain->boxxlo, etc untouched.