forked from pytorch/pytorch
-
Notifications
You must be signed in to change notification settings - Fork 0
/
rmsprop_op_gpu.cu
51 lines (45 loc) · 1.15 KB
/
rmsprop_op_gpu.cu
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
#include "caffe2/sgd/rmsprop_op.h"
#include "caffe2/core/common_gpu.h"
#include "caffe2/core/context_gpu.h"
namespace caffe2 {
__global__ void RmsPropUpdate(
int N,
const float* g,
const float* ms,
const float* mom,
float* ng,
float* nms,
float* nmom,
float decay,
float momentum,
float epsilon,
const float* lr) {
CUDA_1D_KERNEL_LOOP(i, N) {
// Update new mean square estimate
nms[i] = ms[i] + (1.0f - decay) * (g[i] * g[i] - ms[i]);
// Update momentum estimate
nmom[i] =
mom[i] * momentum + lr[0] * g[i] / sqrtf(epsilon + nms[i]);
// New gradient is the momentum
ng[i] = nmom[i];
}
}
template <>
void rmsprop_update<CUDAContext>(
int N,
const float* g,
const float* ms,
const float* mom,
float* ng,
float* nms,
float* nmom,
float decay,
float momentum,
float epsilon,
const float* lr,
CUDAContext* context) {
RmsPropUpdate<<<CAFFE_GET_BLOCKS(N), CAFFE_CUDA_NUM_THREADS, 0, context->cuda_stream()>>>(
N, g, ms, mom, ng, nms, nmom, decay, momentum, epsilon, lr);
}
REGISTER_CUDA_OPERATOR(RmsProp, RmsPropOp<float, CUDAContext>);
}