Causal structure learning in high-dimensional systems remains a fundamental challenge in modern machine learning and statistical inference, particularly when underlying data-generating processes deviate from Gaussian assumptions. Linear Non-Gaussian Acyclic Models (LiNGAM) provide a principled framework for identifying causal directions using non-Gaussianity as an identification condition. However, scalability, computational efficiency, and reproducibility issues continue to limit their practical adoption in large-scale data environments. This study proposes a scalable Python-based architecture implemented through the PyCD-LiNGAM framework to address these limitations by integrating modular computation, optimized matrix operations, and automated causal graph discovery pipelines.

The proposed framework builds upon prior theoretical advancements in causal discovery and graphical modeling, particularly greedy structure learning strategies (Chickering, 2002), probabilistic graphical modeling principles (Drton & Maathuis, 2017), and linear non-Gaussian causal identification theory (Entner & Hoyer, 2011). Furthermore, it incorporates semiparametric inference perspectives for handling latent confounding structures (Bhattacharya et al., 2020). The system is evaluated conceptually for scalability, robustness to noise, and interpretability in non-Gaussian environments.

Results indicate that modular Python-based causal pipelines significantly enhance computational tractability while preserving theoretical identifiability guarantees under non-Gaussian assumptions. The study contributes a unified computational architecture bridging theoretical causal discovery models with practical implementation constraints, enabling reproducible and scalable causal inference workflows.