Skip to content

fix(lammps): distinguish neutral charge and dipole columns#1020

Open
njzjz-bot wants to merge 1 commit into
deepmodeling:masterfrom
njzjz-bot:fix/issue-995
Open

fix(lammps): distinguish neutral charge and dipole columns#1020
njzjz-bot wants to merge 1 commit into
deepmodeling:masterfrom
njzjz-bot:fix/issue-995

Conversation

@njzjz-bot

Copy link
Copy Markdown
Contributor

Fixes #995.

Improve neutral decimal charge-style detection and prevent dipole fields from being reported as spins.

Tests: cd tests && python -m unittest test_lammps_atom_styles test_lammps_spin

Why existing tests missed it: Existing charge tests used non-integral charges, and no dipole-style fixture checked the optional spin key.

Coding agent: Codex
Codex version: codex-cli 0.144.4
Model: gpt-5.6-sol
Reasoning effort: xhigh

Improve automatic six-column detection for neutral decimal charges and restrict legacy spin extraction to atomic-style extension rows so dipole moments are not reported as spins. Add regressions for both cases; existing tests used non-neutral charges and therefore missed the ambiguity.

Coding-Agent: Codex
Codex-Version: codex-cli 0.144.4
Model: gpt-5.6-sol
Reasoning-Effort: xhigh
@dosubot dosubot Bot added size:XS This PR changes 0-9 lines, ignoring generated files. bug Something isn't working lammps labels Jul 16, 2026
@coderabbitai

coderabbitai Bot commented Jul 16, 2026

Copy link
Copy Markdown

Warning

Review limit reached

@njzjz-bot, you've reached your PR review limit, so we couldn't start this review.

Next review available in: 59 minutes

Enable usage-based reviews in Billing to review now. Otherwise, wait until the next included review is available.
You're only billed for reviews past your plan's rate limits ($0.25/file).

How can I continue?

After more reviews become available, a review can be triggered using the @coderabbitai review command as a PR comment. Alternatively, push new commits to this PR.

To avoid repeated limits, reduce automatic review volume by pausing incremental auto-reviews earlier, using label-based review opt-in, excluding WIP or generated PR titles, or requesting reviews manually when the PR is ready. If your team needs uninterrupted high-volume reviews, an organization admin can enable usage-based reviews.

How do review limits work?

CodeRabbit enforces per-developer PR review limits for each organization. Most developers receive the normal plan review availability.

For paid Pro and Pro+ PR reviews, CodeRabbit uses adaptive limits for sustained high-volume activity. When a developer's recent PR review activity reaches the 95th percentile or higher among CodeRabbit users, additional reviews become available more gradually as earlier reviews age out of the rolling window.

Please refer docs for additional details.

Review details
⚙️ Run configuration

Configuration used: Organization UI

Review profile: CHILL

Plan: Pro

Run ID: bb40616f-dd6d-4fb8-bbfe-8cde591bd8e4

📥 Commits

Reviewing files that changed from the base of the PR and between 0416b54 and 27d9e5e.

📒 Files selected for processing (2)
  • dpdata/formats/lammps/lmp.py
  • tests/test_lammps_atom_styles.py

Thanks for using CodeRabbit! It's free for OSS, and your support helps us grow. If you like it, consider giving us a shout-out.

❤️ Share

Comment @coderabbitai help to get the list of available commands.

@codspeed-hq

codspeed-hq Bot commented Jul 16, 2026

Copy link
Copy Markdown

Merging this PR will not alter performance

⚠️ Unknown Walltime execution environment detected

Using the Walltime instrument on standard Hosted Runners will lead to inconsistent data.

For the most accurate results, we recommend using CodSpeed Macro Runners: bare-metal machines fine-tuned for performance measurement consistency.

✅ 2 untouched benchmarks


Comparing njzjz-bot:fix/issue-995 (27d9e5e) with master (0416b54)

Open in CodSpeed

@codecov

codecov Bot commented Jul 16, 2026

Copy link
Copy Markdown

Codecov Report

✅ All modified and coverable lines are covered by tests.
✅ Project coverage is 86.88%. Comparing base (0416b54) to head (27d9e5e).
⚠️ Report is 1 commits behind head on master.

Additional details and impacted files
@@           Coverage Diff           @@
##           master    #1020   +/-   ##
=======================================
  Coverage   86.87%   86.88%           
=======================================
  Files          89       89           
  Lines        8266     8270    +4     
=======================================
+ Hits         7181     7185    +4     
  Misses       1085     1085           

☔ View full report in Codecov by Harness.
📢 Have feedback on the report? Share it here.

🚀 New features to boost your workflow:
  • ❄️ Test Analytics: Detect flaky tests, report on failures, and find test suite problems.

@wanghan-iapcm wanghan-iapcm left a comment

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Both #995 fixes look correct -- the ./e neutral-charge detection inspects the right column, the dipole guard correctly suppresses spurious spins, and the two new tests genuinely fail on pre-fix code. One follow-up on the spin whitelist, see inline.

def get_spins(lines: list[str], atom_style: str = "atomic") -> np.ndarray | None:
# Dipole/sphere/etc. styles use their extra columns for different physical
# quantities. Atomic and spin styles carry magnetic spin vectors.
if atom_style not in {"atomic", "spin"}:

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

"spin" here is unreachable: it isn't a key in ATOM_STYLE_COLUMNS, so detect_atom_style never returns it, and passing atom_style="spin" raises ValueError: Unsupported atom style: spin in _atom_info_style before get_spins runs. So only "atomic" actually enables spins, and the comment implies a "spin" style that doesn't exist.

Side effect: a genuine 9-column spin file not commented Atoms # atomic (e.g. # spin, or no comment) is auto-detected as dipole, so get_spins now returns None and its spins are silently dropped (they were returned before). dpdata's own round-trip output uses # atomic, so this only hits externally-produced spin files -- and those already read coordinates from the dipole column layout, so the ambiguity is pre-existing -- but the spin loss is new.

Consider either adding a real "spin" entry to ATOM_STYLE_COLUMNS (layout id type x y z spx spy spz sp) so the whitelist is honest and # spin files parse, or dropping "spin" from the set and comment to avoid implying unsupported behavior.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

bug Something isn't working lammps size:XS This PR changes 0-9 lines, ignoring generated files.

Projects

None yet

Development

Successfully merging this pull request may close these issues.

[Code scan] LAMMPS atom-style parsing misreads neutral charge and dipole data

2 participants