new 'gradient' filter based on greens function from gradient elasticity theory. work in progress

This commit is contained in:
Pratheek Shanthraj 2013-06-21 17:26:05 +00:00
parent 2e8756b724
commit 7b2e48df68
1 changed files with 3 additions and 0 deletions

View File

@ -954,6 +954,9 @@ real(pReal) function utilities_getFilter(k)
utilities_getFilter = (1.0_pReal + cos(PI*k(3)/grid(3))) &
*(1.0_pReal + cos(PI*k(2)/grid(2))) &
*(1.0_pReal + cos(PI*k(1)/grid(1)))/8.0_pReal
case ('gradient') !< cosine curve with 1 for avg and zero for highest freq
utilities_getFilter = 1.0_pReal/(1.0_pReal + &
(k(1)*k(1) + k(2)*k(2) + k(3)*k(3)))
case default
call IO_error(error_ID = 892_pInt, ext_msg = trim(myfilter))
end select