diff -Nru agate-202111141701+1.3.3/debian/changelog agate-202112111747+1.3.3/debian/changelog
--- agate-202111141701+1.3.3/debian/changelog	2021-11-14 17:28:33.000000000 +0000
+++ agate-202112111747+1.3.3/debian/changelog	2021-12-29 13:32:25.000000000 +0000
@@ -1,8 +1,8 @@
-agate (202111141701+1.3.3-2+421+7388d31~ubuntu21.10.1) impish; urgency=low
+agate (202112111747+1.3.3-2+422+a6874f5~ubuntu21.10.1) impish; urgency=low
 
   * Auto build.
 
- -- Didi Thunder <piti.diablotin@geduld.fr>  Sun, 14 Nov 2021 17:28:33 +0000
+ -- Didi Thunder <piti.diablotin@geduld.fr>  Wed, 29 Dec 2021 13:32:25 +0000
 
 agate (1.3.3-2) unstable; urgency=medium
 
diff -Nru agate-202111141701+1.3.3/debian/git-build-recipe.manifest agate-202112111747+1.3.3/debian/git-build-recipe.manifest
--- agate-202111141701+1.3.3/debian/git-build-recipe.manifest	2021-11-14 17:28:33.000000000 +0000
+++ agate-202112111747+1.3.3/debian/git-build-recipe.manifest	2021-12-29 13:32:25.000000000 +0000
@@ -1,3 +1,3 @@
-# git-build-recipe format 0.4 deb-version 202111141701+{debversion}+421+7388d31
-lp:~piti-diablotin/+git/agate git-commit:7388d310c26b3885d8b88d5e6af6465febb8db2e
+# git-build-recipe format 0.4 deb-version 202112111747+{debversion}+422+a6874f5
+lp:~piti-diablotin/+git/agate git-commit:a6874f58a0930ede32bac86811f5a616fb876934
 nest recipe lp:~piti-diablotin/+git/abiout_debian debian git-commit:84b3e1fa14702e77d51cda4e68d2ea074747899c
diff -Nru agate-202111141701+1.3.3/src/io/eigparserphbst.cpp agate-202112111747+1.3.3/src/io/eigparserphbst.cpp
--- agate-202111141701+1.3.3/src/io/eigparserphbst.cpp	2021-11-14 17:28:31.000000000 +0000
+++ agate-202112111747+1.3.3/src/io/eigparserphbst.cpp	2021-12-29 13:32:22.000000000 +0000
@@ -88,7 +88,7 @@
 
   geometry::vec3d prev_kpt = {{0.,0.,0.}};
   double length = 0.0;
-
+  double initLength =0.0;
   for ( unsigned ikpt = 0 ; ikpt < kpts ; ++ikpt ) {
     int status = nc_inq_varid(ncid, "qpoints", &varid);
     size_t start2[] = {ikpt,0};
@@ -122,7 +122,12 @@
 
     _eigenDisp.push_back(std::move(disp));
     if ( ikpt == 0 ) prev_kpt = kpt;
-    length += geometry::norm(geometry::operator-(kpt,prev_kpt));
+    if ( ikpt == 1 )
+       initLength = geometry::norm(geometry::operator-(kpt,prev_kpt));
+    if (geometry::norm(geometry::operator-(kpt,prev_kpt)) > 10*initLength)
+       length += 0;
+    else
+       length += geometry::norm(geometry::operator-(kpt,prev_kpt));
     _kpts.push_back(kpt);
     _lengths.push_back(length);
     prev_kpt = kpt;