Optimization of the gyroaverage operator based on hermite interpolation