Error when repeating the tutorial for diffusion map in v1.9.1 scanpy

• I have checked that this issue has not already been reported.
• I have confirmed this bug exists on the latest version of scanpy.
• (optional) I have confirmed this bug exists on the master branch of scanpy.

Note: Please read this guide detailing how to provide the necessary information for us to reproduce your bug.

This is the code I'm using which was copied from the tutorial:

import scanpy as sc
adata = sc.datasets.paul15()
sc.pp.log1p(adata)
sc.pp.neighbors(adata, n_neighbors=20, use_rep='X', method='gauss')
sc.tl.diffmap(adata)
sc.tl.dpt(adata, n_branchings=1, n_dcs=10)
sc.pl.diffmap(adata, color=['dpt_pseudotime', 'dpt_groups', 'paul15_clusters'])

The results I got are:

Which is totally different than in the tutorial:

Can anyone run into this problem ever? Thanks for your help.

Versions

Start by rerunning the tutorial code with the reported Scanpy 1.9.1, AnnData 0.8.0, and Python 3.9.7 versions, focusing on sc.tl.diffmap(), sc.tl.dpt(), and sc.pl.diffmap(). Compare the generated plots with the tutorial images and determine whether the discrepancy is reproducible and attributable to the reported environment.