(1) Input the noisy image corrupted by out-of-focus blur and data-dependent |

multiplicative Gamma noise. |

(2) Calculate the regularization weight as where is the noise |

standard deviation and (this was set empirically). |

(3) Specify the stopping threshold value . |

(4) Select a time step of for the PDE evolution. |

(5) Initialize the current restored image and |

(6) **while** (/) **do**, |

(7) Assign |

(8) Evaluate using the PDE defined in (10). |

(9) **end while** |