Skip to content

This issue was moved to a discussion.

You can continue the conversation there. Go to discussion →

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

mspms genome length limitation #1886

Closed
AlisaGU opened this issue Nov 3, 2021 · 5 comments
Closed

mspms genome length limitation #1886

AlisaGU opened this issue Nov 3, 2021 · 5 comments

Comments

@AlisaGU
Copy link

AlisaGU commented Nov 3, 2021

Hey, It's me again!
My friend used mspms to simulate a genome of 100M base pairs. After waiting a long time, the program is still running.
However, when he changed to python, program running quickly.
I am not sure what's the key.

@benjeffery
Copy link
Member

Could you give us some more details, ideally the code snippet and command line?

@AlisaGU
Copy link
Author

AlisaGU commented Nov 3, 2021

Sorry, I don't have these codes... Maybe my friend will give me a copy in tomorrow.
Empirically, will command line and python code run the same time ideally?

@grahamgower
Copy link
Member

Sorry, I don't have these codes... Maybe my friend will give me a copy in tomorrow. Empirically, will command line and python code run the same time ideally?

Well, msp will output a tree sequence, so this should run at a speed on par with using the python API to get a tree sequence (and then writing it to disk). Whereas mspms outputs the ms text format, which is likely to be slower. Without looking at the code/command in question, it's hard to guess if the text output is the cause of the slowness.

@grahamgower
Copy link
Member

Possibly related to #1652.

@AlisaGU
Copy link
Author

AlisaGU commented Nov 3, 2021

Thanks
Before this issue, I have read #1652
I simulate a genome of 10M using mspms and ms. Mspms just took 3 minutes. It's much faster than ms.
Now I am simulating 100M to observe what will happen.

@tskit-dev tskit-dev locked and limited conversation to collaborators Nov 4, 2021

This issue was moved to a discussion.

You can continue the conversation there. Go to discussion →

Labels
None yet
Projects
None yet
Development

No branches or pull requests

4 participants