fixed solver based on Richardson iteration