op_randn_meat.hpp
Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031 template<typename eT>
00032 inline
00033 eT
00034 op_randn::randn()
00035 {
00036 const u32 N = 12;
00037 eT acc = eT(0);
00038
00039 for(u32 i=0; i<N; ++i)
00040 {
00041 acc += eT(std::rand()) / eT(RAND_MAX);
00042 }
00043
00044 return acc/eT(N) - eT(0.5);
00045 }
00046
00047
00048
00049 template<typename eT>
00050 inline
00051 void
00052 op_randn::direct_randn(eT* x, const u32 n_elem)
00053 {
00054 arma_extra_debug_sigprint();
00055
00056 for(u32 i=0; i<n_elem; ++i)
00057 {
00058 x[i] = eT(op_randn::randn<eT>());
00059 }
00060
00061 }
00062
00063
00064
00065 template<typename eT>
00066 inline
00067 void
00068 op_randn::apply(Mat<eT>& out, const Op<Mat<eT>,op_randn>& in)
00069 {
00070 arma_extra_debug_sigprint();
00071
00072 const u32 n_rows = in.aux_u32_a;
00073 const u32 n_cols = (in.aux_u32_b > 0) ? in.aux_u32_b : 1;
00074
00075 out.set_size(n_rows, n_cols);
00076
00077 op_randn::direct_randn(out.memptr(), out.n_elem);
00078 }
00079
00080
00081
00082 template<typename eT>
00083 inline
00084 void
00085 op_randn::apply(Mat<eT>& out, const Op<Col<eT>,op_randn>& in)
00086 {
00087 arma_extra_debug_sigprint();
00088
00089 out.set_size(in.aux_u32_a, 1);
00090
00091 op_randn::direct_randn(out.memptr(), out.n_elem);
00092 }
00093
00094
00095
00096 template<typename eT>
00097 inline
00098 void
00099 op_randn::apply(Mat<eT>& out, const Op<Row<eT>,op_randn>& in)
00100 {
00101 arma_extra_debug_sigprint();
00102
00103 out.set_size(1, in.aux_u32_a);
00104
00105 op_randn::direct_randn(out.memptr(), out.n_elem);
00106 }
00107
00108
00109
00110 template<typename eT>
00111 inline
00112 void
00113 op_randn::apply(Col<eT>& out, const Op<Col<eT>,op_randn>& in)
00114 {
00115 arma_extra_debug_sigprint();
00116
00117 out.set_size(in.aux_u32_a);
00118
00119 op_randn::direct_randn(out.memptr(), out.n_elem);
00120 }
00121
00122
00123
00124 template<typename eT>
00125 inline
00126 void
00127 op_randn::apply(Row<eT>& out, const Op<Row<eT>,op_randn>& in)
00128 {
00129 arma_extra_debug_sigprint();
00130
00131 out.set_size(in.aux_u32_a);
00132
00133 op_randn::direct_randn(out.memptr(), out.n_elem);
00134 }
00135
00136
00137
00138