Dear Users, I am doing constrained relaxation using harmonic springs.
I have been able to successfully modify constr.f. However, I currently have to set the value of each constraint distance D in constr.f hand, recompile siesta and run the relaxation. I was wondering if I can define a block in the input FDF file and simply set the value D. That way I only have to compile siesta once and use different values of D in the FDF input. How may I do that? Which routines do I have to modify and what special declarations should I include in constr.f? Thank you. Regards, Francisco
-- SIESTA is supported by the Spanish Research Agency (AEI) and by the European H2020 MaX Centre of Excellence (http://www.max-centre.eu/)