Lines Matching refs:kin_m_aa
252 kin_mem->kin_m_aa = 0; in KINCreate()
798 if (kin_mem->kin_m_aa) { in KINAllocVectors()
801 … kin_mem->kin_R_aa = (realtype *) malloc((kin_mem->kin_m_aa*kin_mem->kin_m_aa) * sizeof(realtype)); in KINAllocVectors()
816 kin_mem->kin_gamma_aa = (realtype *) malloc(kin_mem->kin_m_aa * sizeof(realtype)); in KINAllocVectors()
832 kin_mem->kin_ipt_map = (long int *) malloc(kin_mem->kin_m_aa * sizeof(long int)); in KINAllocVectors()
849 kin_mem->kin_cv = (realtype *) malloc(2 * (kin_mem->kin_m_aa+1) * sizeof(realtype)); in KINAllocVectors()
867 kin_mem->kin_Xv = (N_Vector *) malloc(2 * (kin_mem->kin_m_aa+1) * sizeof(N_Vector)); in KINAllocVectors()
929 kin_mem->kin_df_aa = N_VCloneVectorArray((int) kin_mem->kin_m_aa,tmpl); in KINAllocVectors()
947 kin_mem->kin_liw += kin_mem->kin_m_aa * kin_mem->kin_liw1; in KINAllocVectors()
948 kin_mem->kin_lrw += kin_mem->kin_m_aa * kin_mem->kin_lrw1; in KINAllocVectors()
952 kin_mem->kin_dg_aa = N_VCloneVectorArray((int) kin_mem->kin_m_aa,tmpl); in KINAllocVectors()
966 N_VDestroyVectorArray(kin_mem->kin_df_aa, (int) kin_mem->kin_m_aa); in KINAllocVectors()
967 kin_mem->kin_liw -= (7 + kin_mem->kin_m_aa) * kin_mem->kin_liw1; in KINAllocVectors()
968 kin_mem->kin_lrw -= (7 + kin_mem->kin_m_aa) * kin_mem->kin_lrw1; in KINAllocVectors()
971 kin_mem->kin_liw += kin_mem->kin_m_aa * kin_mem->kin_liw1; in KINAllocVectors()
972 kin_mem->kin_lrw += kin_mem->kin_m_aa * kin_mem->kin_lrw1; in KINAllocVectors()
976 kin_mem->kin_q_aa = N_VCloneVectorArray((int) kin_mem->kin_m_aa,tmpl); in KINAllocVectors()
990 N_VDestroyVectorArray(kin_mem->kin_df_aa, (int) kin_mem->kin_m_aa); in KINAllocVectors()
991 N_VDestroyVectorArray(kin_mem->kin_dg_aa, (int) kin_mem->kin_m_aa); in KINAllocVectors()
992 kin_mem->kin_liw -= (7 + 2 * kin_mem->kin_m_aa) * kin_mem->kin_liw1; in KINAllocVectors()
993 kin_mem->kin_lrw -= (7 + 2 * kin_mem->kin_m_aa) * kin_mem->kin_lrw1; in KINAllocVectors()
996 kin_mem->kin_liw += kin_mem->kin_m_aa * kin_mem->kin_liw1; in KINAllocVectors()
997 kin_mem->kin_lrw += kin_mem->kin_m_aa * kin_mem->kin_lrw1; in KINAllocVectors()
1095 N_VDestroyVectorArray(kin_mem->kin_df_aa, (int) kin_mem->kin_m_aa); in KINFreeVectors()
1097 kin_mem->kin_lrw -= kin_mem->kin_m_aa * kin_mem->kin_lrw1; in KINFreeVectors()
1098 kin_mem->kin_liw -= kin_mem->kin_m_aa * kin_mem->kin_liw1; in KINFreeVectors()
1102 N_VDestroyVectorArray(kin_mem->kin_dg_aa, (int) kin_mem->kin_m_aa); in KINFreeVectors()
1104 kin_mem->kin_lrw -= kin_mem->kin_m_aa * kin_mem->kin_lrw1; in KINFreeVectors()
1105 kin_mem->kin_liw -= kin_mem->kin_m_aa * kin_mem->kin_liw1; in KINFreeVectors()
1109 N_VDestroyVectorArray(kin_mem->kin_q_aa, (int) kin_mem->kin_m_aa); in KINFreeVectors()
1111 kin_mem->kin_lrw -= kin_mem->kin_m_aa * kin_mem->kin_lrw1; in KINFreeVectors()
1112 kin_mem->kin_liw -= kin_mem->kin_m_aa * kin_mem->kin_liw1; in KINFreeVectors()
2336 if (kin_mem->kin_m_aa == 0) { in KINPicardAA()
2488 if (kin_mem->kin_m_aa == 0) { in KINFP()
2568 i_pt = iter-1 - ((iter-1) / kin_mem->kin_m_aa) * kin_mem->kin_m_aa; in AndersonAcc()
2596 } else if (iter <= kin_mem->kin_m_aa) { in AndersonAcc()
2602 R[(iter-1)*kin_mem->kin_m_aa+j] = N_VDotProd(kin_mem->kin_q_aa[j], kin_mem->kin_vtemp2); in AndersonAcc()
2603 …N_VLinearSum(ONE,kin_mem->kin_vtemp2, -R[(iter-1)*kin_mem->kin_m_aa+j], kin_mem->kin_q_aa[j], kin_… in AndersonAcc()
2605 …R[(iter-1)*kin_mem->kin_m_aa+iter-1] = SUNRsqrt(N_VDotProd(kin_mem->kin_vtemp2, kin_mem->kin_vtemp… in AndersonAcc()
2606 … N_VScale((1/R[(iter-1)*kin_mem->kin_m_aa+iter-1]), kin_mem->kin_vtemp2, kin_mem->kin_q_aa[i_pt]); in AndersonAcc()
2614 for (i=0; i < kin_mem->kin_m_aa-1; i++) { in AndersonAcc()
2615 a = R[(i+1)*kin_mem->kin_m_aa + i]; in AndersonAcc()
2616 b = R[(i+1)*kin_mem->kin_m_aa + i+1]; in AndersonAcc()
2620 R[(i+1)*kin_mem->kin_m_aa + i] = temp; in AndersonAcc()
2621 R[(i+1)*kin_mem->kin_m_aa + i+1] = ZERO; in AndersonAcc()
2623 if (i < kin_mem->kin_m_aa-1) { in AndersonAcc()
2624 for (j = i+2; j < kin_mem->kin_m_aa; j++) { in AndersonAcc()
2625 a = R[j*kin_mem->kin_m_aa + i]; in AndersonAcc()
2626 b = R[j*kin_mem->kin_m_aa + i+1]; in AndersonAcc()
2628 R[j*kin_mem->kin_m_aa + i+1] = -s*a + c*b; in AndersonAcc()
2629 R[j*kin_mem->kin_m_aa + i] = temp; in AndersonAcc()
2638 for (i = 1; i < kin_mem->kin_m_aa; i++) { in AndersonAcc()
2639 for (j = 0; j < kin_mem->kin_m_aa-1; j++) { in AndersonAcc()
2640 R[(i-1)*kin_mem->kin_m_aa + j] = R[i*kin_mem->kin_m_aa + j]; in AndersonAcc()
2646 for (j=0; j < (kin_mem->kin_m_aa-1); j++) { in AndersonAcc()
2647 …R[(kin_mem->kin_m_aa-1)*kin_mem->kin_m_aa+j] = N_VDotProd(kin_mem->kin_q_aa[j], kin_mem->kin_vtemp… in AndersonAcc()
2648 …N_VLinearSum(ONE, kin_mem->kin_vtemp2, -R[(kin_mem->kin_m_aa-1)*kin_mem->kin_m_aa+j], kin_mem->kin… in AndersonAcc()
2650 …R[(kin_mem->kin_m_aa-1)*kin_mem->kin_m_aa+kin_mem->kin_m_aa-1] = SUNRsqrt(N_VDotProd(kin_mem->kin_… in AndersonAcc()
2651 …VScale((1/R[(kin_mem->kin_m_aa-1)*kin_mem->kin_m_aa+kin_mem->kin_m_aa-1]), kin_mem->kin_vtemp2, ki… in AndersonAcc()
2655 for (i=i_pt+1; i < kin_mem->kin_m_aa; i++) in AndersonAcc()
2663 if (kin_mem->kin_m_aa < iter) lAA = kin_mem->kin_m_aa; in AndersonAcc()
2675 gamma[i] = gamma[i]-R[j*kin_mem->kin_m_aa+i]*gamma[j]; in AndersonAcc()
2677 gamma[i] = gamma[i]/R[i*kin_mem->kin_m_aa+i]; in AndersonAcc()