664
664
if (((Respa *) update->integrate)->level_inner >= 0) respa = 1;
665
665
if (((Respa *) update->integrate)->level_middle >= 0) respa = 2;
667
if (respa == 0) irequest = neighbor->request(this);
667
if (respa == 0) irequest = neighbor->request(this,instance_me);
668
668
else if (respa == 1) {
669
irequest = neighbor->request(this);
669
irequest = neighbor->request(this,instance_me);
670
670
neighbor->requests[irequest]->id = 1;
671
671
neighbor->requests[irequest]->half = 0;
672
672
neighbor->requests[irequest]->respainner = 1;
673
irequest = neighbor->request(this);
673
irequest = neighbor->request(this,instance_me);
674
674
neighbor->requests[irequest]->id = 3;
675
675
neighbor->requests[irequest]->half = 0;
676
676
neighbor->requests[irequest]->respaouter = 1;
678
irequest = neighbor->request(this);
678
irequest = neighbor->request(this,instance_me);
679
679
neighbor->requests[irequest]->id = 1;
680
680
neighbor->requests[irequest]->half = 0;
681
681
neighbor->requests[irequest]->respainner = 1;
682
irequest = neighbor->request(this);
682
irequest = neighbor->request(this,instance_me);
683
683
neighbor->requests[irequest]->id = 2;
684
684
neighbor->requests[irequest]->half = 0;
685
685
neighbor->requests[irequest]->respamiddle = 1;
686
irequest = neighbor->request(this);
686
irequest = neighbor->request(this,instance_me);
687
687
neighbor->requests[irequest]->id = 3;
688
688
neighbor->requests[irequest]->half = 0;
689
689
neighbor->requests[irequest]->respaouter = 1;
692
} else irequest = neighbor->request(this);
692
} else irequest = neighbor->request(this,instance_me);
694
694
cut_coulsq = cut_coul * cut_coul;