Setter for Marko-Siggia bond.
U = [-(1/2)*rtilde^2 + 0.25/(1-rtilde)^2 + 0.25*rtilde]*(rmax/lp), r < rmax where rtilde = r/rmax.
User-set parameters:
lp
(persistence length)rmax
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=rp), | intent(inout), | dimension(:) | :: | params | ||
real(kind=rp), | intent(in), | optional | :: | lp | ||
real(kind=rp), | intent(in), | optional | :: | rmax |