Implementation MMA Code (Svanberg) in 88 line code (Andreassen, Sigmund) with Heaviside projection method
Dear members,
I'm strugling with implementing the MMA code (Svanberg) in the 88 line code with Heaviside projection method.
The sensitivity and density filter both work correctly, but when ft == 3 (Heaviside) my volfrac goes to 1 (when it's supposed to be 0.5 or something else)
Can somebody help me with this?
I just paste the code here (without the MMA part, mmasub & subsolv, because its too long and I have not changed anything different than the original of Svanberg)