44   int Cnt_kind,Calc_CntOrbital_ON,spin,spinmax,m;  in DFT()  local
625     if      (SpinP_switch==0)  spinmax = 1;  in DFT()
626     else if (SpinP_switch==1)  spinmax = 2;  in DFT()
627     else if (SpinP_switch==3)  spinmax = 3;  in DFT()
637       ReRhok[m] = (double**)malloc(sizeof(double*)*spinmax);   in DFT()
638       for (spin=0; spin<spinmax; spin++){  in DFT()
646       ImRhok[m] = (double**)malloc(sizeof(double*)*spinmax);   in DFT()
647       for (spin=0; spin<spinmax; spin++){  in DFT()
653     Residual_ReRhok = (double**)malloc(sizeof(double*)*spinmax);   in DFT()
654     for (spin=0; spin<spinmax; spin++){  in DFT()
659     Residual_ImRhok = (double**)malloc(sizeof(double*)*spinmax);   in DFT()
660     for (spin=0; spin<spinmax; spin++){  in DFT()
742       PrintMemory("DFT: ReRhok",sizeof(double)*My_Max_NumGridB*List_YOUSO[38]*spinmax,NULL);  in DFT()
743       PrintMemory("DFT: ImRhok",sizeof(double)*My_Max_NumGridB*List_YOUSO[38]*spinmax,NULL);  in DFT()
744 …    PrintMemory("DFT: Residual_ReRhok", sizeof(double)*My_NumGridB_CB*List_YOUSO[38]*spinmax,NULL);  in DFT()
745 …    PrintMemory("DFT: Residual_ImRhok", sizeof(double)*My_NumGridB_CB*List_YOUSO[38]*spinmax,NULL);  in DFT()
2584     if      (SpinP_switch==0)  spinmax = 1;  in DFT()
2585     else if (SpinP_switch==1)  spinmax = 2;  in DFT()
2586     else if (SpinP_switch==3)  spinmax = 3;  in DFT()
2592       for (spin=0; spin<spinmax; spin++){  in DFT()
2600       for (spin=0; spin<spinmax; spin++){  in DFT()
2607     for (spin=0; spin<spinmax; spin++){  in DFT()
2612     for (spin=0; spin<spinmax; spin++){  in DFT()