Lines Matching refs:subset_size
3388 PetscInt subset_size; in PCBDDCAdaptiveSelection() local
3390 PetscCall(ISGetLocalSize(sub_schurs->is_subs[i], &subset_size)); in PCBDDCAdaptiveSelection()
3391 mss = PetscMax(mss, subset_size); in PCBDDCAdaptiveSelection()
3407 PetscInt n, subset_size; in PCBDDCAdaptiveSelection() local
3409 PetscCall(ISGetLocalSize(sub_schurs->is_subs[i], &subset_size)); in PCBDDCAdaptiveSelection()
3410 n = PetscMin(subset_size, nmax); in PCBDDCAdaptiveSelection()
3411 cum += subset_size; in PCBDDCAdaptiveSelection()
3412 cum2 += subset_size * n; in PCBDDCAdaptiveSelection()
3495 PetscInt j, subset_size, eigs_start = 0; in PCBDDCAdaptiveSelection() local
3513 PetscCall(ISGetLocalSize(sub_schurs->is_subs[i], &subset_size)); in PCBDDCAdaptiveSelection()
3515 PetscCall(PetscBLASIntCast(subset_size, &B_N)); in PCBDDCAdaptiveSelection()
3521 for (j = 0; j < subset_size; j++) { in PCBDDCAdaptiveSelection()
3522 if (PetscRealPart(*(Sarray + cumarray + j * (subset_size + 1))) < 0.0) { in PCBDDCAdaptiveSelection()
3523 … PetscCall(MatCreateSeqDense(PETSC_COMM_SELF, subset_size, subset_size, Sarray + cumarray, &T)); in PCBDDCAdaptiveSelection()
3526 … PetscCall(MatCreateSeqDense(PETSC_COMM_SELF, subset_size, subset_size, Starray + cumarray, &T)); in PCBDDCAdaptiveSelection()
3536 *(Sarray + cumarray + idxs[k] * (subset_size + 1)) *= -1.0; in PCBDDCAdaptiveSelection()
3537 *(Starray + cumarray + idxs[k] * (subset_size + 1)) = 0.0; in PCBDDCAdaptiveSelection()
3551 PetscCall(PetscArrayzero(S, subset_size * subset_size)); in PCBDDCAdaptiveSelection()
3552 PetscCall(PetscArrayzero(St, subset_size * subset_size)); in PCBDDCAdaptiveSelection()
3554 for (j = 0; j < subset_size; j++) { in PCBDDCAdaptiveSelection()
3555 for (k = j; k < subset_size; k++) { in PCBDDCAdaptiveSelection()
3556 S[j * subset_size + k] = Sarray[cumarray + j * subset_size + k]; in PCBDDCAdaptiveSelection()
3557 St[j * subset_size + k] = Starray[cumarray + j * subset_size + k]; in PCBDDCAdaptiveSelection()
3561 PetscCall(PetscArraycpy(S, Sarray + cumarray, subset_size * subset_size)); in PCBDDCAdaptiveSelection()
3562 PetscCall(PetscArraycpy(St, Starray + cumarray, subset_size * subset_size)); in PCBDDCAdaptiveSelection()
3569 …1 && pcbddc->use_deluxe_scaling) PetscCall(PetscArraycmp(S, St, subset_size * subset_size, &same_d… in PCBDDCAdaptiveSelection()
3582 if (nmin >= subset_size) compute_range = PETSC_FALSE; in PCBDDCAdaptiveSelection()
3589 …MT " (range %d) (change %" PetscInt_FMT ").\n", i, sub_schurs->n_subs, subset_size, c, w, compute_… in PCBDDCAdaptiveSelection()
3601 PetscCall(PetscLogFlops((4.0 * subset_size * subset_size * subset_size) / 3.0)); in PCBDDCAdaptiveSelection()
3621 PetscCall(PetscLogFlops((4.0 * subset_size * subset_size * subset_size) / 3.0)); in PCBDDCAdaptiveSelection()
3631 PetscCall(PetscLogFlops((4.0 * subset_size * subset_size * subset_size) / 3.0)); in PCBDDCAdaptiveSelection()
3637 PetscCall(PetscArraycpy(S, Sarray + cumarray, subset_size * subset_size)); in PCBDDCAdaptiveSelection()
3638 PetscCall(PetscArraycpy(St, Starray + cumarray, subset_size * subset_size)); in PCBDDCAdaptiveSelection()
3644 PetscCall(PetscLogFlops((4.0 * subset_size * subset_size * subset_size) / 3.0)); in PCBDDCAdaptiveSelection()
3657 PetscCall(PetscLogFlops((4.0 * subset_size * subset_size * subset_size) / 3.0)); in PCBDDCAdaptiveSelection()
3673 PetscCall(PetscLogFlops((4.0 * subset_size * subset_size * subset_size) / 3.0)); in PCBDDCAdaptiveSelection()
3678 PetscCall(PetscArraycpy(S, Sarray + cumarray, subset_size * subset_size)); in PCBDDCAdaptiveSelection()
3679 PetscCall(PetscArraycpy(St, Starray + cumarray, subset_size * subset_size)); in PCBDDCAdaptiveSelection()
3686 PetscCall(PetscLogFlops((4.0 * subset_size * subset_size * subset_size) / 3.0)); in PCBDDCAdaptiveSelection()
3704 PetscCall(PetscLogFlops((4.0 * subset_size * subset_size * subset_size) / 3.0)); in PCBDDCAdaptiveSelection()
3710 PetscCall(PetscArraycpy(S, Sarray + cumarray, subset_size * subset_size)); in PCBDDCAdaptiveSelection()
3711 PetscCall(PetscArraycpy(St, Starray + cumarray, subset_size * subset_size)); in PCBDDCAdaptiveSelection()
3717 PetscCall(PetscLogFlops((4.0 * subset_size * subset_size * subset_size) / 3.0)); in PCBDDCAdaptiveSelection()
3729 PetscCall(PetscLogFlops((4.0 * subset_size * subset_size * subset_size) / 3.0)); in PCBDDCAdaptiveSelection()
3735 PetscCall(PetscArraycpy(S, Sarray + cumarray, subset_size * subset_size)); in PCBDDCAdaptiveSelection()
3736 PetscCall(PetscArraycpy(St, Starray + cumarray, subset_size * subset_size)); in PCBDDCAdaptiveSelection()
3742 PetscCall(PetscLogFlops((4.0 * subset_size * subset_size * subset_size) / 3.0)); in PCBDDCAdaptiveSelection()
3752 PetscCall(PetscLogFlops((4.0 * subset_size * subset_size * subset_size) / 3.0)); in PCBDDCAdaptiveSelection()
3777 PetscCall(PetscLogFlops((4.0 * subset_size * subset_size * subset_size) / 3.0)); in PCBDDCAdaptiveSelection()
3784 PetscCall(PetscArrayzero(eigv, subset_size * nmax)); in PCBDDCAdaptiveSelection()
3787 eigv[k * (subset_size + 1)] = 1.0; in PCBDDCAdaptiveSelection()
3824 for (j = 0; j < subset_size; j++) { in PCBDDCAdaptiveSelection()
3825 for (k = j; k < subset_size; k++) { in PCBDDCAdaptiveSelection()
3826 S[j * subset_size + k] = Sarray[cumarray + j * subset_size + k]; in PCBDDCAdaptiveSelection()
3827 St[j * subset_size + k] = Starray[cumarray + j * subset_size + k]; in PCBDDCAdaptiveSelection()
3831 PetscCall(PetscArraycpy(S, Sarray + cumarray, subset_size * subset_size)); in PCBDDCAdaptiveSelection()
3832 PetscCall(PetscArraycpy(St, Starray + cumarray, subset_size * subset_size)); in PCBDDCAdaptiveSelection()
3836 …upper, &B_IL, &B_IU, &eps, &B_neigs2, eigs + B_neigs, eigv + B_neigs * subset_size, &B_N, work, &B… in PCBDDCAdaptiveSelection()
3838 …upper, &B_IL, &B_IU, &eps, &B_neigs2, eigs + B_neigs, eigv + B_neigs * subset_size, &B_N, work, &B… in PCBDDCAdaptiveSelection()
3840 PetscCall(PetscLogFlops((4.0 * subset_size * subset_size * subset_size) / 3.0)); in PCBDDCAdaptiveSelection()
3880 PetscReal r = PetscRealPart(eigv[(ii + eigs_start) * subset_size + j]); in PCBDDCAdaptiveSelection()
3881 PetscReal c = PetscImaginaryPart(eigv[(ii + eigs_start) * subset_size + j]); in PCBDDCAdaptiveSelection()
3884 …Printf(pcbddc->dbg_viewer, " %1.4e\n", (double)(eigv[(ii + eigs_start) * subset_size + j]))); in PCBDDCAdaptiveSelection()
3890 …PetscCall(MatCreateSeqDense(PETSC_COMM_SELF, subset_size, B_neigs, eigv + eigs_start * subset_size… in PCBDDCAdaptiveSelection()
3899 …bddc->adaptive_constraints_data_ptr[cum], eigv + eigs_start * subset_size, B_neigs * subset_size)); in PCBDDCAdaptiveSelection()
3907 …PetscReal r = PetscRealPart(pcbddc->adaptive_constraints_data[ii * subset_size + j + pcbddc->adapt… in PCBDDCAdaptiveSelection()
3908 …PetscReal c = PetscImaginaryPart(pcbddc->adaptive_constraints_data[ii * subset_size + j + pcbddc->… in PCBDDCAdaptiveSelection()
3911 …%1.4e\n", (double)PetscRealPart(pcbddc->adaptive_constraints_data[ii * subset_size + j + pcbddc->a… in PCBDDCAdaptiveSelection()
3916 …cbddc->adaptive_constraints_idxs + pcbddc->adaptive_constraints_idxs_ptr[cum], idxs, subset_size)); in PCBDDCAdaptiveSelection()
3917 …>adaptive_constraints_idxs_ptr[cum + 1] = pcbddc->adaptive_constraints_idxs_ptr[cum] + subset_size; in PCBDDCAdaptiveSelection()
3918 …constraints_data_ptr[cum + 1] = pcbddc->adaptive_constraints_data_ptr[cum] + subset_size * B_neigs; in PCBDDCAdaptiveSelection()
3923 cumarray += subset_size * subset_size; in PCBDDCAdaptiveSelection()