There are two definitions of Gorenstein singularities in the literature. Using Grothendieck's (or Serre's) duality, one defines the "dualizing sheaf" an object $\hat K_M$ of derived category of coherent sheaves such that $H^n(\hat K_M)=C$ and the multiplication $$ H^i(M, B) \times H^{n-i}(RHom^*(B, \hat K_M))\to H^n(\hat K_M)=C $$ defines a perfect pairing on the cohomology. Then $M$ is Gorenstein if $\hat K_M$ is a line bundle. This is the definition used in Hartshorne, Wikipedia, Stacks Project, Kollar-Mori and so on.
This is a stronger form of Cauhen-Macaulay property, which assumes that the complex of sheaves $\hat K_M$ is a coherent sheaf, that is, a complex of coherent sheaves quasi-isomorphic to a complex concentrated in one degree.
In their book "Toric varieties" Cox, Little and Schenck define a Gorenstein variety as the one with Cartier canonical bundle. The canonical bundle of a normal variety is defined as the reflexization of the sheaf of top-degree differential forms. In particular, this implies the Gorenstein property for all normal varieties with trivial canonical bundle outside of singularities. This definition is used in some other papers on toric geometry, and in papers on Sasakian geometry, for example Gauntlett, Martelli, Sparks, Yau "Obstructions to the Existence of Sasaki-Einstein Metrics".
Curiously enough, Cox, Little and Schenck use the usual definition of Cohen-Macaulay singularities (in another chapter) and never mention that Gorenstein singularities are a special case of Cohen-Macaulay. They also say (when defining the Gorenstein singularities) that there are other definitions around, and give a reference to Hartshorne. I suppose that they meant to say that their definition is equivalent to the standard, in particular, implies the Cohen-Macaulay property, but they never claim or prove it, as far as I can see.
I want a reference to the equivalence of these two notions. I think I have a vague idea how to prove it, but I am sure it is written up somewhere, maybe not in full generality (I need it only for isolated canonical singularities anyway).