Hello,
You can find details about the implementation of temperature at this
paper. '
Adaptively time stepping the stochastic
Landau-Lifshitz-Gilbert equation at
nonzero temperature: Implementation and
validation in MuMax3' by Leliaert et al.
Yes, it should give temperature dependence of M, as well as the switching field. However, if you are wishing to simulate e.g. hysteresis curves, note that you need to use the run() command instead of minimize() or relax(). It works, but run is much slower than the energy minimization techniques.
For details, I would refer to the paper linked above. Mumax solves the stochastic LLG equation, rather than LLB. There are some limitations/issues (for example it tends to overestimate Curie temperature), that are covered in the paper.
Best,
Josh L.